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Presentación 


Señores miembros del jurado. 

Señor decano de la facultad de ingeniería. 

Señores docentes de la escuela de ingeniería mecánica. 

De conformidad a lo contemplado por la ley universitaria 30220, Art. 45. 
Inc. 45.2, en concordancia con lo dispuesto en el Art. 133 de los Estatutos de 
la Universidad Nacional de Trujillo, presento a vuestra disposición, bajo la 



El presente trabajo dqUnvestigacíón tiene como objetivo principal modelar 
la respuesta de cables de j^^frnpactados transversal mente por una 


masa puntual. Se ha logrado profundizar y ampliar el conocimiento sobre este 
fenómeno y, en consecuencia, se ha contribuido con herramientas analíticas 
para el desarrollo y diseño de barreras de protección viales, barreras de 
prevención de caída de rocas y bases para detener objetos que caen en 
minas ya que estos diseños sirven para disminuir la gran cantidad de 
accidentes mortales que pasan cada año en nuestro país y en todo el planeta. 

Tengo a bien poner a disposición este trabajo, a fin que pueda motivar la 
realización de futuras investigaciones en el tema del impacto transversal de 
cables. 
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Resumen 


Un mayor entendimiento del fenómeno del impacto transversal en 
cables de acero es necesario para poder diseñar más y mejores estructuras 
que eviten la gran cantidad de pérdida de vidas humanas en accidentes 
vehiculares, accidentes por caídas de objetos en minas y por la gran 
cantidad de deslizamientos de tierras y caídas de rocas omnipresentes en 
toda la geografía accidentada que tenemos. Por lo anterior, en esta tesis, se 
presentan los modelos matemáticos de las fuerzas en los anclajes de un 
cable de acero y el movimiento del cuerpo que lo impacta transversalmente 
en su punto medio. Para esfó se considera rígidos a los extremos del cable, 
idealizando su geometría como una cuerda unidimensional de material 
elástico e idealizando al cuerpo que choca como un punto material que 

iquinaria que pueda caer en una 
a una gran desaceleración lo que 
sfuerzos internos por contacto en el 


puede representar a 
barrera de protección sometí 



provoca grandes desplazamientos 
cable. Luego los modelos se^aNi^rfc.usando las ecuaciones desarrolladas 
para calcular la respuesta de los cables tipo 3x7 y 1x7 suponiendo las 
mismas condiciones que se impusieron en estudios previos sobre barreras 
de protección. Al final se concluye que los modelos, a excepción del modelo 
de los esfuerzos internos por contacto, si expresan la física de los 
experimentos y simulaciones referidos de manera que si se ha contribuido a 
aumentar nuestro conocimiento sobre el fenómeno del impacto transversal 
en cables de acero. Además, se ha aportado con un método para poder 
determinar las condiciones de frontera de una cuerda impactada 
transversalmente y que se puede usar como base para solucionar problemas 
más generales como el impacto oblicuo. 

Palabras clave: Impacto transversal, cables de acero, modelo 
matemático, barrera de protección. 
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Abstract 



ck, automobile or any machinery 
achieving this way a high level of 
displacements and internal stresses 


A deeper understanding in the transversal ¡mpact of wire rapes ¡s 
required to be able to design more and improved structures to reduce the 
soaring loss of human Ufes in traffic accident, by fall of object in mines and 
the huge amount of rockfalls and landslides ubiquitous all along the rugged 
geography that we have. By all these reasons, in this thesis are given the 
mathematical models of the torces in the ends of a Steel wire rape and the 
movement of the body that,p|jé^#felp§yersely in its middle point. To achieve 
these objectives, botlvéhds of the wire rape are considered rigidilly fixed, 
reducing its geometry to an unidimensional elastic string and modeling the 
body as a material poi 
that could fall over a protecto 
acceleration wath in turn induces 
by contact in the wire rape. Next the** models are validated using the 
developed equations to calcúlate the reaction of the Steel wire rapes type 3x7 
and 1x7, assuming equal conditions imposed in previous investigations done 
about protection barriers. At the end is concluded that all models, excepting 
the model of internal stresses, express the physics present in the refered 
experiments and simulations thus contributing to augment our knowledge 
about the transversal impact of Steel wire rapes. And additionally, it has been 
contributed with a method to determine the frontier conditions in a 
transversely impacted string that can also be used as basis to solve more 
general problems, as the oblique impact of wire rapes. 

Keywords: Transversal impact, Steel wire rape, mathematical model, 
protection barrier. 


v 


Esta obra ha sido publicada bajo la licencia Creative Commons Reconocimiento-No Comercial-Compartir bajo la misma licencia 2.5 Perú. 
Para ver una copia de dicha licencia, visite http://creativecommons.Org/licences/by-nc-sa/2.5/pe/ 







BIBLIOTECA DIGITAL - DIRECCIÓN DE SISTEMAS DE INFORMÁTICA Y COMUNICACIÓN 


índice general 


Dedicatoria.i 

Presentación.¡i 

Agradecimientos.iii 

Resumen.¡v 

Abstract.v 

índice de tablas.x 

índice de figuras.xiii 


índice de figuras.xiii 



1.4.2 . Justificación tecnológica.7 

1.4.3 . Justificación económico-social.8 

1.5 . Objetivos.8 

1.5.1 . Objetivo general.8 

1.5.2 . Objetivos específicos.8 

1.6 . Alcance de la tesis.9 

1.7 . Organización de la tesis.10 

2 . Marco teórico.11 

2.1 . ¿Cuándo es relevante el impacto ?.11 


vi 


Esta obra ha sido publicada bajo la licencia Creative Commons Reconocimiento-No Comercial-Compartir bajo la misma licencia 2.5 Perú. 
Para ver una copia de dicha licencia, visite http://creativecommons.Org/licences/by-nc-sa/2.5/pe/ 
































BIBLIOTECA DIGITAL - DIRECCIÓN DE SISTEMAS DE INFORMÁTICA Y COMUNICACIÓN 


2.2 . Desplazamientos y deformaciones en cuerdas.13 

2.3 . Ondas de esfuerzo unidimensional.15 

2.4 . Ecuación general del movimiento de una cuerda elástica.16 

2.5 . Módulo de elasticidad de los cables trenzados de acero.19 

2.6 . Física del impacto transversal en cuerdas elásticas.22 

3 . Materiales y métodos.26 

3.1 . Recolección de información.26 

3.2 . Material de estudio.26 


3.3 . Metodología a emplear p^^gj^lanteamiento de los modelos 
matemáticos. 


3.4 . Variables de investig 

3.5 . Procedimiento 



.27 

.28 

.30 


3.5.1 . Restricciones, simplificacior^^ébmétricas y físicas del fenómeno.. 30 

3.5.2 . Procedimiento para ra/%i¡i^| 3 §tdm de los modelos matemáticos.32 

3.5.3 . Procedimiento para la validación de los modelos matemáticos.34 

4 . Planteamiento de los modelos matemáticos.35 

4.1 . Planteamiento del modelo matemático de los esfuerzos axiales.35 

4.2 . Planteamiento del modelo matemático del movimiento del cuerpo que 

impacta transversalmente un cable en su punto medio.57 

4.3 . Planteamiento del modelo matemático de los esfuerzos internos debido 

al contacto.60 

5 . Validación de los modelos matemáticos.65 

5.1 . Validación de los modelos de los esfuerzos axiales y movimiento del 
cuerpo que impacta con los experimentos de C. Stolle.65 

vii 


Esta obra ha sido publicada bajo la licencia Creative Commons Reconocimiento-No Comercial-Compartir bajo la misma licencia 2.5 Perú. 
Para ver una copia de dicha licencia, visite http://creativecommons.Org/licences/by-nc-sa/2.5/pe/ 






























BIBLIOTECA DIGITAL - DIRECCIÓN DE SISTEMAS DE INFORMÁTICA Y COMUNICACIÓN 


5.1.1 . Modelamiento del experimento 1 de C. Stolle. 

5.1.1.1 . Discusión del modelamiento del experimento 1 de C. Stolle. 

5.1.2 . Modelamiento del experimento 2 de C. Stolle. 

5.1.2.1 . Discusión del modelamiento del experimento 2 de C. Stolle. 

5.1.3 . Modelamiento del experimento 4 de C. Stolle. 

5.1.3.1 . Discusión del modelamiento del experimento 4 de C. Stolle. 

5.2 . Validación de los modelos de los esfuerzos axiales y movimiento del 
cuerpo que impacta con las simulaciones de Tytko et al y P. 
Gospodarczyk. 


5.2.1 . Modelamiento de la áimuladóntfé Tytko et al. 


5.2.1.1 . Discusión delpioc 
5.2.2 . Modelamiento de 



simulación de Tytko et al. 


)spodarczyk. 


5.2.2.1 . Discusión del modelamien^de la simulación de Gospodarczyk... 
5.3 . Comparación cualitativa de la aplicación del modelo de los esfuerzos 


internos debido al contacto. 


..68 

..73 

..76 

..81 

..83 

..88 

..89 

..91 

..93 

..94 

..99 

101 


5.3.1 . Aplicación del modelo de esfuerzos internos por contacto al 

experimento 1 de C. Stolle.102 

5.3.2 . Aplicación del modelo de esfuerzos internos por contacto al 

experimento 2 de C. Stolle.105 

5.3.2.1 . Discusión de la aplicación del modelo de esfuerzos internos debido 

al contacto a los experimentos 1 y 2 de C. Stolle.107 

5.3.3 . Aplicación del modelo de esfuerzos internos por contacto al 

experimento 4 de C. Stolle.108 

5.3.3.1 . Discusión de la aplicación del modelo de esfuerzos internos debido 

al contacto al experimento 4 de C. Stolle.111 


viii 


Esta obra ha sido publicada bajo la licencia Creative Commons Reconocimiento-No Comercial-Compartir bajo la misma licencia 2.5 Perú. 
Para ver una copia de dicha licencia, visite http://creativecommons.Org/licences/by-nc-sa/2.5/pe/ 
























BIBLIOTECA DIGITAL - DIRECCIÓN DE SISTEMAS DE INFORMÁTICA Y COMUNICACIÓN 


5.3.4 . Aplicación del modelo de esfuerzos internos debido al contacto a la 
simulación de Gospodarczyk.113 

5.3.4.1 . Discusión de la aplicación del modelo de esfuerzos internos debido 

al contacto a la simulación de Gospodarczyk.115 

6 . Conclusiones y recomendaciones.116 

6.1 . Conclusiones sobre la validación de los modelos de los esfuerzos 

axiales y movimiento del cuerpo que impacta.117 

6.2 . Conclusiones sobre la aplicación del modelo de los esfuerzos internos 

debido al contacto.117 



6.3 . Resumen de resultados. 

6.4 . Sugerencias 

7 . Bibliografía. 

8 . Anexos. 

8.1 . Anexo 1. 


119 

120 
122 
131 
131 


ix 


Esta obra ha sido publicada bajo la licencia Creative Commons Reconocimiento-No Comercial-Compartir bajo la misma licencia 2.5 Perú. 
Para ver una copia de dicha licencia, visite http://creativecommons.Org/licences/by-nc-sa/2.5/pe/ 





















BIBLIOTECA DIGITAL - DIRECCIÓN DE SISTEMAS DE INFORMÁTICA Y COMUNICACIÓN 


índice de tablas 

Tabla 1: Comparación de masas del cable con el cuerpo que lo impacta.37 

Tabla 2: Tiempo de reflexión de las ondas longitudinales en los cables 

estudiados.38 


Tabla 

Tabla 

Tabla 

Tabla 

Tabla 

Tabla 

Tabla 

Tabla 

Tabla 

Tabla 

Tabla 

Tabla 

Tabla 

Tabla 


3: Características físicas del cable de acero tipo 3x7 de los 

experimentos de C. Stolle [31].65 

4: Variables de entrada para el modelamiento matemático del 

experimento 1 de C. Stolle [31].68 


5: Parámetros del modelo matemático del experimento 1 de C. Stolle 
durante el proceso de carga.68 

6: Parámetros del modelo matemático del experimento 1 de C. Stolle 
durante el proceso de descarga.68 

7: Comparación dej^cálculos defsKjensión en los anclajes del cable 
¡mpactado del Exp.1 de C. Stolle.69 

8: Comparación de jfap cálculos del desplazamiento del cuerpo que 
impacta del Exp.1 de C. .70 

9: Comparación de los cálculos de la velocidad del cuerpo que impacta 
del Exp.1 de C. Stolle. wÁ .71 

10: Comparación deláte^l^l^^Efé'la aceleración del cuerpo que 

impacta del Exp.1 de C. Stolle.72 


11: Variables de entrada para el modelamiento matemático del 

experimento 2 de C. Stolle [31].76 

12: Parámetros del modelo matemático del experimento 2 de C. Stolle 
durante el proceso de carga.76 


13: Parámetros del modelo matemático del experimento 2 de C. Stolle 
durante el proceso de descarga.76 

14: Comparación de los cálculos de la tensión en los anclajes del cable 
impactado del Exp.2 de C. Stolle.77 


15: Comparación de los cálculos del desplazamiento del cuerpo que 


impacta del Exp.2 de C. Stolle.78 

16: Comparación de los cálculos de la velocidad del cuerpo que 

impacta del Exp.2 de C. Stolle.79 


x 


Esta obra ha sido publicada bajo la licencia Creative Commons Reconocimiento-No Comercial-Compartir bajo la misma licencia 2.5 Perú. 
Para ver una copia de dicha licencia, visite http://creativecommons.Org/licences/by-nc-sa/2.5/pe/ 





















BIBLIOTECA DIGITAL - DIRECCIÓN DE SISTEMAS DE INFORMÁTICA Y COMUNICACIÓN 


Tabla 

Tabla 

Tabla 

Tabla 

Tabla 

Tabla 

Tabla 

Tabla 

Tabla 

Tabla 

Tabla 

Tabla 

Tabla 

Tabla 

Tabla 


17: 

18: 

19: 

20 : 

21 : 

22 : 

23: 

24: 

25: 

26: 

27: 

28: 

29: 

30: 

31: 


Comparación de los cálculos de la aceleración del cuerpo que 
impacta del Exp.2 de C. Stolle.80 

Variables de entrada para el modelamiento matemático del 
experimento 4 de C. Stolle [31].83 

Parámetros del modelo matemático del experimento 4 de C. Stolle 
durante el proceso de carga.83 

Comparación de los cálculos de la tensión en los anclajes del cable 
impactado del Exp.4 de C. Stolle.84 


Comparación de los cálculos del desplazamiento del cuerpo que 
impacta del Exp.4 de C. Stolle.85 


Comparación de los cálculos de la velocidad del cuerpo que 
impacta del Exp.4 de C. Stolle. 

ace | erac ¡ón del cuerpo que 


.86 


lie. 



.87 


e áfcero tipo 1x7 de las 
Gospodarczyk [33].89 

vrriddelamiento matemático de la 


Comparación de 
impacta del Ex#í de 

Características fí; 
simulaciones üwtmm 

Variables de entrada 
simulación de Tytko et al [2910.91 

Parámetros del modelo matemático de la simulación de Tytko et al 
durante el proceso de carga hasta su rotura.91 

Comparación de los cálculos de la tensión en el cable de la 
simulación de Tytko et al.92 

Variables de entrada para el modelamiento matemático de la 
simulación de P. Gospodarczyk [33].94 


Parámetros del modelo matemático de la simulación de 
Gospodarczyk durante el proceso de carga de los cables de la base 
de seguridad hasta su rotura.96 

Comparación del cálculo de la velocidad final del bloque de 18 TN 
después de pasar por la base de seguridad simulada por 
Gospodarczyk.98 

Variables de entrada para el modelo de esfuerzos internos debido al 
contacto del Exp.1 de C. Stolle.103 


x¡ 


Esta obra ha sido publicada bajo la licencia Creative Commons Reconocimiento-No Comercial-Compartir bajo la misma licencia 2.5 Perú. 
Para ver una copia de dicha licencia, visite http://creativecommons.Org/licences/by-nc-sa/2.5/pe/ 





















BIBLIOTECA DIGITAL - DIRECCIÓN DE SISTEMAS DE INFORMÁTICA Y COMUNICACIÓN 


Tabla 32: Resultado de los cálculos con los datos del Exp. 1 de C. Stolle 
usando el modelo del contacto entre dos cilindricos 
perpendiculares.103 

Tabla 33: Variables de entrada para el modelo de esfuerzos internos debido al 
contacto del Exp.2 de C. Stolle.105 

Tabla 34: Resultado de los cálculos con los datos del Exp. 2 de C. Stolle 

usando el modelo del contacto entre dos cilindricos perpendiculares. 
.105 


Tabla 35: Variables de entrada para el modelo de esfuerzos internos debido al 
contacto del Exp.4 de C. Stolle.108 

Tabla 36: Resultado de los cálculos con los datos del Exp. 4 de C. Stolle 

usando el modelo del contacto entre dos cilindricos perpendiculares. 
.109 

Tabla 37: Resultado de los c os datos del Exp. 4 de C. Stolle 

usando el modelo < entre dos cilindricos perpendiculares 

con dos tensiones ¡^nativas.109 


rczyk para el modelo de 
acto.113 


Tabla 38: Datos de la si 

esfuerzos internos d 

Tabla 39: Cálculos con los datos d 



iulación de Gospodarczyk usando el 
modelo de esfuerzos internos para los cables bajo el cuerpo que 
i m pacta.. ..113 


xii 


Esta obra ha sido publicada bajo la licencia Creative Commons Reconocimiento-No Comercial-Compartir bajo la misma licencia 2.5 Perú. 
Para ver una copia de dicha licencia, visite http://creativecommons.Org/licences/by-nc-sa/2.5/pe/ 














BIBLIOTECA DIGITAL - DIRECCIÓN DE SISTEMAS DE INFORMÁTICA Y COMUNICACIÓN 


índice de figuras 

Figura 1: Ejemplo de una barrera vial de cables de acero [19].2 

Figura 2: Ilustración del funcionamiento de una barrera vial de cables 
CASS™ redireccionando al vehículo que impacta [20].2 

Figura 3: Barrera dinámica usada en la quebrada Carosio-Chosica 
conteniendo detritos de un deslizamiento en el verano de 2017 
detritos [21].2 


Figura 4: Ejemplo de una barrera deteniendo una roca que puede impactar 
autos si cae hacia una carretera [22].2 

Figura 5: Diagrama del cable impactado asumiendo un perfil triangular. , 
velocidad constajttíé' del punto' impactado; , velocidad de la onda 


transversal, ,rífelo# 

/<§5 

/ ^ 

Figura 6: Equilibrio diná 
variable., es 
actual; un elemento in 



da longitudinal [34].5 

r\ 

eformable con sección transversal 
referencia, , es la coordenada 
imal de longitud , secciones y , 


desplazamiento cambio infinitesimal de su área , con un esfuerzo 

. 


uniformemente distriBÉÍÉJtry la velocidad del frente de la onda de 
esfuerzos o del sonido [50], [51].15 


Figura 7: Segmento infinitesimal de una cuerda en el espacio sometida a 


fuerzas internas y externas [37].16 

Figura 8: Elementos de un cable de acero [54].19 


Figura 9: Definiciones del módulo de elasticidad como pendientes de las 
secantes OB, CB, y EF o de las tangentes en O o A. , es el 
esfuerzo, , es la deformación final, , es la deformación durante la 
carga y , es la deformación de descarga; , son los ángulos de las 
secantes o pendientes. [56].20 


xiii 


Esta obra ha sido publicada bajo la licencia Creative Commons Reconocimiento-No Comercial-Compartir bajo la misma licencia 2.5 Perú. 
Para ver una copia de dicha licencia, visite http://creativecommons.Org/licences/by-nc-sa/2.5/pe/ 
















BIBLIOTECA DIGITAL - DIRECCIÓN DE SISTEMAS DE INFORMÁTICA Y COMUNICACIÓN 


Figura 10: Curvas esfuerzo , deformación , de varios tipos de cables de acero, 
nuevos, ante una primera carga y descarga. [31], [55], [57], [58], 
[59], [60], [61].21 

Figura 11: Movimiento de una cuerda impactada transversalmente por una 
masa puntual. M, es la masa del cuerpo, V00, es su velocidad 
inicial, V, es su velocidad actual, O, es el origen de coordenadas de 
ejes x, y, S, es la coordenada curvilínea actual de la cuerda y SO, 
es la coordenada de referencia.25 


Figura 12: Representación como caja gris de los modelos matemáticos 
planteados.30 


Figura 13: Representación 
¡mpactado pofá^a 
es su velocidad 
coordenadas 
la cuerda y 



rica de la idealización de un cable 
puntual. M, es la masa del cuerpo, V00, 
velocidad actual, O, es el origen de 
a coordenada curvilínea actual de 
e referencia. L, es la longitud del 


cable.lo.32 

Figura 14: Modelo de cuétel^ delimjl^tes de la cuerda ideal sometida a 
impacto transversal en su punto medio. M, es la masa del cuerpo, 
V00, es su velocidad inicial, V, es su velocidad actual, O, es el 
origen de coordenadas de ejes x, y, S, es la coordenada curvilínea 
actual de la cuerda y SO, es la coordenada de referencia y 
extremos . YR, desplazamiento del punto medio de la cuerda 
celeste. YM, desplazamiento del punto medio de la cuerda verde o 
del punto material. L, es la longitud del cable.39 


Figura 15 Analogía del modelo delimitante de la tensión en los extremos con 
un modelo de resortes ideales. T, es la tensión en los resortes, a, 
es el ángulo de inclinación de los resortes, F, es la fuerza de 
reacción sobre el cuerpo que choca, M, es la masa del cuerpo, YR, 
es su desplazamiento, x, y son las coordenadas y L, la distancia 
entre ejes.40 


xiv 


Esta obra ha sido publicada bajo la licencia Creative Commons Reconocimiento-No Comercial-Compartir bajo la misma licencia 2.5 Perú. 
Para ver una copia de dicha licencia, visite http://creativecommons.Org/licences/by-nc-sa/2.5/pe/ 













BIBLIOTECA DIGITAL - DIRECCIÓN DE SISTEMAS DE INFORMÁTICA Y COMUNICACIÓN 


Figura 16: Ciclo de carga y descarga de un cable perfectamente elástico, o, 
es el esfuerzo y e, es la deformación .41 

Figura 17: Cambio de la tensión en un cable idealizado como una cuerda 
impactada en su punto medio que toma una forma triangular. M, es 
la masa del cuerpo que choca, a, es el ángulo de inclinación de los 
resortes, F, es la fuerza de reacción sobre el cuerpo que choca, T, 
es la tensión distribuida en la cuerda, C11, es la velocidad del frente 
de la onda de esfuerzos longitudinal. El cable está idealizado como 
las cuerdas de “movimientos delimitantes”, sin masa.42 


Figura 18: Idealización lineal de los procesos de carga y descarga de un 
cable de acero, o, es el esfuerzo, fMáx, es la deformación máxima 
y sf, es la deforpift&ibn final.'^^ .49 


Figura 19: Ilustración ^ara 
procesos de 
elásticas e ideática 

mínimos) e impactadas lp| 

w 



simetría respecto al tiempo de los 
de dos cuerdas perfectamente 
e tensión T, iguales (máximos y 
inalmente por una masa puntual M, 


ti, son instantes de tiempo .Igualmente espaciados, tn, es el 
periodo, g, es la grávéHatl, x,y, son las coordenadas del espacio, S, 
es la coordenada curvilínea.50 


Figura 20: Representación de la Tensión Vs tiempo de los procesos de carga 
y descarga de las cuerdas y correspondiente a lo expresado en la 
figura 19. Tmáx, tensión máxima durante la carga o inicial de la 
descarga, Tmín, tensión mínima al iniciar la carga o finalizar la 
descarga, TD, tensión de descarga, TC.Sim., tensión de descarga 


como función simétrica de la carga.52 

Figura 21: Cilindros de diámetros DI, y D2 con sus ejes a 90°. x,y,z, son los 
ejes coordenados con origen en el punto de contacto.62 

Figura 22: Cilindro de radio R y plano rígido de longitud ñ en contacto.x,y,z, 
son los ejes coordenados con origen en el punto de contacto.63 
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Figura 23: Estado de esfuerzos principales máximos, a lo largo del eje z en el 
interior de un cable transversalmente impactado.Fmáx.es la fuerza 
de reacción que ejerce el cuerpo que choca, Tmáx.es la tensión 
máxima que alcanza el cable en los anclajes, oxx es el esfuerzo 
princiapal en la dirección x, oyy es el esfuerzo princiapal en la 
dirección y, ozz es el esfuerzo princiapal en la dirección z, o es el 
esfuerzo debido a Tmáx.64 


Figura 24: Detalle del anclaje del cable 3x7 usado para los experimentos de 
impacto transversal de C. Stolle [31].66 

Figura 25: Sección transversal del cable 3x7 usado en los experimentos de C. 



Figura 26: Posición de‘impacj^ó entre el carro y el cable 3x7 en los 

'1 Ov 

experimento^ de cables sometidos a impacto transversal de C. 


Stolle [31] 


Figura 27: Idealización de lo 



.66 


tos 1,2 y 4 de C. Stolle. L, es la 
le, MC, es la masa del cuerpo que 
impacta, V00, es ^.^elocid^^Ól cuerpo en el impacto, x, es el eje 


longitud entre anclajes 


coordenado y O, es el punto medio del cable y origen de 
coordenadas.67 


Figura 28: Comparación del Modelo Matemático de la tensión T(KN) en los 
anclajes del cable impactado transversalmente en el experimento 1 
de C. Stolle. Las tensiones indicadas con flechas son las tensiones 
máximas.69 


Figura 29: Comparación del modelo matemático del desplazamiento del 


cuerpo M que impacta en el experimento 1 de C. Stolle.70 

Figura 30: Comparación del modelo matemático de la velocidad del cuerpo M 
que impacta en el experimento 1 de C. Stolle.71 
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Figura 31: Comparación del modelo matemático de la aceleración 
(normalizada con la gravedad 9.81 m/S2) del cuerpo M que 
impacta en el experimento 1 de C. Stolle.72 

Figura 32: Comparación del Modelo Matemático de la tensión T(KN) en los 
anclajes del cable impactado transversalmente en el experimento 2 
de C. Stolle.Las tensiones indicadas con flechas son las tensiones 
máximas.77 


Figura 33: Comparación del modelo matemático del desplazamiento del 
cuerpo M que impacta en el experimento 2 de C. Stolle.78 

Figura 34: Comparación del modelo matemático de la velocidad del cuerpo M 
que impacta en el experimento 2 de C. Stolle.79 

Figura 35: Compararon de^npdelo matemático de la aceleración 
(normalizada o ^ ' 


impacta en el 


Figura 36: Deslizamiento final de 


9r81 m/S2) del cuerpo M que 
. Stolle.80 

del cable 3x7 en la resina de los 


sockets usados como anclajee^en el experimento 2 de C. Stolle 
[31]..82 


Figura 37: Comparación del Modelo Matemático de la tensión T(KN) en los 
anclajes del cable impactado transversalmente en el experimento 4 
de C. Stolle. La tensión indicada con una flecha es la tensión 
máxima.84 


Figura 38: Comparación del modelo matemático del desplazamiento del 


cuerpo M que impacta en el experimento 4 de C. Stolle.85 

Figura 39: Comparación del modelo matemático de la velocidad del cuerpo M 
que impacta en el experimento 4 de C. Stolle.86 


Figura 40: Comparación del modelo matemático de la aceleración 
(normalizada con la gravedad 9.81 m/S2) del cuerpo M que 
impacta en el experimento 4 de C. Stolle.87 
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Figura 41: Sección transversal del cable 1x7.89 

Figura 42: Idealización de la simulación presentada por Tytko et al.L, es la 
longitud entre anclajes del cable, MC, es la masa del cuerpo que 
impacta, V00, es la velocidad del cuerpo en el impacto, x, es el eje 
coordenado y O, es el punto medio del cable y origen de 
coordenadas.90 


Figura 43: Resultado de la simulación de Tytko et al mostrando la ruptura del 
cable 1x7 de 1m.90 


Figura 44: Comparación del Modelo Matemático de la tensión T(KN) en los 
anclajes del cable impactado. transversal mente de la simulación de 
Tytko et al. 246 KN ^^y^ite elástico estándar del cable, 279 KN 
es el límite dé resistencia estándar y 289.63 KN es la tensión 

.92 


máxima o de roti 


Figura 45: Posición inici 

base artificial de 0 

bloque que cae, 2 son k 

% 

. 



echa por P. Gospodarczyk de una 
minas de pozo profundo. 1 es el 
cables y 3 una placa de acero [33], 
.94 


Figura 46: Primera idealización de la simulación de una base de seguridad de 
minas de pozo profundo presentada por Gospodarczyk.95 


Figura 47: Segunda idealización de la simulación presentada por 
Gospodarczyk de una base de seguridad de minas de pozo 
profundo.L, es la longitud entre anclajes del cable, MC, es la masa 
del cuerpo que impacta, V00, es la velocidad del cuerpo en el 
impacto, x, es el eje coordenado y O, es el punto medio del cable y 
origen de coordenadas.95 


Figura 48: Modelo Matemático de la tensión T(KN) en los anclajes de los 
cables impactados transversalmente en una base de protección 
para minas simulada por Gospodarczyk. La tensión indicada con 
una flecha es la tensión máxima.97 
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Figura 49: Comparación del modelo matemático de la velocidad y la 
velocidad final simulada del bloque de 18 TN a través de la base 
de seguridad de Gospodarczyk. La velocidad indicada es la que 
alcanza el bloque luego de traspasar la plataforma de cables 
simulada y que se da al romper el cable a una tensión de 279 KN. 
.98 


Figura 50: Situación final de la plataforma simulada por Gospodarczyk, con el 
bloque desacelerando desde 75m/s hasta 48 m/s [33].100 

Figura 51: Frente de impacto del carro usado por C. Stolle en sus 
experimentos [31].102 

Figura 52: Modelo matemáticcxclgirl^s esfuerzos internos en el cable del Exp. 

1 de C. Stolle^&as temples señaladas son las máximas calculada 

. Máx es el esfuerzo de rotura del 

de rotura nominal y su área 

.104 


y medida en^el 
cable calcula 
metálica.... 



Figura 53: Modelo matemático de 

% 


uerzos internos en el cable del Exp. 
2 de C. Stolle.tensiones señaladas son las máximas calculada 
y medida en el experimento, oR. Máx es el esfuerzo de rotura del 


cable calculado con la tensión de rotura nominal y su área 
metálica.106 


Figura 54: Deformación del cable en la zona del impacto al final del Exp. 2 de 
C. Stolle [31].107 

Figura 55: Deformación del cable en la zona del impacto al final del Exp. 1 de 
C. Stolle [31].107 

Figura 56: Modelo matemático de los esfuerzos internos en el cable del Exp. 

4 de C. Stolle. Las tensiones señaladas son las máximas 
calculada, con dos alternativas, y la tensión máxima medida en el 
experimento, oR. Máx es el esfuerzo de rotura del cable calculado 
con la tensión de rotura nominal y su área metálica.110 


xix 


Esta obra ha sido publicada bajo la licencia Creative Commons Reconocimiento-No Comercial-Compartir bajo la misma licencia 2.5 Perú. 
Para ver una copia de dicha licencia, visite http://creativecommons.Org/licences/by-nc-sa/2.5/pe/ 














BIBLIOTECA DIGITAL - DIRECCIÓN DE SISTEMAS DE INFORMÁTICA Y COMUNICACIÓN 


Figura 57: Situación final del cable 3x7 en el frente del carro que lo impacto 
transversalmente en el experimento 4 de C. Stolle [31].112 

Figura 58: Modelo matemático de los esfuerzos internos en los cables de la 
base de seguridad para minas simulada por Gospodarczyk. Las 
tensiones señaladas son las máximas calculadas para las cuales el 
cable se rompe totalmente según el modelo, oR. Máx es el 
esfuerzo de rotura del cable calculado con la tensión de rotura 
nominal y su área metálica.114 


Figura 59: Modelo del impacto transversal de un cable aproximado con un 
perfil triangular y desde el instante de contacto.M, es la masa del 


cuerpo, V00, es su velocidad inicial, V, es su velocidad actual, O, 
es el origen de, tarden acf¿fcc$e ejes x, y, Yo, desplazamiento del 

Jas. L, es la longitud del cable, F, es la 

r¡'ü% m 

ílocidad del frente de la onda de 


punto medio I 


fuerza de reao 


esfuerzos lo 


onda transversal 



j^tancia recorrida por el frente de 
.131 


Figura 60: Comparación:#,el Modelo Matemático de la tensión T(KN) en los 

'(// - 




anclajes del cable imp&GÍá'd'ó transversalmente en el experimento 1 
de C. Stolle con la simulación de C. Stolle y otros modelos 
numéricos que usó para comparar su simulación [31].136 


xx 


Esta obra ha sido publicada bajo la licencia Creative Commons Reconocimiento-No Comercial-Compartir bajo la misma licencia 2.5 Perú. 
Para ver una copia de dicha licencia, visite http://creativecommons.Org/licences/by-nc-sa/2.5/pe/ 











BIBLIOTECA DIGITAL - DIRECCIÓN DE SISTEMAS DE INFORMÁTICA Y COMUNICACIÓN 


Lista de símbolos 


A 

Área transversal de una barra o área combinada de las 

secciones transversales de los alambres de un cable. 

a 

Radio del área plana que se forma en el contacto de los 

cilindros con ejes perpendiculares. 

b 

Ancho de la franja rectangular plana que se forma en el cilindro 

en contacto con un sólido plano. 

B, B 1 B 2 ... 

Constantes adimensionales. 

c 

Velocidad de la onda de esfuerzo longitudinal o velocidad del 

Cj 

sonido en un mqterfáfAM / 0 

Velocidacy&e propagación deltas ondas transversales en una 

c n 

cuerda deforma^ | E^j^^^onal. 

Velocidad de propagación de las ondas longitudinales en una 

Di y D 2 

cuerda deformable ur^jn§nsional. 

Son los diámetrogj^ jó§‘©mñdros en contacto. 

E 

Módulo elástico 

e d 

Módulo elástico durante la descarga de un cable impactado 

transversalmente. 


E s (cr /ower , a upper ) Módulo de elasticidad de un cable de acero medido 
como la pendiente de la secante entre un esfuerzo inicial y 


E s ( &upper) 

otro final. 

Módulo de elasticidad de un cable de acero medido como la 

pendiente de la secante entre un esfuerzo inicial cero y 

final. 
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E t,up{o) 

Módulo de elasticidad de un cable de acero medido como la 

pendiente de la tangente en cada punto de la curva de carga 

(aumento de la tensión) desde una tensión inicial. 

Et,down ( 

Módulo de elasticidad de un cable de acero medido como la 

pendiente de la tangente en cada punto de la curva de 

descarga (decremento de la tensión) desde un esfuerzo 

máximo alcanzado en un proceso de carga anterior. 

f() 

Función. 

F 

Fuerza de contacto o de reacción entre el cuerpo que impacta y 

el cable impactado. 

Fe 

Fuerza de reacción durarfte la carga entre el cable y el cuerpo 

que impacta. 

f d 

Fuerza de reacción durantes la descarga entre el cable y el 

cuerpo que impaHs¿^jE ; ^ 

ff(0 

Función destiempo. 

: Vru „ - jpF 

¡rNé 

H 1 y H 2 

Parámetros de curvatura para el cálculo del radio equivalente 

de dos cilindros en contacto. 

i 

Número natural diferente de cero. 

m 

Función del tiempo. 

kv 

Velocidad de la onda transveral o del acodado en un perfil 

triangular. 

L 

Longitud del cable o distancia característica de un sólido. 

m 

Coeficiente constante. 

mm' , nn' 

Secciones transversales de una barra 

M c 

Masa del cuerpo que impacta. 
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M 


Cable • 


Masa del cable. 


M 


Masa combinada del cuerpo que impacta y el cable. 


n 


Número natural diferente de cero. 


M 


n 


Conjunto de los números naturales. 

Longitud de contacto del frente plano de un sólido en contacto 
con un cable. 


O 

Po 

Re 

R-i, R-2, 


Origen del sistema coordenado plano. 

Presión de contacto máxima entre dos cuerpos sólidos 
deformables en contacto. 

Radio equivalentejdt^ los cilindros en contacto con sus ejes 
perpendicular! 

/ £ *^l 


extremos 1 y 2 de la cuerda 


Notación 
adimensional 



R, R 3 R 4 R 78 etc Constantes adimensionales. 


Coordenada curvilínea de referencia o Lagrangiana. 


s,S 


Coordenada curvilínea de una cuerda unidimensional. 


At 


Intervalo de tiempo e igual a Ar = r n /n 


Ar 


Intervalo de tiempo finito. 


t l t 2 t 3 ...t¡ Instantes de tiempo espaciados por el intervalo de tiempo Ar 
T A {ti),T B (ti) Tensiones en diferentes instantes de tiempo. 


C.Máx 


D.Máx 


Tensión de carga máxima en los extremos de la cuerda. 

Tensión de carga máxima en los extremos de una cuerda en un 
proceso de descarga simétrico. 
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T 
1 c 

T 

1 D 


T 

||T 


Tensión durante la carga en los extremos del cable calculada 
con el modelo matemático de desplazamientos pequeños. 

Tensión durante la descarga en los extremos del cable 
calculada con el modelo matemático de desplazamientos 
pequeños. 

Tensión en los puntos de una cuerda unidimensional. 

Módulo del vector tensión de una cuerda unidimensional. 


t-C.MÁX 


t 


D.MÁX 


t 


n 


Ír 


t 


total 


t 


c 


Tiempo total del proceso de carga o tiempo en que le toma al 
cable alcanzar una tensión máxima desde el instante del 
impacto. 

¿fe. 

;o de descarga o tiempo en que le toma 
fciónlmáxima desde el instante del 
tao de carga. 

Tiempo de carga' i &^ife de una cuerda ideal sin rozamiento 
interno y perfectamen 

Tiempo de reflexión o tiempo de ida y vuelta del frente de onda 
longitudinal desde el punto de impacto. 


Tiempo total del 

w 

al cablejálc 
impacto e 



Tiempo total de desaceleración y aceleración de vuelta del 
cuerpo que impacta o tiempo de carga y descarga del cable 
entre su tensión inicial mínima y tensión final cuando el cuerpo 
que impacta regresa y se desprende del cable o lo rompe. 

Tiempo característico de una colisión hasta que la fuerza 
alcanza valor máximo o vuelve a uno mínimo. 


t Tiempo. 

u Desplazamientos en dirección x. 
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u h u 2 

Parámetros característicos del movimiento de la cuerda 

delimitante con desplazamientos pequeños. 

V* 

Velocidad inicial del punto material o en el instante que impacta. 

V 

Velocidad de la partícula material durante el impacto 


Velocidad durante la descarga del punto medio de la cuerda con 

desplazamientos pequeños. 

w(q) 

Función de la variable q . 

x , y , z 

Coordenadas de referencia. 

X 

Coordenada actual o fuerza por unidad de longitud en la 

( X,Y,Z) 

dirección x d£%ra cgerda. 

Fuerzas externas en lo^pWto^e una cuerda. 

y 0 

Ordenadcy^^i meclio 0 P unt0 de impacto de la 

Y R. Corr. D 

Desplazamiento durante\j$r descarga de la cuerda con 
desplazamientos fjé^iM’ños. 

y R.Corr.S(ED) 

Desplazamiento durante la descarga expresado como función 

simétrica del desplazamiento durante la carga calculado con el 

módulo de elasticidad E D de descarga de la cuerda y con el 

modelo de desplazamientos pequeños 


yR.corr.s(ED.Máx) Máximo el desplazamiento durante la descarga expresado 
como función simétrica del desplazamiento durante la carga 
calculado con el módulo de elasticidad E D de descarga de la 
cuerda y con el modelo de desplazamientos pequeños. 
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y R. aprox. 

Desplazamiento, solución aproximada, de la ecuación de 

movimiento del punto material que impacta en la cuerda 

delimitante de desplazamientos pequeños. 

y R. aprox. 

valor corregido del desplazamiento, solución aproximada, de la 

ecuación de movimiento del punto material que impacta en la 

cuerda delimitante de desplazamientos pequeños. 

y M 

Desplazamiento de la partícula material (modelo de 

desplazamientos grandes). 

y m.c 

Desplazamiento de la partícula material durante la carga 

y m.d 

(modelo de desplazamientos grandes). 

Desplazami^fwO de la partícula material durante la descarga 

(modelo de desplazamientos grandes). 

y" M 

Aceleración ||e la wamtícula material (modelo de 
desplazamientos^randes). 

y"m.c 

Aceleracióq;de la partfcula’m| l terial durante la carga (modelo de 
desplazamiento^ grandes). 

y "m.d 

Aceleración de la partícula material durante la descarga 

(modelo de desplazamientos grandes). 

y ’m 

Velocidad de la partícula material (modelo de desplazamientos 

grandes). 

y m.c 

Velocidad de la partícula material durante la carga (modelo de 

desplazamientos grandes). 

y m.d 

Velocidad de la partícula material durante la descarga (modelo 

de desplazamientos grandes). 

y R 

Desplazamiento (modelo de pequeños desplazamientos) del 

punto medio de la cuerda delimitante. 
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( dyldx ) 2 

Cuadrado de la pendiente del eje neutro de una barra. 

a 

Ángulo que forma el cable con su posición inicial al suponerse 
que toma un perfil triangular durante el choque. 

d 

Desplazamiento finito lineal. 

ÓA 

Infinitesimal de área. 

ÓX 

Infinitesimal de la coordenada actual de los puntos de un sólido 

deformable o de una cuerda. 

£ 

Deformación, pequeña, lineal. 

£ 

Ratio de deformación. 

tPcorr. 

Factor de cpífécción del lí^pplazamiento aproximado (modelo 
de pequ^os desplazamientos) del punto medio de la cuerda 

Tí 

celeste delimitantel^j 1 ^^^ 

Constante 3. 141 í 9^ 

Po 

Densidad lipeal de una cuérda en su estado inicial. 

. t 

IBS Si?©!®** 

p 

Densidad lineal actual de una cuerda. 

Pmat 

Densidad del material. 

o 

Esfuerzo uniforme sobre un área. 


Esfuerzo axial o longitudinal promedio compartido por todos los 

alambres del cable en sus extremos durante el proceso de 

carga. 

a xx 

Esfuerzos en la dirección x en cualquier punto a lo largo del eje 

z con origen en el punto de impacto sobre el cable. 

°yy 

Esfuerzos en la dirección y en cualquier punto a lo largo del eje 

y con origen en el punto de impacto sobre el cable. 
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&L.MÁX 
&ÍR .MAX 
&V.M 

^R , Máx 

V 

V 


Esfuerzos en la dirección z en cualquier punto a lo largo del eje 
z con origen en el punto de impacto sobre el cable. 

Esfuerzo longitudinal máximo en el cable. 

Esfuerzo de resistencia máxima del cable. 

Esfuerzo de Von Mises equivalente en el eje z con origen en el 
punto de contacto entre el cable y el sólido que lo impacta. 

Tensión de rotura del cable, 

Módulo de corte del material de los cilindros en contacto. 
Coeficiente de Poisson. 


&L.MÁX 
&ÍR .MAX 
&V.M 

^R , Máx 

V 

V 
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1 . Introducción. 

1.1 . Realidad problemática. 

Alrededor del mundo se usan barreras de protección para disminuir las 
consecuencias fatales de los accidentes de tránsito, accidentes por caída de 
objetos en minería y de deslizamientos de suelos [1], [2], [3], [4], Las razones 
para el diseño e instalación de estas estructuras están en lo que muestran 
las estadísticas cada año como que más de 1.25 millones de personas 
mueren alrededor del mundo en colisiones de tráfico, que la minería es 
responsable del 8% de todos los accidentes laborales; o que los 
deslizamientos de tierras causaron 4626 muertes a nivel global en el bienio 
de 2015 y 2016 [5], [6], [7], [§^ppfc¡^Juestro país no es una excepción, es 

a usar barreras de protección ya que un 
50% de nuestra geografí as muy susceptible a deslizamientos de tierras y 

den ser mayores a 38.9 m/s) por 
os que se dan con frecuencia aquí 
con 59 muertes solo entre 2015 y 
2016. También hay que tener en cuente^ gran cantidad de accidentes de 
tráfico: 2253 muertes en accidentes viales en 2016 (con velocidades que 
pueden ir hasta 250 Km/hr o 69.4 m/s) y los accidentes laborales en minas 
que se dan cada año: 115 accidentes de trabajo por caídas de objetos en 
minas en 2017 (objetos con 140 m/s o que caen desde 1000 m de altura) 
[10] [11], [12], [13], [14] [15], [16], Por estas razones, aquí ya se han 
empezado a usar barreras de protección, por ejemplo, las barreras de mallas 
y cables de acero instaladas en el distrito de Lurigancho en Chosica que 
detuvieron rocas de grandes proporciones y derrumbes durante el 
Fenómeno del Niño del 2017, también están las geomallas de la Costa 
Verde en Lima que protegen de la caída de rocas a los motoristas que 
circulan en la autovía de ese malecón [17], [18], 

En las figuras 1, 2, 3 y 4 se muestran ejemplos de las barreras de 
protección en algunas de sus aplicaciones [19], [20], [21], [22], 


caídas de rocas (Con 
los terremotos y fenómenos 
como los 591 deslizamientos de 
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Figura 1: Ejemplo de una Figura 2: Ilustración del funcionamiento de una barrera 
barrera vial de cables de vial de cables CASS™ redireccionando al vehículo que 
acero [19], impacta [20], 



Figura 3: Barrera dinámica usada en la 
quebrada Carosio-Chosica conteniendo detritos 
de un deslizamiento en el verano de 2017 
detritos [21]. 



Figura 4: Ejemplo de una barrera deteniendo una roca 
que puede impactar autos si cae hacia una carretera [22]. 
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Sin embargo, estas acciones aún no son suficientes para resolver 
nuestra situación porque no solo debemos limitarnos a esperar o usar estas 
soluciones si no que se deben hacer mejores y nuevos diseños. Pero para 
diseñar y mejorar se debe estudiar el comportamiento de los elementos 
estructurales de las barreras como los cables de acero ya que estos 
elementos son parte principal de varios tipos de barreras de seguridad que 
existen, como: de las barreras viales flexibles, de las barreras para detener 
rocas y deslizamientos de suelos y de las bases artificiales que se usan en 
minas de pozo profundo. En todos estos diseños los cables son impactados 
oblicua o transversalmente por vehículos, piedras o maquinaria los que a a 
su vez son redirigidos o detenidos para evitar o aminorar daños económicos 
o mortales [1], [2], [3], [4], [23],[2$, [2^26], [27], [28],[29], [30], Debido a 



nftará en 


estudiar el impacto 
itir dar soluciones a los problemas 


esta problemática, este 
transversal en cables d 
mencionados antes. 

1.2 . Antecedentes 

El desarrollo de las ^g^otección se basa en investigaciones 

como la de C. Stolle quien experimentó y simuló con precisión un cable tipo 
3x7 sometido a impacto transversal que se usa en barreras viales flexibles 
aunque limita su modelo numérico solo a este cable y a velocidades bajas 
(6.6 m/s) [31], En el caso de barreras de protección para caída de rocas 
diseñadas con cables y malla de acero se tiene el trabajo de , Koo et al, 
quienes hicieron el estudio numérico de dos tipos de estas estructuras para 
predecir su respuesta al impacto de esferas con hasta 4000 Kg de masa y 
velocidades de hasta 36.6m/s. Validaron su modelo con datos 
experimentales aunque los resultados de su simulación difirieron de los 
medidos en los postes y cables [32], También se han hecho trabajos 
similares para las bases de protección en minas de pozo profundo como el 
de Tytko et al y P. Gospodarczyk quienes simularon un diseño de estas 
estructuras con cables de acero tipo 1x7 entrecruzados unos con otros para 
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soportar la caída de maquinaria desde hasta 800m de altura. Debido a la 
limitación del conocimiento teórico y de recursos para la experimentación en 
escala 1:1 lo que hicieron primero fue llevar acabo experimentos bajo 
impacto en el acero de la estructura y bajo impacto transversal en los cables 
en forma separada. Luego ingresaron los datos en la simulación numérica y 
así llegan a optimizar su diseño [29], [33], El trabajo mas reciente ( 2017) 
sobre cables impactados transversalmente usados en barreras de protección 
es el de Dqbek et al, en su estudio también tuvieron el objetivo de diseñar 
una base para minas de pozo profundo para detener masas con una caída 
desde hasta 1000 m, simularon una estructura con cables tipo 6x36 con 
alma de fibra probando de antemano las propiedades físicas en laboratorio. 
Al final, afirman que los cables absorben solo una pequeña porción de la 
energía y no llegan a optimizar su diseño [4f^ 

Lo que se observa en todas las investigaciones mencionadas es que 
limitan la extensión de su trab ajó a los diseños y cables simulados usando el 
programa comercial LS-DYNA. Tarrjpqpo hacen uso de modelos analíticos 
del impacto en los cable&JnvolucradoíHii hacen referencia a desarrollos 
teóricos que han venido haciéndose por muchos años como el de G. W. 
Housner (1945) quien tuvo que modelar matemáticamente el impacto de 
aviones con cables que sostenían globos usados como barrera protectora de 
una instalación nazi durante la Segunda Guerra Mundial. Su modelo predijo 
que los cables se romperían en su base sin dañar las alas de los aviones y 
eso pasó efectivamente en la práctica [34], Después, durante la primera 
mitad del siglo pasado, otros autores también hicieron formulaciones 
similares como F. O Ringleb, que modeló y probó experimentalmente los 
cables tipo 6x19 con alma de cáñamo para diseñar sistemas de frenado para 
el aterrizaje de aeronaves en porta-aviones [35], Otro investigador de esta 
época fue K. A. Rakhmatulin, que estudió teóricamente los casos de impacto 
oblicuo y con rozamiento de un cable con un punto material a velocidad 
constante y llegó a desarrollar un modelo para cuando hay una 
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desaceleración muy pequeña de la masa que impacta una cuerda con 

desplazamientos pequeños [36], Pero es N. Cristescu, quien da mayores 

alcances sobre otros estudios y de como el mismo modelo se aplica a fibras 

plásticas, además, presenta un método numérico para resolver el caso del 

impacto transversal en una cuerda con un extremo fijo asumiendo de 

antemano las condiciones de frontera asumiendo una función sinusoidal 

para el movimiento del punto material que impacta [37], Después, en la 

década de 1990, J. Lombard, se basó en el trabajo de G. W. Housner, y 

modela accidentes aéreos entre aviones militares y cables de teleféricos, la 

diferencia es que asume un perfil parabólico y establece que el cable se 

deforma solo plásticamente durante el choque, sus resultados coincidieron 

con los de G. W. Housner,_jaunque tampoco presenta algún dato 

experimental [38], [39], 

/ ^ 

plantearon una solución 
dinámicas de los cabléSTáSum 


como constante en todo instante, 



tico de todos los trabajos que 
simplificación de las ecuaciones 
idad de la partícula que impacta 
ién establecieron de antemano una 


forma (triangular o parabólica) para el q^Jale o asumieron de antemano las 
condiciones de frontera como definir la función que describe el movimiento 
de la partícula durante el choque. En la figura 5 se muestra el perfil triangular 
asumido por varios autores para el cable durante el impacto así como las 
velocidades características de sus puntos donde el punto de impacto tiene 
velocidad constante v . 


> 



Figura 5: Diagrama del cable impactado asumiendo un perfil 
triangular. v , velocidad constante del punto impactado; kv , 
velocidad de la onda transversal, c , velocidad de la onda 
longitudinal [34], 
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Y no fue hasta 2005 que X. Teng, en su modelo de la fractura de cintas 
sometidas a impacto transversal por un bloque rectangular, donde si toma en 
cuenta el cambio de la velocidad del proyectil durante el choque y no 
supone una forma para la cinta (asumida como cuerda unidimensional) pero 
también simplifica las ecuaciones de movimiento suponiendo que el 
comportamiento del material es perfectamente plástico lo que reduce las 
expresiones a solo el intercambio de cantidad de movimiento en dirección 
transversal entre la masa y la cinta que se va fracturando en forma continua 
durante su movimiento y además no hace algún contraste con datos 
experimentales [40], 


Es así que en la aGtt3^Rtíé«d^tenemos que todos los modelos 


experimentales o numéricos 
de acero, limitan su éfrnu, 
hacen uso de los estu, 
estudiar el impacto oblicuo 
constante la velocidad de impact 



arreras de protección que usan cables 
tipos de cables que estudiaron, no 
s y estos, a su vez, se limitan a 
I de cables de acero manteniendo 
asa puntual que choca en el punto 


medio del cable o con una desaceleración muy pequeña o asumiendo a 
priori la forma del cable impactado como triangular, con pequeños 
desplazamientos, además suponen una función para el movimiento de la 
partícula que impacta o establecen que se da la fractura instantánea del 
material del cable. Pero todo esto no pasa y no debe pasar en las 
aplicaciones reales de las barreras de protección mencionadas ya que ahí el 
objetivo es detener a algún automóvil que se accidenta o a rocas o 
maquinaria que caen y prevenir que alcancen a personas en las carreteras, 
pueblos donde viven o lugares de trabajo lo que implica desaceleraciones 
grandes y aún más, la estructura debe seguir manteniendo su integridad 
estructural después del choque sin llegar al punto de fractura o rotura de sus 
partes. 


Por tanto, debido a la necesidad de diseñar más y mejores barreras de 
protección y sabiendo las limitaciones que tienen las simulaciones, 
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experimentos y teorías sobre el impacto en cables; este trabajo se centra en 
modelar matemáticamente el impacto transversal entre cables de acero y 
una masa que en la realidad puede ser un automóvil, alguna maquinaria o 
una roca para determinar, principalmente, los esfuerzos que surgen teniendo 
en cuenta la desaceleración de la masa que impacta y asumiendo que el 
cable se comporta elásticamente de manera que se pueda diseñar algo que 
no llegue a romperse o fracturarse ya que la función principal de las 
estructuras de protección es justamente ese, proteger. 

Ahora, de acuerdo a todo lo dicho se plantea un problema de 
investigación. 



determinar los esfuerzos que 
pactado transversalmente por una 


1.4 . Justificaciones. 

Las justificaciones para trabajar en el problema planteado son las 


siguientes: 

1.4.1 . Justificación científica. 

Porque es necesario entender mejor y conocer más sobre el fenómeno 
del impacto transversal en cables. Esto, como tarea pendiente por investigar 
y como alternativa de análisis de los modelos numéricos actuales. 

1.4.2 . Justificación tecnológica. 

La necesidad de conocer el valor de los esfuerzos que surgen en un 
cable de acero que es impactado transversalmente por una masa que podría 
ser un vehículo, roca o maquinaria que choca con alguna estructura de 
protección vial, de protección ante caída de rocas o de detención de objetos 
que caen en minas y contar con una herramienta más para el 
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dimensionamiento y selección correcta de los cables a usar y de las demás 
partes de las estructuras mencionadas. 

1.4.3 . Justificación económico-social. 

La necesidad de dar más y mejores alternativas de solución para 
disminuir o eliminar los efectos mortales y monetarios de los deslizamientos 
de tierras, caídas de rocas, accidentes automovilísticos y caídas de objetos 
en minas. Y en el desarrollo de estas soluciones para que se pueda ahorrar 
recursos de tiempo, dinero y humanos cuando se decida hacer experimentos 
o simulaciones con prototipos de diseños que incluyan cables bajo impacto 
transversal. 


íatemático de los esfuerzos 
la zona de contacto de un 
modelo del movimiento del 

cuerpo sólido que impacta en el punto medio del cable. 

'* 101 ¡ iflSs-' 

1.5.2 . Objetivos específicos. 

Luego de plantear y solucionar el modelo se busca hacer lo siguiente: 

• Analizar la dependencia de la aceleración de la masa que impacta en 
relación a las características del cable y de la velocidad de impacto. 

• Analizar el comportamiento de los desplazamientos de los puntos del 
cable durante el impacto en relación a sus propiedades y de la masa y 
velocidad del cuerpo que impacta. 

• Analizar el valor de los esfuerzos longitudinales en el cable durante el 
proceso de carga (desaceleración) o descarga (aceleración de vuelta) 
del cable impactado. 


1.5 . Objetivos ¿ 

1.5.1 . Objetivé 

El objetivo princip 
longitudinales, otro de los es 
cable de acero impactado transa 



rnos en 
nte y un 
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• Validar el modelo matemático de los esfuerzos axiales (fuerzas) con los 
resultados experimentales del estudio que hizo C. Stolle 1 sobre el 
impacto transversal en el cable tipo 3x7 que se usa en barreras de 
protección vial y también con los resultados de la simulación numérica 
del impacto transversal en el cable tipo 1x7 hecho por Tytko et al y P. 
Gospodarczyk y que es usado en bases de protección de minas de pozo 
profundo [29], [31], [33], 


• Calcular los esfuerzos de contacto entre el cable y el cuerpo que lo 
impacta para cada experimento de C. Stolle y la simulación de Tytko et al 
y comparar los resultados con las observaciones cualitativas de la zona 
de impacto en el cable trabajos mencionados [31], [29], 

1.6 . Alcance de 

En esta tesis se 
mecánicos que se dan en el i 



matemáticos de los fenómenos 
svérsal de cables de acero por una 
masa puntual. Este estudio se limita al impacto transversal en el punto medio 
del cable, con un comportamiento elástico , con su geometría como una 
cuerda unidimensional deformá'ble y sin resistencia a la flexión pero si 
incluye la deformación permanente por la histéresis del cable. El sólido que 
impacta y se idealiza: lo que puede ser un auto, piedra o maquinaria que se 
accidenta o cae se representa como un punto material de manera que no se 
tiene en cuenta otro efecto más que su desaceleración en una trayectoria 
recta despreciando el rozamiento con la superficie del cable. El 
planteamiento analítico comprende la formulación de modelos específicos 
según las medidas tomadas en los experimentos y en la simulaciones, es 
decir, un modelo para la fuerza en los extremos de la cuerda, un modelo del 
movimiento del objeto que impacta y también se incluye un modelo de los 
esfuerzos internos que surgen dentro del cable en la zona de impacto por el 
contacto con el cuerpo que lo impacta. Después se validan los modelos 


1 C. Stolle lisa la nomenclatura DBC-1,2,3... para denominar a sus experimentos con cables 
impactados transversalmente y DTC-1,2,3.. para los experimentos de impacto longitudinal [31], 
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aplicando las ecuaciones desarrolladas e idealizando cada uno de los 
experimentos y simulaciones hechas por 3 autores cuyas investigaciones 
fueron sobre barreras de protección usando cables 1x7 y 3x7. Después 
vienen las conclusiones, recomendaciones y al final dos anexos. 


1.7 . Organización de la tesis. 

Se tienen 8 capítulos estructurados de la siguiente manera: 


El capítulo I, comprende la realidad problemática, antecedentes, 
planteamiento del problema, justificaciones, objetivos, el alcance y la 
organización de la tesis. 


En el capítulo II, se d t a el íña/^teórico que abarca conceptos básicos 


sobre el impacto entre 

/<§3? r 

cuerdas, las ondas de 
de una cuerda elástic 
esfuerzos longitudinales y 



esplazamientos y deformaciones en 

- m 

al, las ecuaciones de movimiento 
de un cable de acero sometido a 
ripción de la física del impacto 


transversal en cuerdas elásticas. 

\o (jrtx 

Luego en el capítulo III, se.'éSíSSlece los materiales, métodos a seguir, 
se señala las variables de los modelos matemáticos desarrollados y se dan 
las restricciones, simplificaciones geométricas y físicas del fenómeno 
estudiado así como los procedimientos para la formulación y la validación de 
los modelos matemáticos planteados. El capítulo IV, comprende el 
planteamiento de los modelos matemáticos para la tensión en los extremos 
del cable, para el movimiento del objeto que impacta y para los esfuerzos 
internos provocados por la fuerza de contacto durante el choque. 


Después en el capítulo V, se validan los modelos usando datos 
experimentales y simulados de tres trabajos anteriores hechos sobre 
barreras de protección que usan cables 3x7 y 1x7, luego en VI se dan las 
conclusiones y recomendaciones. Y para finalizar en el capítulo Vil, se da la 
bibliografía seguida del capítulo VIII, donde hay dos anexos. 
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2 . Marco teórico. 

Los temas relevantes en este trabajo son el impacto entre sólidos 
deformables, la deformación, desplazamientos, ondas de esfuerzo y 
dinámica de barras y cuerdas sometidas a esfuerzos externos y el 
comportamiento y propiedades de los cables de acero. 


2.1 . ¿Cuándo es relevante el impacto ? 

Antes de pasar al estudio de cables bajo impacto, propiamente dicho, 
es necesario dar algunos alcances sobre el fenómeno del impacto en 
general ya que su estudio no es muy común en nuestro medio. 


El impacto entre sóli$áé v Heforf 

♦¿r A k 

discontinuidades en la¿ prc 

/i§5r 

otras cantidades termodir 
por algún intercambie 
estas discontinuidades se 
discontinuos se pueden dar por 



tes se observa cuando se presentan 

% 

de los materiales como la densidad u 
los grandes cambios en la presión 
ntidad de movimiento, además, 
en los cuerpos. Estos cambios 


is percusiones y los percibimos con 
nuestros sentidos como eventos ,^^fiñíbs o repentinos, sin embargo, no 
existe una definición absoluta de lo que “discontinuidad” o “evento violento” 
significa en términos de duración de tiempo por lo que cualquier concepto 
sería algo solamente subjetivo. Por ejemplo, en el mundo intelectual de la 
física de choques una discontinuidad típica sería el frente de onda que viaja 
desde la zona de contacto hacia todos los puntos de un sólido como cuando 
dos barras de metal colisionan entre sí, aquí la presión de choque en el área 
de contacto puede aumentar en décimas de nanosegundo. Otro ejemplo se 
da en el ámbito de la geología, donde los terremotos pueden provocar 
cambios drásticos e irreversibles en varios segundos. Curiosamente estos 
casos, aunque se dan en escalas de tiempo y dimensiones muy diferentes, 
lidian con los mismos órdenes de magnitud cuando se los relaciona con sus 
intervalos de vida tomados como sus cantidades características de la 
siguiente manera: Supongamos que un medio sólido tenga una escala 
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característica de l (medida espacial, longitud de una barra metálica por 
ejemplo) y la velocidad de propagación de las onda elásticas dentro del 
material, o velocidad del sonido, c . Luego se somete al cuerpo a una carga 
dinámica externa que tiene un tiempo característico t c que es, por ejemplo, 
el periodo de tiempo de la carga para alcanzar su valor máximo o el tiempo 
de duración del impulso; si t c cL/c , entonces la distribución del esfuerzo y 
la deformación en este sólido no son uniformes, por tanto, el efecto de la 
propagación de la onda de esfuerzos debe de considerarse. Intervalos de 
tiempo de estas características son, por ejemplo, 10 nanosegundos. que 
corresponde a un 4.2 x 10' 18 del intervalo de vida de un humano de 75 años 
(o una generación o la estructura) mientras que un periodo 

de 1s corresponde a un 6.9 x 10^|de la edad de la tierra (unos 4.6 billones 
de años o puede serjína era geológica. coit¡,o el Cretáceo o el Mesozoico). 
Son estos tiempos tan pequeños la razón por la que los modelos de impacto 
en sólidos como las cuerdas 'asume 
entre un cuerpo con cierta velocidad 


»ro|üe luego del instante de contacto 
I y algún punto en la superficie de un 

cable de acero adquiérala misma vejQfeidad, o lo equivalente, se asume 

__ _ 


que el esfuerzo en la cuerda adéjliífere un valor no nulo en forma instantánea 
modelado como una función de Heaviside [35], [40], [41]. Sin embargo, lo 
anterior solo nos dice que al menos todo inicio de contacto entre sólidos 
genera una onda de esfuerzos que se propaga y según su tiempo de 
duración se puede decir que es un impacto. Pero no nos dice si los 
esfuerzos son lo suficientemente relevantes como para considerarlos 
determinantes en la integridad estructural del cuerpos deformables que 
impactan. Para esto consideramos otro criterio, la llamada “Caracterización 
de Weirauch”, cuyas dos primeras escalas indican que para ratios de 
deformación (velocidad de deformación) del orden de 10° y velocidades de 
choque menores a 50 m/s se tiene un efecto plástico localizado y a partir de 
50 m/s hasta 500 m/s con ratios de deformación del orden de 10 2 se tienen 
efectos puramente plásticos. Claro que estos parámetros solo son puntos de 
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referencia flexibles ya que no se consideran muchos otros factores como la 
geometría, el calor generado, rozamiento o la similitud de resistencias entre 
los materiales del proyectil y su objetivo, la respuesta de los materiales a 
diferentes ratios de deformación etc. un ejemplo de esto es el choque de dos 
cuerpos dentro de la primera escala de Weirauch, si uno es de acero de alta 
resistencia y el otro de aluminio 1100 este último desarrollará un flujo 
plástico típico de los impactos de hiper-velocidad mientras que el acero solo 
sufrirá una perturbación elástica [42], [43]. Por tanto, es imperante analizar 
cada caso y establecer cuales son los parámetros importantes. Un ejemplo 
es el estudio de de fibras de nylon hecha por N. Cristescu, donde llega a 
establecer que cuando la velocidad de impacto excede 1/20 de la velocidad 
de propagación de las o^jiás transversal^ o ratios de deformación de 33.33 

ción de estas ondas o mejor dicho el 


m/m/seg, entonces el^íec 

/ <5 ^r 

efecto local de las fuer; 
estos casos no se p 
cuerdas deformables 



tomarse en cuenta, es decir, en 
uaciones de movimiento de las 
priori como la forma triangular del 


cable asumida por otros^utores [37], [34J r [38], En el caso de los cables de 

’íijj» _ 

acero no hay un criterio para*íésfóbfécer cuando tomar en cuenta la inercia 
local según la velocidad de impacto o el ratio de deformación. Esto 
representa un campo poco explorado. 


2.2 . Desplazamientos y deformaciones en cuerdas. 

Para estudiar el impacto en cuerdas es necesario también aclarar los 
conceptos sobre los desplazamientos y deformaciones en este tipo de 
elementos estructurales ya que un cable puede tener grandes 
desplazamientos junto a pequeñas deformaciones. Sobre esto A. E. H. Love 
indica que para barras, y por extensión a las cuerdas, es posible tener 
desplazamientos relativos de las partes de barras o cuerdas largas no 
pequeños y aún así las deformaciones ser lo suficientemente pequeñas 
como para satisfacer los requerimientos de la teoría matemática tanto en 
Resistencia de Materiales como en la Teoría de la Elasticidad [44], Entonces, 
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por una parte tenemos la aproximación de “deformaciones pequeñas” que 
nos permite considerar la respuesta elástica de un material como lineal y por 
otro lado está la aproximación de los “desplazamientos pequeños” que tiene 
que ver con las rotaciones del material: los elementos de volumen pueden 
tener rotaciones de cuerpo rígido y ser irrelevantes para calcular las 
deformaciones elásticas, pero, si se quiere considerar las rotaciones se tiene 
que introducir la no-linealidad en nuestras ecuaciones [45], La condición de 
“desplazamientos pequeños” es usual encontrarla en las aplicaciones 
convencionales de ingeniería, el ejemplo clásico es la linealización de la 
relación del momento flector y la curvatura en una viga o barra despreciando 


el cuadrado de la inclinación ( dy/dx ) 2 al compararla con la unidad que se 

. 

puede hacer siempre y.jpíando^g^considere deflexiones muy pequeñas con 
una curva final de lajO'igq 


demás casos de barr; 
comparables con su propi 
elemental. Por tanto, para mode 



lai*J£ [46], Por lo que en todos los 
bres donde las deflexiones son 
se puede aplicar esta teoría 
rras, alambres o cuerdas se debe 


saber de antemano que va a haber grabes desplazamientos lo que lleva a 






la no-linealidad de las ecuaciones diéDido a la geometría, mientras que la ley 
constitutiva del material nos lleva a poder restringir el estudio solo a 
deformaciones pequeñas. 


Lo dicho se debe complementar con la definición, en términos 
matemáticos, de la deformación en un cuerpo unidimensional como la 
cuerda, que es el sólido elástico objeto de este estudio. Se puede elegir 
entre varias definiciones de deformación, pero para un régimen de 
deformaciones pequeñas (en el caso elástico de los metales) se puede 
comprobar que todas las demás definiciones resultan ser iguales a la de 
Cauchy [47], También se debe aclarar que aunque el impacto involucra 
esfuerzos dinámicos esto no necesariamente implica que la respuesta del 
material difiera de cuando se considera una aplicación estática. Von Karman, 
notó esto en su estudio del impacto longitudinal con esfuerzo plásticos en 
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una barra señalando que la dependencia de los esfuerzos del ratio de la 
deformación se da en velocidades muy altas pero como no especifica 
exactamente cuan altas uno puede tomar la referencia de Weirauch, para 
establecer un criterio de referencia para la dependencia de la velocidad de 
deformación [43], [48], [49] . 


2.3 . Ondas de esfuerzo unidimensional. 

Las ondas de esfuerzo unidimensional aparecen sin importar si se 
considera impacto o no a la interacción con otro sólido, solo es cuestión que 
haya un cambio de posición relativa entre los puntos del sólido. Para hallar 
una expresión 


matemática sa 

representa el modelo de 

cí 

las barras prismáticas 
cuerdas rectas en figt 
6 donde se muestra una 


ByU 

u + ^=ÓX 



barra 


de 


a rea 

transversal variable con 
su eje axial coincidente 
con el eje coordenado 
de referencia x, con 


A a 



{A+SA){a+p^dX) 


8yX 


m 


n 


Figura 6: Equilibrio dinámico de una barra deformable 
con sección transversal variable, x , es la coordenada de 
origen en O [50], [51]. referencia, X , es la coordenada actual; un elemento 

Luego en el equilibrio inft ™ tesimo1 de lon S' wd sx ■ acciones mm' y nn', 

desplazamiento u , cambio infinitesimal de su área ó A 

dinámico y con un ) con un esfuerzo uniformemente distribuido a y la 

proceso isocórico se velocidad del frente de la onda de esfuerzos o del sonido 

^ ’ c [ 50 ], [ 51 ]. 

tiene: 


A(l+e) = A 0 ; A 0 p„ 


df u _ d x ( rA d x u \_ / A 0 d x u v 
d t t 2 d x x d x x d x x l+sd x x 


( 1 ) 


En la ecuación (1) se tiene en cuenta el cambio del área de la sección 
transversal debido al estiramiento en la configuración actual pero como lo 
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indica Timoshenko, esta no toma en cuenta la inercia de las partículas que 
se mueven en dirección lateral debido a la contracción por lo que la ecuación 
es válida solamente cuando la amplitud de la onda de esfuerzo no sea tan 
pequeña como para compararse con las dimensiones transversales de la 
barra [52], Aún así, la ecuación se hace mas simple si es que se considera 
una barra sólida de sección transversal constante, resultando en: 


afu 2 

— C 


dtu 


•-'¡up. 


d t t d x x 
Siendo su solución general la de D’Almabert: 


( 2 ) 


u(x,t) = f 1 (x—ct)+f 2 {x+ct) (3) 

2.4 . Ecuación generaMdel movimiento de una cuerda 
elástica,^ 

La teoría del mock^ a desarrollar en este trabajo tiene su base en lo 
hecho por N. Cñstescu, q|uien derivo las ecuaciones de movimiento de una 
cuerda extensible que se mueve en un espacio tridimensional [37], [53], para 
esto usó una coordena(%,curvilínea S ^Ho largo de la cuerda empezando 
desde un origen arbitrario, cómtePefíW figura 7 , a esta coordenada Euleriana 


T(s 2 ) 


la llamaremos “actual” y corresponderá a la cuerda deformada en el 
momento actual t. La fuerza total en la cuerda es denotada por T , (se 
toma la fuerza ya que la definición 
de cuerda unidimensional implica 
que no tiene área), varía a lo largo 
de la cuerda y en el tiempo siendo 
tangente en todos sus puntos al ser 
perfectamente flexible. En la misma 
figura se considera una porción de 
la cuerda entre las coordenadas r J, S| ) 



A 


curvilíneas y s 2 , sometida a 


Figura 7: Segmento infinitesimal de una 
cuerda en el espacio sometida a fuerzas 


una tensión en cada extremo y internas y externas [37]. 
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también a las fuerzas externas (x, Y, z ), las cuales son proporcionales a la 
longitud de la porción considerada. Entonces el equilibrio dinámico se 
expresa como sigue: 


z 

/ 


d,s 


d c x\ d 2 x 

+ X-p 1 


d'S 


dX 


ds = 0 


o s 


1 .. ..3 s x\ d]x 

\\T\\f^\+X=p- t 


dX 


(4) 


Junto a dos ecuaciones similares para y y z . Aquí, p es la densidad 
lineal actual de la cuerda y es función de s y t. Las proyecciones de T en 
los tres ejes considerados son ||T||(d s x/d s s), ||T||(d Si y/0 s s), 

||r||(S s z/d s s), siendo ||T|| el módulo de T teniendo en cuenta que y 
s 2 son arbitrarios. Tambió^ri'é^ijec/^eemplazar la coordenada curvilínea s 


por la coordenada mati 


coordenadas x , y , z 
tiene esta relación: x 



ngiana s 0 y así expresar las 
última dado que para cada una se 
manera: 


d s x 5 s x d s s 

-— , ídem para y, z. 
d s„ . 


d s*0 


s n 0 

i® 


(5) 


Igual aquí se puede escribir ecuaciones similares para y y z . Luego 
se define la deformación en la cuerda de la misma manera que se ha 
definido anteriormente, es decir, según Cauchy: 


£ — 


d s-d.s 0 


d s „ s « 


J(M 

2 

+ 


2 

+ 

o 

Z 



S s S 0 

• 7 






- 1 


( 6 ) 


Idem para y, z 

Y por la ley de conservación de la masa: pd s s=p 0 d s s 0 => p 0 =p(l+e) 
aquí, p 0 , es la densidad de la cuerda en el estado inicial no deformado. Con 
esto la ecuación del movimiento (4) queda así para x (similar para y y z ): 


¿Ls 


llrll d *x 

l+ed.soj 


-Po 


d 2 t x 

dX 


+(l + £)x = 0 


(7) 
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Ahora se puede agregar la ley constitutiva del material, que usualmente 
es una función monótona creciente, se define como ||r||=||r||(e) . Entonces 
con la relación (6) y las ecuaciones para x, y y z de (7) se tiene el 
siguiente sistema: 


d So x d]x _^ d s j d\y d s z d)z 


a s» s o d s / 0 d s s 0 e s / 0 d So s 0 d s / 0 

d 2 t x d s x (l+ e )d s \\T\\ld s e- 


/ N d s e 

(l-e)^- = 0 


d] x 


(l+e)0 So s o Po dt 2 'd So s 0 (l+e) 2 

Idem para las direcciones y y z. Este sistema de ecuaciones 
diferenciales parciales es hiperbólico de orden 2 y tiene cuatro familias 


5 So s 0 

d.s . , 

—+ (l + f-)X = 0 


( 8 ) 


8 s 0 S 0 


distintas de curvas caracterísíiGaSPl!/gs 4 relaciones diferenciales son: 


d,s 0 _ 

d t t~ 

Por esto es qu 
cuya inclinación depende 



i MUI 

Po d e £ 


1/2 


(9) 


ticas son generalmente curvas 
Desde el punto de vista mecánico 


estas líneas representapJos frentes de or^la por lo que se dice que hay dos 
que se propagan con las sigútente^fzéf^cidades: 


•i(c) = 


Po( 1+£ f 


1/2 


: iAe) = 


1 d, 


P o d c s 


1/2 


( 10 ) 


Siendo 


■ i ’ 


la velocidad de las ondas que se propagan 


transversalmente y afectan solamente a la forma de la cuerda mas no a la 
deformación longitudinal, es decir, la velocidad del punto material alcanzado 
por la onda transversal será perpendicular a su tangente. Mientras que c n 
es la velocidad de las ondas que se propagan longitudinalmente y solo 
producen elongaciones que se propagan a lo largo de la cuerda sin modificar 
la pendiente de la tangente de sus puntos. 
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2.5 . Módulo de elasticidad de los cables trenzados 
de acero 

Como la teoría de la sección anterior se va a aplicar a cables de acero, 
es necesario saber sobre las propiedades de estos 
elementos mecánicos. La elongación de un cable Alma 
depende, por supuesto, del arreglo y de las 
propiedades de los materiales de sus partes (Figura 
8) especialmente de sus alambres (si el alma es 
sintética no se toma en cuenta). En general el 
módulo de elasticidad del cable es diferente al del 
material de sus partes. La curva esfuerzo- 
deformación de todos lo^Bables es no- lineal, por cable I 

tanto, dado un cierto ti 

/ o 


elástico no será consta 
tensión. Por esto en 




Cordón 


su ^módulo 

i del nivel de Figura 8: Elementos de un 
e habla cable de acero [54], 


de un módulo elástico medio. Sin embargo, si se trata de cuerdas simples 
formadas solo por alambres trenzados éfi espiral se da una no-linealidad 
mínima y puede ser desprecrátíSE ©mía mayoría de casos [55] [56], Para el 
caso de los cordones se puede calcular aproximadamente el módulo de 
elasticidad usando solamente métodos analíticos, pero esto no es tan exacto 
para los cables compuestos o de varios cordones para los que se necesita 
datos experimentales ya que la curva no-lineal de esfuerzo-deformación de 
los mismos hace que el módulo resultante pueda darse solo especificando 
su carga actual. También el límite elástico se puede definir como aquella 
tensión hasta la cual un cable no recupera su forma inicial luego de 
destensarlo. Las definiciones del módulo de elasticidad que son de 
importancia práctica son [55][56]: 


• E s (a ¡ower , a upper ) , es la pendiente de la secante entre los esfuerzos, 
final e inicial. Se aplica en el caso de esfuerzos fluctuantes, un caso 
especial es cuando el esfuerzo inicial es cero E s (o,a upper ) . 
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• E tup (o),E t down (a), serían, respectivamente, la pendiente de la 
tangente cuando aumenta la deformación y cuando baja la tensión en 
algún punto. 



Los esfuerzos en un cable de acero se definen como, a=TlA, aquí 
A es la sección transversal metálica de todos los alambres. Es importante 
recalcar que el módulo de elasticidad calculado con las mediciones no es de 
importancia práctica si es que para su cálculo no se destensa el cable hasta 
el esfuerzo inicial a partir del cual se quiere hacer el cálculo. Esto se debe a 
los fenómenos internos que suceden en el cable como el acomodo de los 
alambres entre sí o el rozamiento y no necesariamente a la deformación 
elástica de los mismos, pp|&Í^^ÍÉ^s se ha establecido que para medir 
experimentalmente el | 
de un cable este debe 
10 ciclos para estabili 
lado también se puede usar 
la curva para hallar el módulo 
práctica se usa solo como parám^tR,^ de 
evaluación y comparación. En la figura 9 
se puede apreciar la interpretación 
geométrica de las definiciones de los 
módulos de elasticidad para los cables, la 
pendiente de la línea OB sería el módulo de 
elasticidad E s (o,a upper ) entre un esfuerzo 0 

Figura 9: Definiciones del módulo 
y el esfuerzo al que trabaja el cable. Si es de elastícidad como pendientes de 

que hay un esfuerzo inicial (punto E de la las secantes OB, CB, y EF o de las 

, i , . , i , tangentes en O o A. o , es el 

cuerva) y uno final (punto F de la curva), ~ , , ’ 

esfuerzo, e t , es la deformación 

entonces el módulo elástico sería igual a la fjna/> £c > es ¡ a deformación 

pendiente de la línea secante E s (a E ,o F ) . durante la carga y e s , es la 

deformación de descarga; cp , son 
En este caso el módulo correcto o práctico /os dngulos de /flS secantes Q 

sería aquel medido como el de la pendientes. [56]. 
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pendiente de la línea CB donde se ha descargado el cable y se mide el 
módulo desde la deformación permanente que tiene este. Luego están las 
líneas E t 0 (a),E t A (a) en el punto Ot y A que son tangentes y solo sirven 
para describir la curva o hacer comparaciones [55], [56], En suma, la 
respuesta de un cable al esfuerzo axial cambia de acuerdo al tipo o arreglo 
de sus alambres y a sus demás componentes como el tipo de alma: de fibra 
u otro cable, también influye si tiene plástico entre los alambres o la forma de 
estos. 


En la figura 10 se ¡lustra, en modo general, el comportamiento no-lineal 


descarga de una fuerza 
son los que tienen un comp 
los cables de una solé c 
en la mayor parte d 
cuenta que luego de la pr 



de varios tipos de cables nuevos cuando son sometidos a la carga y 

r ‘' s plastificados o con alma de fibras 
completamente no-lineal mientras que 
e tienen un comportamiento lineal 
simismo, se tiene que tener en 
a todos los cables presentan un 
comportamiento diferente con iM^ftíddulo elástico mayor, por esto se 
recomienda someter al 10 ciclos hasta que se tenga una respuesta 

estable. En la misma figura se presenta la idealización lineal de un ciclo de 
carga y descarga de un cable nuevo con las líneas rectas [31], [55], [57], 
[58], [59], [60], [61], 

O t 

Barra de acero sin 

trefilar. .......... 

Idealización lineal de 
una cable de acero. “ ” 

Cables de la series 1x, 

3x. (de una sola trenza _____ 
y de 3). 

Cables de mas de 3 
trenzas: 4x, 6x, 

8x,..18x,.., plastificados. 

y alma de fibras. 

►S 

Figura 10: Curvas esfuerzo o, deformación e, de varios tipos de cables de acero, 
nuevos, ante una primera carga y descarga. [31], [55], [57], [58], [59], [60], [61]. 
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2.6. Física del impacto transversal en cuerdas 
elásticas. 

En esta parte se usa la geometría y el castellano, evitando a propósito 
la matemática para poder comunicar la física involucrada que se dá cuando 
un cuerpo cualquiera impacta transversalmente en un cable en un cable de 
acero y también de lo que pasa en la idealización de este fenómeno, o su 
analogía cuasi-isomorfa, que viene a ser una cuerda unidimensional sin 
resistencia a la flexión e impactada transversalmente por un punto material. 


Primero, la física del impacto transversal en un cable real se describe 
como sigue: 


1 


. Cuando un cuerpo chocEL transversal mente 


diámetro es varias 


entonces el cont 


varios de los al 
fuerzas al resto de la é 

velocidad relativa.entre los 

. 

t 



con un cable (cuyo 


lor a cualquier dimensión del cuerpo) 
sobre la superficie de alguno o 
se deformarán y transmitirán las 
a. Desde el primer instante habrá una 
tos en contacto del cuerpo que impacta 


y la superficie del%^@ £ ^s^^contacto se transmite a todos los 
alambres en cuyas secciones transversales surgen esfuerzos 
adicionales a los esfuerzos longitudinales. 


2. La deformación longitudinal que empieza a darse en los alambres 
corresponde a la de un cuerpo con hélice doble o simple (como un 
resorte común), dependiendo del tipo y ubicación en el cable. Es 
decir, no solo hay esfuerzos de tensión, sino también de torsión, 
flexión y cortantes en cada alambre [55], 


3. Por lo anterior, las ondas de esfuerzos que surgen, son también de 
torsión, flexión y cortantes Todas estas perturbaciones viajan a lo 
largo del eje curvo de los alambres y se suman entre si para dar el 
desplazamiento en conjunto de todo el cable. Además, en la superficie 
de cada alambre aparecen ondas superficiales provocadas por el 
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rozamiento. En conjunto todas las ondas llegan los extremos del cable 
donde se refractan y reflejan. Las ondas reflejadas vuelven y se 
suman a los esfuerzos actuales de cada alambre que a su vez van 
cambiando mientras cambia la fuerza de reacción que se da entre el 
objeto que choca y el cable. 

4 . Las fuerzas rozamiento internas entre alambres no dependen de un 
coeficiente constante; son más bien son funciones de estado de los 
materiales en contacto [62], 


5. Además se tiene que tener en cuenta el rozamiento del medio y la 
gravedad que también cambia la geometría del cable. 


ar 

Por tanto, si en una inves 
de un cable real, el ffíod 
externas e internas me 
tener presente la interacción 



ión no se quiere simplificar la geometría 
indlpir: las deformaciones, fuerzas 
ión y refracción en los extremos y 
que choca por que su superficie 
del cable. Pero, este planteamiento 


puede rozar o abarcar una cierta I 
aún está fuera del alcance del autor de este trabajo por lo que aquí se 
modela al cable como una cuerda unidimensional y al cuerpo que impacta 
como un punto material. 


Entonces para idealizar el movimiento de un cable se hacen las 
siguientes simplificaciones generales: 

1. El cable se reduce a su eje geométrico cuyo estado inicial es una 
curva recta unidimensional con densidad uniforme en un plano y en 
reposo. 

2. Solo se toma en cuenta la deformación axial y se desprecia los demás 
esfuerzos internos en los alambres del cable, también se establece 
que está relacionada con los esfuerzos axiales por la ley de Hooke 
con un módulo de elasticidad constante del cable . 
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3. Se reduce al cuerpo que ¡mpacta a un punto material con una 
velocidad dada que empieza a desacelerar cuando toca el cable. 

4 . Se considera que el punto de la cuerda, donde impacta el punto 
material, adquiere una velocidad igual a la de este en forma 
instantánea. 


La última simplificación implica es una consecuencia de que la cuerda 
es deformable y no tiene resistencia a la flexión. Este cambio discontinuo de 
la velocidad también significa que la geometría de la cuerda presente un 
punto singular, acodado o discontinuidad de su tangente en el punto de 
impacto. Esta singularidad es el frente de una onda transversal y su 
velocidad se modela con el término ojéele la ecuación (10) y provoca que 


los puntos contiguos al del i 
transversal lo que signifi 

Por otro lado la propa 
longitudinal y la velocidad con 



paulatinamente un movimiento 
inuo de la masa en movimiento. 

estiramiento se da en forma 
e se propaga depende de cuán 


deformable sea la cuerda, en este caspf'de cuán elástica se modela con 
C n de la ecuación (10) y provoca que los puntos se separen entre sí pero 


no alteran la dirección de su tangente como la onda transversal. 


En la figura 11 se muestra el cambio de forma de una cable idealizado 
(antes que se de la reflexión de cualquier onda) cuando es impactado 
transversalmente por una masa puntual. El frente de onda transversal tiene 
una velocidad C¡ que varía con la tensión, la velocidad del frente de onda 
longitudinal depende del cambio de la tensión respecto de la deformación, 
en el caso elástico, es constante C u . La velocidad inicial, o en el instante 
del impacto, del punto con masa M es V 00 , V es la velocidad actual, S 0 , 
es la coordenada inicial (Lagrangiana) de los puntos y S es la actual. El 
espacio donde se da el movimiento no ofrece ninguna resistencia ni hay 
fricción. 
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Frente de onda longitudinal 


M? V oo Frente de onda transversal 
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X 

Idealización del cable 
V Como cuerda unidimensional 


Figura 11: Movimiento de una cuerda impactada transversalmente por una masa 
puntual. M, es la masa del cuerpo, V 00 , es su velocidad inicial, V, es su velocidad 
actual, O, es el origen de coordenadas de ejes x, y, S, es la coordenada curvilínea 
actual de la cuerda y S 0 , es la coordenada de referencia. 

Entonces, idealmente, la física del impacto transversal de un punto 


material con una cuerda^ 
etapas: 


W 


1. Etapa 1: En el 

punto material y ün pu 
la misma velocidad. Así er 
tensión. 



Unidimensional tiene las siguientes 


duce el contacto inicial entre el 
a e instantáneamente adquieren 
la etapa de carga o aumento de la 


2 . Etapa 2: en la segunda etapa el cuerpo empieza a desacelerarse y 
desde el punto de contacto una onda longitudinal viaja a la velocidad 
del sonido (según el material) en la cuerda, a la vez una onda 
transversal, que solo cambia de dirección a la tangente de los puntos 
por los que pasa, se propaga hacia los extremos. Esta última onda es 
la manifestación de la inercia de la cuerda que es acelerada aunque 
la onda longitudinal también contribuye, en menor medida, a esto 
debido a que estira el cable y cambia la densidad en cada punto. 


3 . Etapa 3: Si el cuerpo sigue desacelerándose y el frente de la onda 
longitudinal iniciado en el punto de contacto llega a un extremo se 
refleja o refracta de acuerdo con lo que tope, luego la onda reflejada 
se añade al estado de esfuerzos que también va cambiando según la 
fuerza ejercida por el objeto que impacta, es decir, la onda de 
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esfuerzos producida inicialmente en el contacto vuelve e interfiere en 
el movimiento. La onda transversal también se refleja o refracta y 
empieza a influir en el movimiento de los puntos de la cuerda y del 
cuerpo debido a que esta onda interactúa con la inercia del cable por 
donde pase. 


4 . Etapa 4: El cuerpo adquiere una velocidad cero y empieza a cambiar 
el sentido de su movimiento. Así finaliza la etapa de carga y empieza 
la descarga. Entonces se tiene que tener en cuenta si es que el punto 
material se ha pegado o no a la cuerda, por que si esta ha tenido 
alguna deformación plástica, entonces en algún momento, durante el 
retorno, el punto mat^ri^féíAQlfg, una velocidad mayor a la velocidad 
del punto de contacto 


se disipe. Acá¿?se 
cable vuelva a s 



vimiento seguirá hasta que su energía 
viiíriento hasta que la tensión del 
e finaliza la etapa de descarga. 


. Materiales y métodos 
3.1 . Recolección de infom>§€ión 


La información necesanéf^ffe^'plantear el problema y solucionar las 
ecuaciones se encuentra, principalmente, en los trabajos hechos por F. O. 
Ringleb, N. Cristescu y K. Feyrer [35], [37], [53], [55], La validación con datos 
experimentales y simulados se hace con los trabajos de C. Stolle, Tytko et al 
y P. Gospodarczyk [31], [29], [33], 


3.2 . Material de estudio. 

El objeto de estudio de este trabajo es un cable anclado en ambos 
extremos, hecho completamente de acero con torcido helicoidal y alambres 
de perfil circular. Adquiere movimiento debido a que es impactado 
transversalmente por una masa puntual. Para su análisis matemático se 
idealiza su geometría y movimiento aplicando varias restricciones y 
simplificaciones. Y para validar los modelos matemáticos planteados se usan 
los datos del comportamiento bajo impacto transversal de los cables tipo 3x7 
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(especificado por ASTM A741 -11 (2016)) y 1x7 (Especificado en BS 5896- 
2012) usados en las investigaciones de C. Stolle, Tytko et al y P. 
Gospodarczyk respectivamente [31], [29], [33] . 


3.3 . Metodología a emplear para el planteamiento de 
los modelos matemáticos. 


El método a seguir en esta investigación forma parte de las 
metodologías del modelamiento matemático y estas a su vez del 
modelamiento en general el cual consiste en relacionar el mundo real con 
alguna representación simplificada del mismo tanto física como conceptual 
yendo y viniendo entre lo real e ideal en forma cíclica. Específicamente está 
enmarcado en la filosofí§$ét cambio de Newton, y en lo desarrollado por J. 
Olivencia 2 . 


Entonces, el prirn 
establecer restricciones 



ntear un modelo matemático es 
del fenómeno real y luego se 
¡enta o lenguaje para expresar lo que 


puede usar la matemática como he 

hay en el mundo real id%lizado [63], En-Jfc primera parte, que consiste en la 

- 

o indicado por P. Teodorescu, aquí se 


idealización de la realidad se 


traslada lo observado en el impacto transversal en cables de acero en los 
varios tipos de barreras de protección a una realidad ideal dentro de la 
mecánica de sólidos deformables modelando los cuerpos físicos en su forma 
más simplificada: las barras son reducidas a su ejes o las placas a sus 
superficies medias, o un cable de acero reducido a su eje longitudinal. 
Correspondiéndose así con lo enmarcado en la Resistencia de Materiales 
[64], Luego en el mundo ideal y con la necesidad de establecer una teoría se 
puede seguir el paradigma de la Mecánica Racional y así plantear el método 
para el modelo matemático de un cable sometido a impacto transversal por 
una masa puntual como lo indica A. Klarbring de manera que el método a 
seguir tiene los siguientes pasos [63]: 


2 Profesor de matemática de la Universidad Nacional de Trujillo. 
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1. Primero se idealiza el fenómeno real al definir las restricciones, 
simplificaciones geométricas y físicas. Para esto se describe con 
lenguaje geométrico y matemático el cuerpo a estudiar y sus 
elementos. Los puntos materiales: Se pueden tener 4 modelos de 
cuerpos materiales: discretos, unidimensionales, bi-dimensionales 
y tridimensionales. 


2. Segundo, incluimos las leyes universales que son: conservación 
de la masa, energía, momentum lineal y angular; con esto se 
pueden dar expresiones en términos de ecuaciones diferenciales 
de las leyes universales según el cuerpo material definido. 


ti 

3. Y el tercer paso para obtener una expresión matemática 
solucionable es la intfóducción de las leyes particulares: 
restricciones leyes de fuerzas y las leyes 

constitutivas de los materiales (estas pueden ser experimentales 

y/ 

o idealizadas). 

Por otro lado, la metodología para validar el modelo consiste en 
comparar cuantitativa o cualitativamente los resultados experimentales o 
simulados con los teóricos y luego analizar las razones de sus diferencias y 
coincidencias de acuerdo a la influencia de los parámetros escogidos para el 
modelo, características del cable y objeto que choca. 


3.4 . Variables de investigación. 

Si se clasifica los modelos matemáticos de este proyecto de acuerdo a 
K. Velten y a F. E. Cellier, se tiene que de acuerdo al sistema en estudio 
cada modelo es: técnico, determinista, continuo, bi-dimensional y físico. 
Tomando en cuenta lo que se quiere responder se tiene que cada modelo es 
mecanista, no-estacionario (en el tiempo), distribuido (en el espacio), 
inverso, de investigación, de análisis y en macro-mecánica. Por último, la 
matemática es: no-lineal, numérica, continua y con ecuaciones diferenciales 
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[65], [66], En suma, los modelos son tipo caja gris y se representan en la 
figura 12. Las variables de entrada son las siguientes: 

1. Masa del punto material que impacta, M c . 

2. Velocidad del punto material al momento del impacto, V 00 . 

3. Ley constitutiva del cable a modelar: ley de Hooke para la carga y 
descarga, o-Ee. 

4. Módulo de Poisson del cable, v. 

5 . Área transversal o metálica del cable a modelar, A . 

6. Diámetro equiy^^íl^^'tfeWe, D . 



1. Límite elástico del cable. 

2. Tensión de rotura del cable, o R Máx . 

Las variables no controlables son. 

1. Margen de error del método de solución de las ecuaciones.. 
Las variables de salida son. 

1. Esfuerzos o fuerzas en la cuerda, o,T . 

2. Aceleración de la partícula material, y " M . 

3. Velocidad de la partícula material, y ' M . 

4 . Desplazamiento de la partícula material, y M . 

5. Esfuerzos internos en el cable en la zona de impacto, o VM . 
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Figura 12: Representación como caja gris de los modelos 
matemáticos planteados. 


3.5 . Procedimien 

/Sf 

El procedimiento 
primero se establece I 
Luego se formula un procedí 
su solución y para la validación coi 



expuesto antes. De esta manera 
jlificaciones geométricas y físicas, 
establecer el modelo matemático, 
atos experimentales y simulados. 

4# 


3.5.1 . Restricciones, simplificaciones geométricas y 
físicas del fenómeno. 


Este estudio tiene como marco de referencia a la primera escala de 
impacto de la caracterización de Weirauch, es decir, a deformaciones 
plásticas locales (zona de contacto), velocidades de impacto menores a 50 
m/s y ratio de deformaciones de los cuerpos que impactan del orden 10° 
mm/mm/seg [43], 


Para el cable de acero bajo impacto se toma en cuenta las siguientes 
restricciones y simplificaciones: 


1. La geometría del cable se simplifica a la de una cuerda 
unidimensional finita con tensión inicial cero y densidad uniforme. Se 
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asigna el valor del área transversal como una propiedad constante en 
cada punto de la cuerda. 

2. Solo se modela cables hechos con acero, alambres trenzados 
helicoidalmente y con perfil circular sin alma de fibras naturales o 
sintéticas. 

3. Únicamente se toma en cuenta los efectos lineales expresados en la 
curva esfuerzo-deformación lineal del cable y hasta el límite de 
resistencia o de ruptura del cable. 

4 . No se toma en cuenta los efectos de los esfuerzos de torsión, 
cortantes o rozamien^ v pu!^¡^ambres del cable que surgen durante 

longitudinal. El efecto del 



reSpcomodo de estos durante la 
Irá en cuenta usando un módulo 
¡I cable. 


5. Se supone deformaciones pequeñas que solo dependen del tiempo y 

.Ñr 

no de la posición e independientes de su ratio temporal y de los 
cambios térmicos. Además, se considera un material isótropo. 

6. También se supone que la tensión aplicada al cable no provoca un 
cambio de volumen, es decir, el acero sufre un cambio de estado 
isocórico y no se presenta ninguna onda de choque que resulte en 
cambios de fase de los materiales considerados. 

7. El impacto se produce en el punto medio del cable o en la coordenada 
cero de la cuerda. 

Para la masa que impacta se tiene lo siguiente: 

1. Los cuerpos que puedan impactar como rocas, vehículos o 
maquinaria se asumen como puntos materiales y con velocidad dada, 
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aunque en el modelo de los esfuerzos internos debido al contacto si 
se asigna una forma. 

2. Su dirección siempre es recta y transversal al cable al inicio. 


Para las condiciones del movimiento del cable y la masa durante el 
choque se tienen las siguientes restricciones: 

1. No se toma en cuenta ni la gravedad ni la resistencia del aire (u otro 
medio) sobre el cable o masa que impacta 


2. Se supone que el punto de contacto entre la masa y el cable adquiere 


instantáneamente la misma velocidad. 

irj, 


El cable idealizadO;.ifnpacta 
figura 13, la masa del p 
iniciales S 0 , coordena 



pof-^n punto material se presenta en 
elocidad inicial V 00 , coordenadas 
velocidad actual V . 


Anclajes rígidos 

¥— 

Idealización del cable ^ 


Figura 13: Representación geométrica de la idealización de un cable 
impactado por una masa puntual. M, es la masa del cuerpo, V 00 , es su 
velocidad inicial, V, es su velocidad actual, O, es el origen de coordenadas de 
ejes x, y, S, es la coordenada curvilínea actual de la cuerda y S 0 , es la 
coordenada de referencia. L, es la longitud del cable. 


3.5.2. Procedimiento para la formulación de los 
modelos matemáticos. 

1. Primero se define un tiempo continuo t, simétrico en cada momento 
(lo que permite las leyes de conservación de la energía y masa) e 
independiente del espacio o los cuerpos que hay en este. También se 
define un espacio a partir de una referencia ¡nercial con coordenadas 
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cartesianas Oxy , propiedades euclidianas y bi-dimensional, continuo, 
simétrico respecto a las acciones entre los cuerpos dentro de este 
(esto permite la conservación de la cantidad de movimiento) 

2. En el espacio definido se encuentra la idealización de un cable de 
acero como un sólido unidimensional finito, deformable, con densidad 
uniforme, sin resistencia a la flexión, de material elástico, inicialmente 
recto y coincidente con el eje x de la referencia ya definida. 


También se define un punto material de masa finita M , que es la 
idealización de un cuerpo sólido, en el mismo espacio que llega a tocar 
al cable con una velocidad inicial, V nn , en la dirección negativa el eje 
Y hacia el punto rj^aTo dej^cüófcá^que es el punto de contacto. Se 
supone que esta 

/ «f 

durante el movimi 


4 . Luego de la geometría 


para plantear matemáticam 



a tener una gran desaceleración 


s leyes generales de la mecánica 
modelo general, esto se encuentra 
en el marco teórico por lo que queda el planteamiento de las leyes 


particulares. Aquí solo se considera la ley constitutiva que relaciona los 
esfuerzos dentro del cable idealizado y su deformación axial a=Ee o 
su equivalente en la cuerda con la tensión T=AEe que se idealiza 
como la ley de Hooke para la etapa de carga (aumento de la tensión 
en el cable) o desaceleración y para la etapa de descarga (disminución 
de la tensión en el cable) o aceleración de vuelta del cuerpo que 
impacta. 


5. Después se plantean las ecuaciones del movimiento del cuerpo que 
choca, para lo cual se puede hacer simplificaciones matemáticas de 
acuerdo a las condiciones físicas del fenómeno. 


6. Con los resultados de las fuerzas axiales se plantea las ecuaciones 
para los esfuerzos internos en la zona de contacto entre el cable y el 
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objeto que lo ¡mpacta transversalmente. Se usa la teoría de Hertz 
asignando una geometría y material al cuerpo que impacta o bien se lo 
considera rígido. 


3.5.3. Procedimiento para la validación de los 
modelos matemáticos. 

La validación cuantitativa o con los números resultado de los 
experimentos o simulaciones solo se hace para el caso del modelo de los 
esfuerzos axiales o tensiones axiales y del movimiento del cuerpo. 


Mientras que para el caso de los esfuerzos de contacto estos se 
calculan para cada caso y de r aQuerdo a lo que predigan se analiza lo 
observado en forma cualitativa cíalos otros autores en la zona de contacto 
entre el cable real y el cuera 


Además se toma 
Si el cable usado 



intos es un cable nuevo entonces su 


curva fuerza-defe^mación tendrá una no-linealidad inicial que 
corresponde al acomod©üy ¡ rozamiento de los alambres y cuerdas, 
seguido de una parte lineal y luego una pequeña transición a la zona 
plástica hasta su rotura. Acá, como se ha establecido en la parte 3.5.1, 
solo se considera la parte lineal para determinar el módulo elástico 
obviando la primera parte no lineal y plástica pero manteniendo el 
límite de rotura dado o establecido por los autores de los estudios 
experimentales y simulados. 


2. Solo se toma en cuenta el tiempo desde el instante de contacto o 
desde que la tensión se detecta en los sensores o es mayor a cero en 
la simulación y hasta la rotura del cable o si la tensión vuelve a cero. 


Con las consideraciones anteriores se validará el modelo de los 
esfuerzo axiales con los resultados de los experimentos del impacto 
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transversal en cables tipo 3x7 realizados por C. Stolle y las simulaciones que 
usan el cable tipo 1x7 de Tytko et al y P. Gospodarczyk [31], [29], [33], 

4 . Planteamiento de los modelos matemáticos. 

4.1 . Planteamiento del modelo matemático de los 
esfuerzos axiales. 

Para afrontar el problema matemático que representa el solucionar las 
ecuaciones diferenciales (6) y (7) que modelan el movimiento de una cuerda 
se puede proseguir usando varios métodos matemáticos que hay. Siempre 
teniendo en cuenta las restricciones físicas reales del fenómeno en estudio. 
Siguiendo esta línea en 1993 Yu. A. Dem’yanov presentó un método 
asintótico (sin dar ejemplo^^araa&i©p con datos experimentales) basado 
en el método de caracÉérístjcas para resolver problemas de propagación de 
ondas en cuerdas aplicadp^impact^transversal con un punto material [67], 
Para esto hace var¡^^P%)lificacpnes como que la deformación y el 
desplazamiento del punto de la cuerdá donde impacta el punto material 
depende solo del tiempo, como sorpresa a continuación (En la figura 13 se 
puede ver las variables dé ^ófíes): 



s{S 0 t)=Bt m ; B,m >0 


( 11 ) 


y(S 0 =0,t) = y 0 (t) = B 1 t n ;B 1 >0 


( 12 ) 


Pero, si se trata con una cuerda perfectamente elástica con 
deformaciones pequeñas (para un cable de acero y comparado con la 
unidad, 1, esto no limita el tamaño de los desplazamientos) y se usa la 
ecuación de la deformación (11) reemplazándola en la expresión T=EAs, 
se tiene que las ecuaciones del movimiento bi-dimensional de una cuerda 
resultan en: 



(13) 
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d l X _ Po . d \y _ Po 

_ O n * O * _ C n * 


< s „ ^ a t r 


Sn 


EA 


6 ( t 2 


( 14 ) 


Este par de ecuaciones son diferenciales parciales no lineales de 
segundo orden y se resuelve con un cambio de variable, para x o y [68]: 


x,y=w(q);q = \a(2 + m) 2 (S 0 +B 2 Y~4t 


2+m|l/2 


;a = 


Po 

EA 


(15) 


Con lo que las ecuaciones diferenciales tienen la siguiente solución: 


w = B 3 e [filn(,) “ amBj -B 5 


R = 


— mi 2 
2+m 


B,—— 1 /am 


(16) 


Entonces resulta que la solución para el desplazamiento del punto de la 
cuerda impactado y(s 0 ^0;t)=wG& () =0;'tZ es diferente a lo que supone, Yu. 
A. Dem’yanov, en su.-fhétodo asintótico. Aa^más si se ingresa la condición 


de frontera d s xld t t{s 0 
impacto se mueve hacia 



e cumple. Predice que el punto de 
n dirección del eje x durante el 
movimiento lo que se supone no p*^¿!Por lo que tenemos que por un lado 
el método propuesto lleva a contradicciones al usar las ecuaciones del 
movimiento, de una cuerda y ño lleva a una solución donde se pueda 
expresar la desaceleración y posterior cambio de dirección o rebote de la 
partícula que impacta. Este error aparente de la teoría tiene su raíz en la 
idealización geométrica de pasar un sólido a un punto material y de un cable 
a una cuerda unidimensional y se explica en el anexo 1, donde se presenta 
una aparente solución al problema en donde se cumplen con las leyes de 
Newton y el procedimiento matemático pero se llega a una contradicción 
física por que la idealización geométrica lleva a una singularidad matemática 
en el punto de impacto. N. Cristescu, evitó tratar esta discontinuidad 
considerando una función sinusoidal a priori para el desplazamiento de ese 
punto en vez de relacionarlo con la inercia del punto material que impacta. Él 
también lo indica la solución de un problema de este tipo requiere que se 
conozcan de antemano las condiciones de frontera como la velocidad del 
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punto de impacto o los desplazamientos de los extremos y con esto hallar el 
movimiento y estado de tensiones de los puntos interiores de la cuerda [37], 
Por otro lado los experimentos y resultados de simulaciones que se citan acá 
solo dan los datos para las fuerzas en extremos del cable y el movimiento 
del punto material que impacta. Entonces acá se determina primero estas 
condiciones de frontera y así aproximar el estado de esfuerzos en los puntos 
críticos que son los anclajes y el punto de impacto. En suma, el objetivo de 
conocer a profundidad el movimiento de un cable impactado se debe limitar 
a lo medido experimentalmente o simulado y en general a como se da la 
aplicación ingenieril y con esto validar los modelos. Esto quiere decir que lo 
que se busca en específico es el valor de la tensión en los anclajes o 
extremos del cable y para el punto material se busca su desplazamiento, 
velocidad, aceleración^ la fuerza de regccíóh a la que es sometido por parte 


del cable. Luego de 
esfuerzos dentro del 
el modelo de los esfuerzos de con 



uede aproximar el valor de los 
pacto añadiendo lo que predice 
Pero para llegar a estas respuestas 


es necesario observar limero las características físicas propias y comunes 

_ 

de los casos experimentaleé'yíSireiuliááos que son [29], [31], [33]: 


1. La masa del cuerpo que choca es muy grande comparada con la del 
cable: 


Investigación de: 

Tipo 

de 

cable. 

Masa (Kg) del 
cable en el 
experimento o 
simulación: M .. 

cable 

Masa (Kg) 
del cuerpo 
que 

i m pacta: M c 

Ratio 

M /M hl 

c cable 

C. Stolle 

3x7 

1.06E+01 

7.76E+02 

7.34E+01 

Tytko et al. 

1x7 

1.17E+00 

4.00E+03 

3.42E+03 

P. Gospodarczyk 

1x7 

8.78E+00 

1.80E+04 

2.05E+03 


Tabla 1: Comparación de masas del cable con el cuerpo que lo impacta. 


2. El tiempo de reflexión o intervalo de tiempo que le toma al frente de la 
onda de esfuerzos de tensión, en recorrer la distancia desde el punto 
de contacto hasta el anclaje: t R = 0.5L/C n , es muy pequeño en 
comparación con duración del movimiento total (desaceleración y 
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aceleración de vuelta, 2da, 3ra, 4ta etapa del movimiento de una 


cuerda impactada). Esta se muestra en el ratio t t otai/t R de la tabla 2: 


Investigación de: 

Tipo 

de 

cable. 

c,„ 

velocidad 
del sonido 
en el cable 
(m/s). 

Semi- 

Longitud 0.5L 
del cable 
(entre 

anclajes) (m). 

Tiempo de 
reflección 
calculado t R 
(Seg) 

Tiempo total 
del 

movimiento 

U, ( Se g)- 

Ratio 

*tota/*R 

C. Stolle (exp. 
1,2-4) 

3x7 

3.17E+03 

4.31E+00 

1.36E-03 

3.00E-01 

2.21E+02 

Tytko et al. 

1x7 

5.00E+03 

5.00E-01 

1.00E-04 

3.80E-02 

3.80E+02 

P. Gospodarczyk 

1x7 

5.00E+03 

3.75E+00 

7.50E-04 

1.00E-01 

1.33E+02 


Tabla 2: Tiempo de reflexión de las ondas longitudinales en los cables 

estudiados. 


Entonces de acuerdo a los ratios de las tablas 1 y 2 se puede aseverar 
que en los movimientos descri|p,s¡Sl 4 ^de lo siguiente: 


La inercia del 
cuerpo que choca 
transversales n 
y tampoco en el movim 



eña comparada con la inercia del 
e en que la reflexión de las ondas 
en las fuerzas dentro del cable y 
erpo. 


Como el tiempo dte reflexión es mpy pequeño comparado con lo que 


duran las etapas de carga y descarga, t to tai, del cable entonces cada 
tiempo, t R , se puede considerar que las ondas reflejadas uniformisan 
la tensión en el cable. 


El punto 1 permite plantear las ecuaciones considerando una masa 
concentrada del cable en el punto de impacto en vez de una masa 
distribuida debido a que la inercia distribuida no influye significativamente en 
el movimiento del cable. El punto 2 dice que el tiempo que toma en tener una 
tensión mas o menos uniforme a lo largo de la cuerda es tan pequeño que la 
tensión en los extremos es similar o no se diferencia mucho de la tensión en 
el punto de impacto, sin embargo, debido a que en la realidad la forma no es 
triangular entonces en el punto de impacto la fuerza de reacción que actúa 
entre los cuerpos tiene una dirección diferente dependiente del ángulo del 
cable en ese punto. Por tanto, como los experimentos y simulaciones se han 
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hecho con cables impactados en su punto medio y reportan la cinemática del 
cuerpo que impacta así como la tensión en los extremos; entonces los 
modelos se plantean para calcular esas variables. Se empieza con un 



Figura 14: Modelo de cuerdas delimitantes de la cuerda Ideal sometida a Impacto 
transversal en su punto medio. M, es la masa del cuerpo, V 00 , es su velocidad inicial, 
V, es su velocidad actual, O, es el origen de coordenadas de ejes x, y, S, es la 
coordenada curvilínea actual de la cuerda y S 0 , es la coordenada de referencia y 
extremos R iy R 2 . V R , desplazamiento del punto medio de la cuerda celeste. V M , 
desplazamiento del punto medio de la cuerda verde o del punto material. L, es la 
longitud del cable. 

En la figura 14 se^xisa un concepto análogo al de las condiciones de 
frontera pero en este ca^j^^ presenta como unos “movimientos 
delimitantes”. Se que el movimiento del cable ideal impactado se encuentra 
delimitado entre los movimientos de dos cuerdas sin masa distribuida, 
impactadas en su punto medio por el mismo cuerpo y con las siguientes 
características: la primera cuerda de color celeste es tangente a la cuerda 
ideal 3 en sus extremos R lf R 2 de manera que la tensión en estos puntos es 
la misma para ambas cuerdas. Luego se tiene a la cuerda de color verde, 
que no es tangente a la cuerda ideal en ningún punto, pero coincide en sus 
extremos R lf R 2 y en el punto medio o de impacto S 0 = 0 de manera que el 
movimiento de este último punto es igual para ambas. La características de 
estas dos cuerdas delimitantes es que no tienen masa distribuida y la forma 
que adquieren durante el impacto es triangular, además, para tener en 


3 La cuerda ideal puede tener ángulos diferentes a cero en sus extremos ya que no tiene resistencia a 
la flexión. 
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cuenta la pequeña contribución de la inercia del cable se re-define a M 
como la masa combinada del cuerpo, M c , y la masa del cable M Cabte . A lo 
anterior se añade que, por un lado, la cuerda de color celeste tendrá 
desplazamientos pequeños y la cuerda verde tendrá desplazamientos 
grandes para tener en cuenta el movimiento del punto material que 
impacta,sometido a la tensión del punto medio en la cuerda ideal. Esto 
último no quiere decir que la tensión de la cuerda de color verde modelará la 
tensión en el punto medio de la cuerda ideal ya que su forma triangular 
implica que su deformación en el punto medio es mucho mayor a la de la 
cuerda ideal. La cuerda verde solo sirve como referencia para modelar el 


movimiento del cuerpo que im^ejt^ 

,ér u . 

Entonces tenemos qug^e^fnp^ geométrico de la cuerda celeste se 


puede escribir a partir d 
de una masa M sujet 



p, &n este caso, el del movimiento 
lineales como en la figura. 


Nivel de tensión en el resorte izquierdo. 

I 1 Nivel de tensión en el resorte derecho. F 



Figura 15 Analogía del modelo delimitante de la tensión en los extremos con un modelo 
de resortes ideales. T, es la tensión en los resortes, a, es el ángulo de inclinación de los 
resortes, F, es la fuerza de reacción sobre el cuerpo que choca, M, es la masa del 
cuerpo, Y r , es su desplazamiento, x, y son las coordenadas y L, la distancia entre ejes. 

En este caso, la idealización con resortes lineales no toma en cuenta la 
geometría de un resorte, en cambio, el modelo matemático consiste en decir 
que en los anclajes R t ,R 2 y en la masa M se ejerce una sola fuerza T , 
proporcional al desplazamiento del punto de impacto y que solo depende del 
tiempo. Entonces se tiene que la transmisión de la fuerza entre el punto de 
impacto y los anclajes es instantánea. El comportamiento del esfuerzo será 
lineal como el que se muestra en la figura 16, para la carga y la descarga. 
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Figura 16: Ciclo de carga y descarga de un cable 
perfectamente elástico, a, es el esfuerzo y e, es la 
deformación. 

Luego el modelo matemático de este movimiento al considerar 
desplazamientos pequeños es el siguiente: 



uniforme en todo el cable o cuerda. Esta tensión es la que “siente” la masa 
M y que en este caso coincide con la tensión en sus anclajes o extremos. 
Ahora, si examinamos los niveles de tensión en la cuerda impactada 
tenemos lo siguiente: 

1. La fuerza en el cable es una función de la posición y del tiempo: 
T(S 0 ,t) , aquí S 0 e [~L/2,L/2\. 

2. Como los anclajes son perfectamente rígidos entonces en estos 
puntos la tensión se duplica (reflexión en una frontera rígida). 

A continuación se presenta como cambia la tensión y la fuerza en el 
punto material según múltiplos del tiempo de reflexión t R =0.5L/C n con 

C„=4(e7p) : 


41 


Esta obra ha sido publicada bajo la licencia Creative Commons Reconocimiento-No Comercial-Compartir bajo la misma licencia 2.5 Perú. 
Para ver una copia de dicha licencia, visite http://creativecommons.Org/licences/by-nc-sa/2.5/pe/ 

















BIBLIOTECA DIGITAL - DIRECCIÓN DE SISTEMAS DE INFORMÁTICA Y COMUNICACIÓN 






Figura 17: Cambio de la tensión en un cable idealizado como una cuerda 
impactada en su punto medio que toma una forma triangular. M, es la masa del 
cuerpo que choca, a, es el ángulo de inclinación de los resortes, F, es la fuerza 
de reacción sobre el cuerpo que choca, T, es la tensión distribuida en la cuerda, 
Cu, es la velocidad del frente de la onda de esfuerzos longitudinal. El cable está 
idealizado como las cuerdas de "movimientos delimitantes”, sin masa. 
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Según la figura 17, recién en el primer intervalo de tiempo lt R hay 
esfuerzos de tensión en toda la cuerda y simétricos respecto al punto medio, 
pero es recién desde el instante 2t R que la tensión se puede considerar 
uniforme en cada mitad, es decir, esto pasa luego que el ciclo de reflexión de 
la onda se completa o mejor dicho cuando el frente de onda, con nivel cero 
de fuerza, regresa a su origen. Después, esto va ha pasar en 4t R y luego 
en St R y así sucesivamente. Entonces se puede aseverar que a partir del 
tiempo t = 2t R cada intervalo de tiempo At = 2t R , la tensión en el cable se 
puede considerar uniforme en todos sus puntos. Por tanto, la ecuación (19) 
se cumple cada intervalo 2it R ; i = min(r (otoí /r R ), Vi e IN — { 0 } . Pero como 
el tiempo t R equivale ajp^rrer LI2 en este caso se tiene que debido a la 

punto medio la uniformidad sucede 
va a la siguiente expresión: 


reflexión y a que el impac 

/W 

cada una distancia 2 C 


F=2AE 



= AE 


y R 


( 20 ) 


2(;[L/2]) 10 kVJ c„(¡v) 

Otra forma de verlo es que el modelo del movimiento transversal de un 

.. . 


punto material unido a dos resortes como en figura 15 se cumple cada un 
tiempo 2t R . Es decir, cada vez que pasa este intervalo de tiempo se cumple 
lo siguiente: 


T(-L/2,it R )~T(0,it R ) ; Vi € IN-{0} (21) 

Y como el tiempo que toma uniformarse la tensión At = 2t R sigue 

siendo pequeño, entonces se puede considerar que se cumple durante 
todo el movimiento del cable: etapas de carga y descarga. Y como it R , 
abarca todo el intervalo de tiempo que duran estas dos etapas; entonces 
variando i se puede aproximar lo que sucede en todo tiempo t, es decir, 
hacer it R **t y escribir lo siguiente: 


y„ y R 

F=AE^^ = AE R 


Cu ít R 


C„t 


( 22 ) 
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El hecho físico modelado con la geometría anterior también se puede 
traducir a un modelo en lenguaje natural usando términos del método de los 
cambios introducido por Newton: 


“Cada cambio de la fuerza de contacto F , es decir, F+d t F(At) 
provoca un cambio en la tensión de la cuerda, que se traduce en 
T+dj(At) . Este valor viaja desde el punto del impacto a través de la 
cuerda a la velocidad del sonido, C n , llega y se refleja en los extremos, que 
al considerarse perfectamente rígidos, provocan que: se duplique su valor, 
en esos puntos, invierte su sentido de avance pero no la naturaleza del 
esfuerzo (la onda sigue siendo de tensión 4 ) lo que se traduce en que las 
ondas de ambos lados de la cuerda,'en este caso iguales al ser el impacto 

sumen al encontrarse. Pero como 
T en recorrer la cuerda es muy 
luego del análisis de modelo 
s uniforme luego de cada intervalo 


en el punto central, se sup 

/W 

el tiempo que demora 
pequeño entonces s< 
geométrico, que el valor de /, 


At = 2t R y por aproximación en 



momento y como este tiempo sigue 


siendo pequeño se dice J qO&üter ecuación de movimiento es válida cada 
instante t. Además, como no se considera rozamiento interno o algún 
amortiguamiento esto resulta en un movimiento ondulatorio resonante” 


La división por C n t significa que el punto 0, o punto medio o de 
impacto, “ve” que los cambios de tensión viajan una distancia C u t y no 
simplemente Ll 2 . Mientras que la división entre 2 se convierte en la 
multiplicación por 2 de la fuerza lo que se interpreta, geométricamente, como 
el efecto duplicador de la reflexión en un extremo perfectamente rígido es 
equivalente a que una onda recorra la mitad de la longitud del cable dos 


4 La reflexión de una onda se puede ver como si dos ondas de iguales, de tensión o compresión, 
viajasen hacia un punto a la misma velocidad pero en sentido contrario, esto provoca que en el 
punto de encuentro una onda estire (o comprima) desde un lado y la otra también estire (o 
comprima) desde el otro lado del punto lo que duplica la amplitud del esfuerzo en ese punto 
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veces o toda la longitud de la cuerda una vez ya que. Teniendo en cuenta 
que el impacto se da en el punto medio. 

Ahora se puede aplicar las leyes de Newton, para modelar el 
movimiento del punto material, obteniéndose la siguiente ecuación 
diferencial: 


d? y H y R 

M- lZ ^-=-AE- r 


df 


C n t 


(23) 


Pero como se observó en la primera parte con el uso del método de 
Yu. A. Dem’yanov [67] y el trabajo de N. D. Cristescu [37] las ecuaciones no 
son integrables desde el tiempo cero por lo que acá consideramos el 
movimiento desde el tiempo 2 t R yja ecuación queda como sigue: 


(24) 

orden, (24), es del tipo Emden- 
L. Leach, este tipo tiene una forma 


La ecuación diferencial 


Fowler y según K. S. Govinder y 



estándar y es integrable [69]: 


l£>ÜS 


d]y 

d.t 2 


=r(')s 


(25) 


En el caso de la ec. (24), se tiene que: 




AE 


n= 1 


M (L+C n t) ’ " " ^ 

Entonces, sabiendo que si existe una solución, se puede dar una 
solución aproximada análoga a la solución lineal de la ec. 19. Se parte de 
que la solución aproximada debe tener la siguiente forma: 


y R. aprox 

Las condiciones iniciales son: 


= B 3 sen (Vb 4 0 +b scos(Vb 6 0 


(27) 
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; ^( í =°)=' , « < 28 ) 

La masa a considerar es M=M c +M cable . Y una solución aproximada 

del desplazamiento (modelo de pequeños desplazamientos) del punto medio 

de la cuerda celeste delimitante es: 


y r. 


aprox 


= -y^ -j^Jl + CjjtsenL^ 


EA 


(29) 


M -J L+C n t 

La ecuación (29) describe un sistema resonante cuya amplitud va 
creciendo junto al decremento de su frecuencia en forma continua. Para 

comprobar que la solución es aproximada se sustituye en la ecuación (24), lo 

TjV i¿\aoai|, 

que resulta en lo siguiente: 


d*y R . 


aprox. 


U¡y R . 


d t t 


Donde: 



U,= 


j(0= 


L + C n t 


4UÍ 


C u t . 

C 2 

u, ii 

r 2 \ 

2 

4 

L + C n t 1 


(30) 


(31) 


4 AE U 2 U x C n t __ 

—;-rflví J = —-77 eos 

M(L+C n t) 4 (L+C n t) 2 

y R . apr 0X . j(t)+g(t) 


U 2 t 


V E + C ¡j I 


(32) 


(33) 


n ( t \— V »» C nt 
Ul(L+C„t) 


eos 


U 2 t 


Vi+Cjjí , 


(34) 


Si se analiza el valor de las constantes de j(t) y g(t) para tclse 
tiene lo siguiente: L^:M/Ap, como E tiene sus unidades en el orden de 
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10 n [Pa], M en 10 3 [Kg] , A en 10 _4 [m 2 ] , p en 10 3 [Kg/m 3 ], C„ en 10 3 
y V hasta el orden de 10 1 . Entonces se tiene que: 


J'(0>1 ; 


4 AE 


g(í)«° 


(35) 


M(L + C n t ) 

La característica de j(t) , con las dimensiones experimentadas y los 


supuestos hechos es que en cualquier situación su valor tiende a ser 1.1091. 
Por lo que se establece un factor de corrección para y R aprox 5 : 


tcorr-K^m 
Entonces, el valor corre 
cuerda con movimiento d 

y R. corr. 0 Corr. y R .aprox. 



(36) 


= 1.1091 



desplazamiento del punto medio de la 
picamientos pequeños es: 


tPcorr.y R. 


aprox. 


y R.Corr. TCorr. 


00 




EA t 
M <J L+C n t\ 


(37) 


(38) 


Luego la ecuación que modela la tensión o fuerza axial en los anclajes 
o extremos durante el proceso de carga o des-aceleración de la masa que 
impacta es la siguiente: 


T C =AE 


Voo lM_ rj—^- 

<t> c ° rr . 2 y EA V L + C n 


tsen 






(ff 

2 

+ 1-1 

V L+C n t 


(39) 


Esta ecuación es válida desde el inicio del contacto hasta que alcanza 
una tensión máxima: 0<f<í r 


'■C.MAX. ■ 


5 Este cálculo se puede comprobar en los archivos de cálculo adjuntos. 
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Por consiguiente, el modelo del esfuerzo axial (promedio o compartido 
por todos los alambres del cable) en los extremos de la cuerda es: 



(40) 


con 0<t<t 


C.MÁX. 


Las dos últimas ecuaciones representan un sistema resonante que solo 
modela al proceso de carga desde el instante de contacto. Por esto se 
restringe el tiempo hasta que se alcanza una tensión máxima donde, de 
acuerdo al movimiento de^ljr®^^ de desplazamientos pequeños y 
deformación elástica, la masa M, sé<$etiene y empieza el proceso de 
aceleración de vueltáyb descar^Jp|^\tanto, para el siguiente proceso se 
tiene que hacer otro modelo porque las ecuaciones (38) y (39) predicen que 
el periodo y la tensión seguirán creciendo indefinidamente en el tiempo 
debido a la reflexión, lo que en la realidad no pasa. 

'$/ .¿Ér 

Para tener en cuenta la descarga se tiene que revisar, otra vez, el 
concepto de módulo elástico de un cable de acero. Este parámetro cambia 
de acuerdo a los ciclos y nivel de carga y descarga a los que se someta el 
cable, de acuerdo a esto y en general, el módulo elástico durante la 
descarga siempre es mayor al que inicialmente se midió en la carga. Por 
esta dependencia de la carga y con el objetivo de tener propiedades 
constantes y por tanto una diseño predecible, K. Feyrer, recomienda 10 
ciclos de carga y descarga de los cables nuevos antes de medir y dar un 
módulo elástico promedio. Pero todo esto se da en el contexto de 
aplicaciones estáticas o con velocidades de deformación ¿ pequeñas 
e: lO^alO 1 mm/mm/s . Con esto no alcanza para estar siquiera en la 
primera escala de la Caracterización de Weirauch ( ¿ en el orden de 
10 s ) [43], [70], Pero como estos ratios de deformación altos si se 
dan en el caso del impacto transversal en los cables donde, como se 
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comprueba en la verificación de estos modelos más adelante, la deformación 
puede alcanzar valores mayores a e:2mmlmmls . Entonces si 
consideramos esta variable se puede predecir que el efecto de la velocidad 
sobre el cable afectará su módulo elástico de la misma manera que el 
aumento de la velocidad de deformación provoca un aumento en la 
resistencia a la tensión de un acero al carbono. Considerando lo anterior se 
puede deducir que un cable nuevo, o no estabilizado, cambiará su módulo 
elástico en forma continua durante el impacto transversal y alcanzará un 
valor mayor a lo medido estáticamente. 


La descarga de un cable ¡mpactado es un fenómeno que podría ser 
descrito usando los cambfo^vJáélódiagrama de esfuerzo-deformación de la 


figura 16. Pero todo cable, 
permanente cuando es c 
alambres y más si se 
aproximación se da en la 
deformación máxima alcanzada, 
esta etapa el cable presi 



te uno nuevo, tiene una deformación 
ado, debido al re-acomodo de sus 
ción plástica. Por esto, una mejor 
de el cable se descarga desde la 
asta una final, e f , mayor a cero. En 
i elástico mayor porque su estructura 


jn móduli 

DE 

es más estable lo que implica que su módulo elástico será más próximo a un 
valor constante. 



Figura 18: Idealización lineal de los procesos de carga y 
descarga de un cable de acero, a, es el esfuerzo, £ M áx, es la 
deformación máxima y s f , es la deformación final. 

En el fenómeno que se estudia acá, la descarga se produce luego de 

que el objeto que impacta está en un estado de velocidad nula y toda su 
cantidad de movimiento se ha transferido al cable que ha alcanzado una 
tensión máxima. La aceleración en el cuerpo también es máxima y empieza 
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el proceso de disminución de la tensión en la cuerda o proceso de descarga. 
Entonces sucede lo inverso a lo que sucedió desde el instante del contacto 
entre los cuerpos hasta que se alcanza la tensión máxima, es decir, todo el 
esfuerzo o tensión acumulado por la reflexión de las ondas de esfuerzo 
empieza a disminuir de la misma manera que aumentó pero con más 
velocidad, ya que el módulo elástico es mayor. Con geometría, lo dicho se 
da en la figura 19, donde se representa un espacio sin gravedad y sin un 
medio que ofrezca resistencia. Ahí dos cuerdas perfectamente elásticas con 
propiedades uniformes e idénticas en todos sus puntos son sometidas al 
impacto longitudinal de masa puntual con una velocidad inicial Voo lo que 


resulta en ciclos de carga y descarga entre una tensión máxima y mínima 
(inicial), de manera quejando la puefá^A está en proceso de descarga la 
cuerda B está en proceso de carga ^^^nd^ una alcanza su tensión máxima 
la otra alcanza una tensión mínima. ¡T/como las propiedades son uniformes 
todos sus puntos tendrán la misma deformación y tensión. 


Descarga de la cuerda A Carga de la cuerda B 



Figura 19: Ilustración para demostrar la simetría respecto al tiempo de los procesos 
de carga y descarga de dos cuerdas perfectamente elásticas e idénticas con niveles de 
tensión T, iguales (máximos y mínimos) e impactadas longitudinalmente por una 
masa puntual M, t¡, son instantes de tiempo igualmente espaciados, t„, es el periodo, 
g, es la gravedad, x,y, son las coordenadas del espacio, S, es la coordenada curvilínea. 
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Y la interpretación matemática de la figura 19 es: primero la descarga 
de la cuerda A y la carga de la cuerda B se expresan en un mismo 
intervalo de tiempo re[0; íj , donde t n es el tiempo en que la cuerda A o 
B alcanzan su tensión máxima en cualquier punto S 0A durante su descarga 
o la tensión mínima en cualquier punto S 0B durante su carga 
respectivamente, En este caso t n , resulta ser igual al tiempo que toma la 
carga t c Máx , que aquí también es igual al de la descarga t D Máx . Si se divide 
el intervalo tiempo en partes iguales A t=t n /n, luego 
At = f,-t i _ 1 ; l<í<n, Vi,nelN—{0} entonces se puede escribir la tensión 
durante la descarga de la cuerda A en cualquiera de sus puntos S 0A en 
función de la tensión d0>*1a cueríte B , en cualquiera de sus puntos S 0B 
durante la carga para dada tiempo í|>p 



T A ( *2 ) — ^ H (¿n "2 ) > 

"k|£> ü£ \Í\&^ 


(41) 


T a{íí) — T bÍ^h-í) > 

Como se tiene que t—t^j+At; luego t n =t n _¡+iAt y como iAt=t¡ 
entonces t n _ i =t n -t¡, luego la tensión de la cuerda A en proceso de 
descarga en función de la cuerda en proceso de carga se expresa como 
T A (t¡) = T B (t n -t¡) y teniendo en cuenta que el proceso es continuo, en el 
límite se tiene: 


T A (t) = T B (t n -t)=T B (-(t-t n )) (42) 

Es decir, son funciones simétricas respecto al tiempo con los valores de 

la función T B trasladados un tiempo t n para que se den al mismo tiempo 
que la carga. Esto es compatible con el principio de simetría del tiempo de 
Emmy Noether, con el cual se deduce que “lo que fue antes es después” y 
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así se llega a plantear la conservación de la energía que en el caso de las 
cuerdas en carga y descarga viene a ser la energía potencial elástica [71]. 

La tensión en función del tiempo durante la carga y su proceso 
simétrico de descarga de las cuerdas A y B se indica en la figura 20 que 
sigue. 



Figura 20: Representación de la Tensión Vs tiempo de los procesos de carga y 
descarga de las cuerdas y correspondiente a lo expresado en la figura 19. T máx , 
tensión máxima durante la carga o inicial de la descarga, T min . tensión mínima al 
¡nielar la carga o finalizar la descarga, T D tensión de descarga, T c .s¡m., tensión de 
descarga como función simétrica de la carga. 

Con lo desarrollado, se puede expresar la tensión durante el proceso 

de descarga T D (t) de cualquier cuerda ideal elástica y con propiedades 

constantes como la función simétrica del proceso de carga de la misma: 

TD (0 = -^C.Sim. ( ( ^ ^C. Máx) ) (43) 

Esta ecuación es válida en el intervalo t c Máx <t<t c Máx +t D Máx siendo 

el instante t c Máx , donde se termina el proceso de carga o se alcanza el 

valor de tensión máximo. Pero aquí, debido a que las cuerdas son iguales 

también t c Máx = t D Máx . En la figura 20, se indica como la función simétrica se 
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puede trasladar a cualquier instante de tiempo de manera que puede darse 
de manera simultánea en las dos cuerdas usando la expresión ~{t-t CMáx ) 
o se puede representar como continuación de la carga, en forma continua, 
usando la expresión -(t-2t C Máx ). 

Como las cuerdas A y B no tienen masa ni resistencia a la flexión, 

entonces, su estado en cualquier momento es equivalente al de la cuerda 

del “movimiento delimitante” con perfil triangular, desplazamientos pequeños 

y tensión igual a la tensión en los extremos de la cuerda ideal impactada de 

la figura 14. Y así se vuelve al caso del impacto en cables, acá se tiene que 

durante la carga el módulo elástico no es el mismo que el de la descarga 

. 



de descarga), se de hfi *ajrnniir m 



(44) 


Luego de la carga se produce el cambio de módulo elástico, en este 
caso se asume que es discontinuo aunque en la realidad el cambio se de en 
forma progresiva durante el cambio de la tensión. Pero acá se puede 
establecer la relación E D =R 78 E , donde i? 78 >l es una constante. Para 
demostrar que la ecuación (44) se cumple en el caso del impacto transversal 
se tiene que expresar las tensiones en función de los desplazamientos 
máximos que serán diferentes para la carga y descarga ya que se calculan 
con módulos elásticos diferentes: 



(45) 
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Como E d =R 78 E la ecuación (46) se puede expresar como sigue: 


( 46 ) 


T D.Máx.— AE 




V \ 

co 

Qí 

1 1 

yR.Corr. (C. Máx) | 

-^78 

L/2 1 


+ 1-R 


78 


(47) 


con y R Co rr.S( ED . M á X )~ yR.C 0 rr.( C .Máx)^ R 78 

En esta ecuación y R . Corr . S ( ED ) vendría a ser el desplazamiento durante 
la descarga expresado como función simétrica del desplazamiento durante la 
carga calculado con el módulo de elasticidad de descarga de la cuerda, E n , 
y con el modelo de desplazamientos pequeños de la ecuación (38). Y como 

diente la carga (desplazamientos 


y R. Corr.(C . Máx) > el 

pequeños) es una con 
una constante y para que se 



¿o r ,.(c.M áx / L / 2 ) 2 =- R 79 I también es 


lación (44) se tiene que dar: 


(48) 

raíces de las sumas, con los 


Si en esta relación 

primeros términos de una serie de McCIauñn, se llega a probar que para 
cualquier aumento del módulo elástico, en la descarga, la tensión máxima es 
la misma para el final de la carga y el inicio de la descarga, como sigue: 


1 ^79 

l + ^^79 — ^78 + — 2 ^ 78 : 


’ ^79 — ^79 


(49) 


Sabiendo esto se puede expresar el modelo de la tensión en los 
anclajes durante la descarga como la función simétrica respecto al tiempo 
de la carga del cable con su módulo elástico de descarga E D : 


54 


Esta obra ha sido publicada bajo la licencia Creative Commons Reconocimiento-No Comercial-Compartir bajo la misma licencia 2.5 Perú. 
Para ver una copia de dicha licencia, visite http://creativecommons.Org/licences/by-nc-sa/2.5/pe/ 


















BIBLIOTECA DIGITAL - DIRECCIÓN DE SISTEMAS DE INFORMÁTICA Y COMUNICACIÓN 


T D —AE D 


I M 

2 t E d A 


<Í>co ^J—JL + Cnt 


sen 


1 M DESC 


V 


(f) 


+ 1-1 


(50) 


Válida 


en 


el 


» 11 ÍJ Cj JU ^ 

tiempO tDESC = { — {t~{tD.MÁX + tc.MÁX ))) ■> 


tc.MÁx < t^tD.MÁx + tc.MÁx ■ Siendo, t DESC . 


tDESC — { {t ( ÍD.MÁX + tc.MÁX ))) i t C .MÁX<t^t DM Á X +t CM Á X (51) 

En (53) las funciones simétricas del desplazamiento y tensión durante 

la descarga expresadas en función de (-t) se tienen que trasladar, al 


instante t 


C.MAX 


por que a partir de ese momento empieza a disminuir la 

'Si „ 

tensión y el cuerpo empieza su movinftepto de vuelta. Por eso se tiene que 

% 

t _ sí 

iento de la partícula que impacta 
carga del desplazamiento y R , Co rr.D 
correspondiente al movimiento de^^s^lazamientos pequeños “delimitante”. 
De acuerdo a la ecuación (50) el desplazamiento durante la descarga como 
función simétrica de un proceso de carga con módulo E D es: 


añadir el valor, t c 

Antes de pasar al 
también se debe dar el mode 



V 00 j M i 

y R.Corr.S(ED)~ tficorr. 2 J ^ ^ ^ + ^ II t DESC Sen 

Debido a que el módulo E D es siempre mayor en un cable durante la 
descarga el máximo o valor inicial del desplazamiento de la ecuación (52) es 
menor al calculado durante la carga y Jí . Corr . s(£D . Máx) <>' R . Corr . (c . Má;() , Y se 

corresponde con la menor deformación durante la descarga como en la 
figura 18. Sin embargo, esta ecuación (52) no es el modelo del 
desplazamiento del punto medio de la cuerda de movimiento delimitante con 
desplazamientos pequeños. Solo indica cuanto cambia el estiramiento de la 

6 La determinación de estos tiempos, en la práctica, se puede hacer con mayor facilidad con el uso de 
la computadora. 
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cuerda durante la descarga. Entonces para poder hallar el desplazamiento 
actual durante la descarga y que sea una función continua a partir del 
desplazamiento máximo durante la carga se tiene que sumar a este valor los 
cambios de coordenada provocados por la velocidad a partir del instante en 
que empieza la descarga. La velocidad será la derivada de la ecuación (52), 
y ' R Corr s(£D j, luego el desplazamiento en la descarga según el modelo de los 
desplazamientos pequeños: 


y R.Corr.D —y R.Corr.(C.Máx)^ 


f V D d,, 


(53) 


En la ecuación, V D , ^avtóPIMiftcidad durante la descarga y va desde 

.^ 

cero hasta un valor máximo que corresponde a cuando la tensión se vuelve 
cero durante el mismQiprocesIj^La^f^siórTJde esta velocidad es: 


(54) 



Donde: 

¡ \\\ 
y R.Corr.S(ED) — y R.Corr.(ED)\\t— \tD.MÁX + tc.MÁx))) 

Entonces el desplazamiento durante la descarga será: 


(55) 


I L 

J* ^ d — y R .Cor r .S(ED) d t t — y R. Corr. S(ED){t) ~ X R. Corr. S(ED){^C .Máx) (^6) 

1-C.Máx ^C.Máx 

Si se reemplaza el tiempo t c Máx en la ecuación (55) se halla que el 
termino luego del signo menos de la ecuación (56) es igual al máximo que se 
alcanza durante el proceso simétrico de carga con módulo de elasticidad de 
la descarga, E D , como sigue: 

y R.Corr.S(ED){^C.Máx) —y R.Corr.(ED){^D.Máx)~yR.Corr.(ED). MÁX (57) 

Y reemplazando este resultado en la ecuación (53) se halla el 
desplazamiento durante la descarga: 
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y R.Corr.D~ y R.Corr.(C.Máx) ^ R .Corr .(ED). MÁX + ^ R.Corr ,(ED){ (* ( ^D. MÁX + . MÁx) ) ) (58) 

Válido en el tiempo: 


^C.MÁX < ^^^D.MÁX + ^C.MÁX (59) 

De esta manera hay continuidad en el desplazamiento de la partícula 

que impacta según el modelo de desplazamientos pequeños. 


En suma, las ecuaciones (38) y (58) son los modelos de los 
desplazamientos pequeños del punto medio de la cuerda que tiene un 
“movimiento delimitante” de la cuerda ideal impactada en su punto medio y 
con tensión igual a la que esta cuerda tiene en sus extremos. 


4.2. Planteamiento del ^modelo matemático del 
movimiento Oel cuerpo que impacta 
transverialmente un cable en su punto medio. 

Ahora se puede modelar el jn^^niento del cuerpo idealizado como 
punto material que impacta en la cuerda ideal a partir del movimiento 
delimitante de desplazamientos grades representado como la cuerda de 
perfil triangular de color verde de la figura 14. Primero hay que considerar 
que en el caso lineal tiene la siguiente relación: 


F= 2 


AEy, 


y R _ 

con -z-rzr ^ 


L/2 ’ —' Ll 2~ 1 / L V4 + y ! í <60) 

En esta ecuación se asume al ángulo a como pequeño y así, el 
estiramiento de S 0 (L/2,t) , se puede aproximar con el valor de y R pero acá 
ya no se puede hacer esta simplificación por lo que una mejor aproximación 


del estiramiento del punto S 0 (Ll2,t) es V L 2 14+y 2 M -L/2 . Aquí tampoco se 
puede aproximar la tensión en el punto medio de la cuerda ideal con la 
tensión de la cuerda de color verde, de desplazamientos grandes, ya que la 
deformación de esta última es mucho mayor debido a que la forma triangular 
resulta en una deformación mucho mayor que en la idealización curva del 
cable de acero (negro). Entonces, una alternativa para hallar un modelo 
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matemático para el movimiento del punto material se puede usar un método 
aproximado a partir del valor de, y RCorn , es decir, se establece que el 
estiramiento de la cuerda con desplazamientos grandes, y M , se calcula, 
para la carga, y M c , y descarga, y M D , en función de los desplazamientos 
de la cuerda con desplazamientos pequeños: 

V L 2 /4+y 2 M — L12<^ jL 2 /4 + y R Corr — L12 (61) 

Con esto la aceleración del punto Ll 2 o de la masa M durante la 
carga es: 


C ^t3 / M.c _ r, AE I y R.Corr. 

d, t 2 


n 

Para 


-AElyR .Corr. _ yR . Corr. | 

' % 1 


Y la aceleración dure 


d t y m.d 

b ae b 

\Y i 

y R.Corr. D y R.Corr .D 

df ~ 

“ M 

¿fe-, ' 

¿/J J L 2 /A+ y 2 

^ VV ^ J R.Corr.D 


.MAX 


La velocidad de M durante el proceso de carga es: 


i 

d t y M ,C _ .r _ r r-| AE 

d t t ~ V °° J M 


0 




Corr. 


y R . 


Corr. 


1,2 h 2 l^y\. 


d t t 


Corr. J 


Para el intervalo 0<f<í 
Y la velocidad durante la descarga es: 


C .MAX 


(62) 


(63) 


(64) 


d t y m.d. _ 

d t t 


r 2 ae 

y R.Corr.D 3^ R.Corr.D | 

J 1 M 

MÁX 

Ll 2 Jl 2 /4+ v 2 

V ^ 7 J R.Corr .D j 


(65) 


Para el intervalo íc.máy — ^^^d.máy"*'^c.máx 
L a velocidad V F c no necesariamente es cero en este modelo debido 


a que es aproximado y no necesariamente se da al mismo tiempo que la 
fuerza máxima en los anclajes, su valor es: 
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V C.MAX 


=V - f 

<.C. v 00 I 


AE I YR. Corr. 


M 


y R . 


Corr. 


0 


L12 VY/4+/ fi . 


Corr. 


d.t 


( 66 ) 


Para el intervalo 0 <t<t CMÁX 

El desplazamiento de M durante la carga será: 


t 


y M.C.~V 00^ 


-J 


o 


t 


J 


o 


AE 


M 


y R.Corr. 

y R. Corr. j 

L/2 

Vl 2 /4+j4 Corr | 


d.t 


Para el intervalo 0<r<t 


d.t 


(67) 


C.MAX 


El cálculo del desplazamiento en el movimiento de vuelta se hace con: 


t 


y M.D. — V F.C 


-J 


t, 


1 C.MAX 

Para e 

Por último la fuerza de 
impacta en la carga se expresa co 


F C =-2AE 



y r. 


Corr. D 


2 <jL 2 /4+y 2 R ' Corr . D 


I d t t 


D.MÁX + tc.MÁX 


d,í (68) 


F entre el cable y el cuerpo que 


Y en la descarga es: 


F d =-2AE d 


y R. Corr. 

y R. Corr . 

L/2 

Vi 2 /4 + >4. C orr. ( 

intervalo 

0 ^íc.MÁX 

y R.Corr .D 

y R.Corr.D 

L/2 

VY /4+ y R.Corr.¡ 


(69) 


Para el intervalo t 


C.MAX 


— MÁX + t C . MAX 


(70) 


Las ecuaciones (62) a (70) pueden integrarse numéricamente con el 
con el programa GNU Octave o de manera más sencilla usando el programa 
QtiPlot (software de licencia libre para las distribuciones Debían). 

El cambio brusco o discontinuo del módulo elástico que se considera al 
pasar de un estado de carga a uno de descarga provoca una discontinuidad 
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en el cálculo de su aceleración, velocidad y desplazamiento. Aún así, todo lo 
construido se puede usar, directamente, como las condiciones de frontera en 
un método numérico para solucionar las ecuaciones parciales (8) y (9) 
aplicados a una cuerda finita impactada transversal mente en su punto 
medio, con la única restricción de que sea elástica. Aunque implementar el 
método numérico está, por ahora, fuera del alcance del autor. 

4.3 . Planteamiento del modelo matemático de los 
esfuerzos internos debido al contacto. 

Los esfuerzos internos surgen debido al contacto entre el cable y el 
objeto que lo impacta. Ahí aparece la fuerza de reacción, F, que causa 
esfuerzos adicionales en toé*^Yaparte de la tensión lonqitudinal. El 



choca 


cable y así tener un criterlQfj^^íJ^j^hThar un factor de seguridad. 


Para modelar estos esfuerzos se tienen vahas opciones, como: 

1. Tomar en cuenta la geometría completa del cable y el cuerpo que 
impacta: la curvatura de las superficies en contacto, de los alambres 
en el cable, la fuerza entre los alambres y otros constituyentes 
interiores. Pero esto llevaría a tener en cuenta una geometría muy 
elaborada y representaría otro trabajo de investigación. 

2. Otra forma es modelar la superficie del cuerpo que choca y solo tomar 
en cuenta los elementos en la superficie, sin embargo, esto no 
tomaría en cuenta que durante el impacto los alambres y las cuerdas 
de los cables se re-acomodan entre sí interiormente. 
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Entonces la opción mas viable es tener en cuenta el trabajo de G. 
Rebel y R. Verreet, quienes estudiaron el desgaste superficial por la presión 
a la que se someten los cables de acero que se enrollan en un tambor. 
Modelaron los esfuerzos de contacto con la teoría de Hertz tomando en 
cuenta el diámetro nominal del cable de acero, el coeficiente de Poisson 
igual a 0.3 y el módulo elástico del cable (no del acero) [12]. Con este 
método se puede aproximar los esfuerzos a lo largo del eje z (figuras 21 y 
22), ya que son los de mayor valor en la zona alrededor del punto de 
contacto y representan los esfuerzos adicionales que surgen en los alambres 
de un cable de acero impactado. Se considerarán 2 situaciones: Cilindro de 
acero, con diámetro D 1 , en contacto con el equivalente sólido cilindrico del 
cable, con diámetro D 2 , con sy^ejes áy%0° y un cuerpo con frente plano 


rígido, con longitud |f, 
cable, con diámetro D 
y simulaciones con los que s 
equivalente del cable se halla usa 



el equivalente sólido cilindrico del 
rresponden con los experimentos 
parar las predicciones. El diámetro 
área metálica del cable: 


r ^ShJil^Aht (71) 

El caso de dos cilindros perpendiculares se hace suponiendo que el 

área de contacto es circular. Primero se halla un radio equivalente de los 
cilindros en contacto de la siguiente manera [73], [74], [75]: 




(72) 


Donde H 1 y H 2 se calculan usando los diámetros de los cilindros en 
contacto D l y D 2 : 


H 1+ H 2 = |[2/D 1+ 2/D 2 ] 


(73) 


H 2 -H 1 =|[4/Di + 2/ D 2 -8/(DjD 2 )] 


1/2 


(74) 


Con lo anterior se calcula el radio del circulo de contacto: 
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k 1/3 


( 75 ) 


En (75), n es el módulo de corte n-E! 2(l+v). Y los esfuerzos a lo 
largo del eje z son: 


<7 ZZ =-P, 


1 + 


a 


(76) 


°xx — Gyy — p 0 


(l+ v)[ — arctan(a/z)-l]+^(l + ^) 
a ¿a 


(77) 


p 0 , es la presión máxi^igj^el círculo de contacto y depende de la 


fuerza de reacción F j^fre 
solo es necesario usaj'la 


La disposición de los eje 



o que impacta y el cable, en este caso 


(78) 


ver en la figura 21: 



Figura 21: Cilindros de diámetros D h y D 2 con sus ejes a 90 °. x,y,z, son 
los ejes coordenados con origen en el punto de contacto. 


Para el caso de un cilindro de diámetro D en contacto con un cuerpo 
de frente plano rígido se tiene una franja rectangular que se forma en el 
cilindro cuyo ancho, b , es: 
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b=(-^FD[(l-v)/E]) 


1/2 


( 79 ) 


Los esfuerzos a lo largo del eje z son: 


o 


p 0 ((l+2z 2 /b 2 )(l+z 2 /b 2 ) 1/2 —2 z/b) 


(80) 


cr yy =-2p 0 v((l +z 2 /b 2 ) ll2 -zlb) 


(81) 



(82) 


Aquí, p 0 , también es a presión máxima en la superficie de contacto y 
depende de la fuerza de reacción F y de la longitud de contacto del frente 
plano: 



(83) 


En la figura 22 se ilustra la disposición de los ejes y sólidos. 

ñ 



Figura 22: Cilindro de radio R y plano rígido de longitud ñ en contacto.x, y, z, 
son los ejes coordenados con origen en el punto de contacto. 

Con lo desarrollado se puede dar un modelo para el estado de los 
esfuerzos en la zona de contacto durante el impacto. Este se representa en 
la figura 23, acá se incluyen los esfuerzos de compresión, esfuerzos 
principales, o xx ,o yy ,o 2Z calculados con la fuerza de impacto máxima 
F Máx .=Ma Máx • que aparecen debido al contacto que se dan en la superficie 
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cilindrica idealizada del cable y se añade el esfuerzo de tensión máximo 


axial en ese punto ct l , m áx{ s o= 0,t) calculado o medido. 



Figura 23: Estado de esfuerzos principales máximos, a lo largo del eje z 
en el interior de un cable transversalmente impactado.F máx es la fuerza 
de reacción que ejerce el cuerpo que choca, T máx es la tensión máxima 
que alcanza el cable en los anclajes, a xx es el esfuerzo princiapal en la 


dirección x, 


Jyy 


es el esfuerzo princiapal en la dirección y, a zz es el 


esfuerzo princiapal en la dir ección z , o es el esfuerzo debido a T máx . 

Ahora se puede usar algún criterio de falla de materiales como el de 

Von Mises para relacionar los esfuerzp^ internos debido al contacto y el 

esfuerzo axial y así tener unaüSártfBad para compararla con el límite de 


resistencia del material del cable para así determinar un factor de seguridad 


y cuan cerca de la ruptura del cable se está durante el impacto. También se 
debe considerar que no se ha determinado específicamente el valor del 
esfuerzo axial en el punto de impacto. Pero como lo que se necesita saber 
es el valor máximo entonces se lo aproximará con el esfuerzo máximo 


alcanzado en los anclajes. De esta manera el esfuerzo de Von Mises a lo 
largo de eje z con origen en el punto de contacto se determinará de la 
siguiente manera: 


°Vm — y 2 (( & yy Gl.MÁx) + ( ° yy + &L . MÁX &zz) + {&zz &xx) ) (&4) 

Y esto se debe comparar con el esfuerzo de resistencia máximo del 
cable cr m MAX . 
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5 . Validación de los modelos matemáticos. 

Los experimentos con los que se valida el modelo matemático de 
fuerzas axiales son los de C. Stolle [31] y las simulaciones son las de Tytko 
et al y P. Gospodarczyk [29], [33], 


5.1 . Validación de los modelos de los esfuerzos 
axiales y movimiento del cuerpo que impacta 
con los experimentos de C. Stolle. 


En esta parte se aplica el modelo de los esfuerzos axiales según los 
datos de los cables usados así como las condiciones de los experimentos y 
simulaciones hechas. La tensión en los extremos del cable, en el proceso de 
carga y descarga, se halla^Bn las eé&áteiones (39) y (50); la aceleración con 
las ecuaciones (62) y (63); Ap^velocidgdcoleas ecuaciones (64), (65) y (66); 
el desplazamiento con las ecuaciones (67^, (68). 

C. Stolle realizó 


impacto transversal en su punto rrr 



cables de acero sometidos a 
acá se tendrán en cuenta el n.° 1, 2 


y 4 7 . Las características, del cable necesarias para hacer los cálculos 

'y/>, _ 

matemáticos son los de la tabÜ^ lfilQ£- 


Tipo. 

Número de 
cordones. 

Número total de 
alambres. 

Diámetro 

nominal. 

3x7 

3 

21 

19 mm 

Área metálica o 
real del cable. 

1.54508E-04 m 2 

Peso unitario. 

1.22803E+00 Kq/m 

Densidad del 
material 

7.80000E+03 Kq/m 3 

Resistencia a la 
rotura. 

1.80100E+05N 

Límite elástico. 

Módulo elástico en 
prueba estática. 

Módulo elástico de 
cable pretensado. 

Módulo elástico 
recomendado. 

1.10400E+05N 

7.99500E+10 Pa 

1.16000E+11 Pa 

1.04000E+11 Pa 


Tabla 3: Características físicas del cable de acero tipo 3x7 de los 
experimentos de C. Stolle [31]. 

El cable de acero se ancló usando sockets cónicos con resina epóxica 


y se unieron a un soporte de acero que permitía la rotación del anclaje en un 


7 C. Stolle, menciona que en sus experimentos 3, 5 y 6 hubo pérdida de datos y hubo problemas 
técnicos, por lo que no fueron apropiados para su simulación y tampoco lo son para este trabajo. 
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eje vertical. El soporte, a su vez, estaba anclado a una base de concreto 


armado hecha para los experimentos, como en la figura 24. 




Figura 24: Detalle del anclaje del cable 3x7 Figura 25: Sección transversal del cable 
usado para los experimentos de Impacto 3x7 usado en los experimentos de C. 
transversal de C. Stolle [31]. Stolle. 


Para sus experim 
incluyó en su parte fron 


51) para que sea lo único que con 



también construyó un carro que 
de 076mm de acero A36 (figura 
cable durante el movimiento y así 


evitar la interferencia defyresto del carroñe usó un riel que hizo de guía de 
dos de las llantas mientras ér^Frwfúe impulsado por una camioneta hasta 


la velocidad deseada. En la figura 26 se muestra el carro en contacto con el 
cable: 



Figura 26: Posición de impacto entre el carro y el cable 3x7 en los 
experimentos de cables sometidos a impacto transversal de C. Stolle [31]. 


66 


Esta obra ha sido publicada bajo la licencia Creative Commons Reconocimiento-No Comercial-Compartir bajo la misma licencia 2.5 Perú. 
Para ver una copia de dicha licencia, visite http://creativecommons.Org/licences/by-nc-sa/2.5/pe/ 












BIBLIOTECA DIGITAL - DIRECCIÓN DE SISTEMAS DE INFORMÁTICA Y COMUNICACIÓN 


En esta sección se aplica el modelo desarrollado para las fuerzas 
axiales en cuerdas y se comparará con los resultados experimentales y 
simulados de C. Stolle. Es necesario mencionar que en estos experimentos 
la longitud del cable no fue igual a la distancia entre anclajes: en los 
experimentos 1 y 2 la longitud adicional fue de 46.2 mm y en el experimento 
4 fue de 50.8 mm. El cable que se usó fue nuevo y su forma antes del 
impacto tenía curvaturas debido a su enrollamiento y a la gravedad. Además 
C. Stolle, en sus simulaciones si incluyó en sus simulaciones la longitud y 
forma curvada inicial del cable además de incluir factores para modelar el 
rozamiento interno, la resistencia a la flexión, la torsión y los datos de la 
deformación plástica del cablean el modelo matemático planteado acá, no 

r * n — !x J adfqjpnal, ni el rozamiento interno o 


se toma en cuenta esíl^longit 
resistencia la flexión, 
velocidades de deforma 
parámetros a tener e 
dadas en la tabla 3. La velocidad 



ación plástica o endurecimiento por 
a curvada del cable. Los datos y 
lo son las propiedades del cable 
en el instante del impacto, V 00 , 


su masa, M c , y la longitud, L , entre g^ófajes del cable como se presentan 

■^/Sl.' 

en la figura 27, y las tablas 8, 


Anclajes rígidos 


M c = [772.0; 772; 769.3] [ Kg] 

V 00 = [6.625; 6.75; 10.63] [m/s] 
Cable 3x7. 



Figura 27: Idealización de los experimentos 1,2 y 4de C. Stolle. L, es la longitud 
entre anclajes del cable, M c , es la masa del cuerpo que impacta, V 00 , es la 
velocidad del cuerpo en el impacto, x, es el eje coordenado y O, es el punto 
medio del cable y origen de coordenadas. 
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5.1.1 . Modelamiento del experimento 1 de C. Stolle. 

El valor de las variables de entrada y los parámetros para el modelo 
matemático del experimento 1 de C. Stolle se dan en la tablas 4,5 y 6, los 
resultados de los cálculos se presentan en las figuras 28, 29, 30 y 31 y 
tablas 7,8,9 y 10. 


Velocidad de 
impacto. 

Masa del carro. 

-. . . . Tensión inicial del 

Distancia entre . , . , , 

. . cable (cable nuevo sin 

anclajes. . v , ., . 

1 ninguna pre-tension) 

6.62500E+00 m/s 

7.72000E+02 Kg 8.61100E+00m 

0.00000E+00 N 

Tabla 4: Variables de entrada para el modelamiento matemático del 
experimento 1 de C. Stolle [31]. 

Masa del 
cuerpo M c 

7.72000E+02 Kq 

Velocidad de 
impacto del cuerpo 

Veo. 

6.62500E+00 m/s 

Longitud del cable 
idealizado L. 

8.61100E+00m 

Módulo elástico E 

1.04000E+11 Pa 

Área del cable 

Densidad del 

Velocidad del 

Masa del cable 

material del cable. 

sonido C M 

M H 

cable 

1.54508E-04 m 2 

7.80000E+03 Kg/m 3 

3.65148E+03 m/s 

1.05746E+01 Kg 

M= M + M H 

c cable. 

u, 

u 2 

^c.MÁX 

7.82575E+02 Kg 

2.31167E-02 m 

2.86589E+02 m° 5 /s 

1.84057E-01 Seg 


Tabla 5: Parámetros del modelo matemático del experimento 1 de C. Stolle 

durante el proceso de carga. 


Masa del 
cuerpo M c 

Velocidad de 
impacto del cuerpo 

v 00 

Longitud del cable 
idealizado L. 

Módulo elástico E 

7.72000E+02 Kg 

6.62500E+00 m/s 

8.61100E+00m 

1.68411 E+11 Pa 

Área del cable A 

Densidad del 
material del cable. 

Velocidad del 
sonido C M 

Masa del cable 

^cable 

1.54508E-04 m 2 

7.80000E+03 Kg/m 3 

4.64662E+03 m/s 

1.05746E+01 Kg 

M= M + M hl 

c cable. 

u 1 

u 2 

^c.MÁX 

7.82575E+02 Kg 

1.81660E-02 m 

3.64693E+02 m 05 /s 

1.84057E-01 Seg 


^D.MÁX 

^c.MÁX + ^D.MÁX 

1.44601E-01 Seg 

3.28658E-01 Seg 


Tabla 6: Parámetros del modelo matemático del experimento 1 de C. Stolle 

durante el proceso de descarga. 
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Tabla 7; Comparación de los cálculos de la tensión en los anclajes del cable 

impactado del Exp. 1 de C. Stolle. 
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Tabla 9: Comparación de los cálculos de la velocidad del cuerpo 
que impacta del Exp. 1 de C. Stolle. 
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9.13274E+00 g _ 1.02500E+01 g _ 9.49000E+00g 

1.84095E-01 Seg. 1.71030E-01 Seg. 1.93400E-01 Seg. 

Tabla 10: Comparación de los cálculos de la aceleración del 
cuerpo que impacta del Exp. 1 de C. Stolle. 
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5.1.1.1 .Discusión del modelamiento dei experimento 
1 de C. Stolle. 

La problemática que surge al querer usar las ecuaciones desarrolladas 
es determinar el valor de los parámetros que se va a usar. En este caso la 
magnitud de la masa, M , y velocidad de impacto, V 00 , no son problema ya, 
sin embargo, detalles como que en el experimento no se pudo tener un cable 
completamente recto (o lo más posible) requiere que se tome una decisión 
sobre que longitud usar. Aquí se usa la longitud, L , igual a la distancia entre 
anclajes y se calcula la masa del cable con esta medida teniendo en cuenta 
que la diferencia de longitud con el cable real es de menos de 50mm. Otra 
decisión a tomar es sobre eLq^yJ^elástico a usar ya que se usa un cable 
nuevo y se espera que £é$/alor c^hnbie al variar la tensión. En cuanto a esto, 
C. Stolle, menciona que durante la desqarga el módulo será mayor a lo que 

lo de la velocidad del sonido, C n 
ulo a usarse para un cable nuevo 
pa que se mide en la prueba de un 


se da en la carga, agre 
, en su simulación recomiend 



sea 104 Gpa. A diferencia de los 
cable nuevo tipo 3x7. Además, agregare si el cable ha tenido una pre- 
tensión entonces se debe usar“116 Gpa y aclara que estos no son los 
módulos que usa para simular la carga o descarga, pero no precisa que 
valores usa para el caso de la carga o descarga. Entonces la situación acá 
es decidir que valor tomar para el módulo de carga y descarga dado que en 
los modelos planteados el módulo elástico es constante y con este mismo se 
calcula la velocidad del sonido. Por lo que conociendo lo anterior se 
establece el valor de 104 Gpa como módulo promedio para el tiempo de 
carga por que el cable usado en el experimento es nuevo. Para el caso de la 
descarga se sigue este razonamiento: al finalizar la etapa de carga se puede 
considerar que el cable ha sido sometido a una pretensión y su modulo sería 
aproximadamente 116 Gpa como lo recomienda C. Stolle, para un cable pre¬ 
tensado. Pero como en la descarga debe aumentar más este incremento se 
debe dar, aproximadamente con la misma proporción que el aumento 
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durante la carga, es decir, si al aplicar carga para a modo de pretensión 
(con tensión final cero) el módulo elástico cambia de 79.9 Gpa. a 116 Gpa.; 
entonces en la descarga aumentará la misma proporción, es decir pasará de 
116 Gpa. a 168.46 Gpa. y así se usa en las tablas 6 y 13. 


En la figura 28, el modelo matemático predice una tensión superior 
durante la carga hasta poco después del límite elástico de referencia (según 
la prueba estática) durante la etapa de carga, en los primeros 0.12 seg, esto 
sucede porque el módulo escogido es constante y mayor al que inicialmente 
tiene el cable real que es menor a 79.9 Gpa. en el inicio. Por eso se necesita 
una tensión menor para lograr estirarlo y así lo muestran los datos 
experimentales, con un rápj^Ü®IÉé^te de la tensión luego de los 0.05 seg., 

Ax 

- • hasta que su estructura se estabiliza 


luego sobrepasa poco avpocoi 
y la tensión sobrepasa' el 
que influye en esta difi 
los sockets de los anclajes e 
tensión medida al inicio del movi 



e 140 KN. de referencia. Otro factor 
miento del cable en la resina de 
parte de la energía y disminuye la 
Por otra parte en la zona plástica y 
también en parte de la £%parga la aproximación es mejor, pero se aleja al 
aproximarse al límite elástico y sigue esta tendencia hasta que la fuerza es 
cero. Esto se explica, de nuevo, por que el módulo elástico escogido para las 
ecuaciones matemáticas es constante durante el movimiento mientras que 
en la práctica el módulo disminuye mientras disminuye la tensión como 
expresa en la figura 10. 


La aceleración calculada, de la figura 31 muestra como la 
discontinuidad que representa el asumir un cambio instantáneo de 104 Gpa. 
a 168 Gpa. para la descarga supone un salto en la aceleración calculada. 
Antes de este cambio brusco los valores de la línea celeste se pueden 
considerar como una aproximación promedio de todas las aceleraciones 
durante el proceso de carga mientras que en la descarga, como sucede con 
la aproximación de la tensión de los anclajes, acá también es mejor, lo que 
permite predecir un máximo cercano al del medido experimentalmente. El 
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máximo de la aceleración es una cantidad importante al diseñar barreras de 
protección viales porque esa cantidad será la referencia para examinar las 
fuerzas a las que se someten, por ejemplo, los pasajeros de un automóvil 
durante un choque. La velocidad calculada de la figura 30, muestra que la 
aproximación a la medición experimental es muy próxima tanto en la carga y 
descarga e incluso la velocidad final de vuelta. Esto pasa porque, aunque 
hay una discontinuidad grande en la aceleración calculada, la aceleración 
real no es una curva suave, tiene picos de unas centésimas de segundos 
que al momento de integrarse resultan en una curva más suave de la 
velocidad. Y debido a que la aceleración calculada en la etapa de carga es 
como un promedio y se aproxima mejor durante la descarga; resulta en 
valores de la velocidad muy cercanos a lo real. También, por lo anterior, el 
desplazamiento calculado en la ré?Bulta ser prácticamente el mismo 


a lo medido, a excepció 
que se da la velocid 



en general, aunque el tiempo en 
a a lo real, no coincide con el 
, esto pasa por la diferencia entre la 


tiempo calculado para la tensión m 

% 

tensión de los anclajes y la del punto ,^e impacto, algo que los modelos 
matemáticos desarrollados rfeoraan en cuenta. 


En cuanto a la comparación con la simulación de C. Stolle, es evidente 
que esta predice con exactitud tanto la carga como la descarga y la zona 
plástica, la desventaja es que este autor no brinda todos los datos 
necesarios para replicarla en forma independiente. Sin embargo, al observar 
los modelos matemáticos de la aceleración y tensión durante la descarga 
también modelan con exactitud esta etapa por que el cable tiene un módulo 
elástico casi constante. Por lo que en la práctica se recomienda, para tener 
un diseño predecible, estabilizar el cable antes de usarlo. Esto disminuiría la 
capacidad de absorción de energía del cable, pero en las barreras de 
seguridad es mejor usar absorbentes de energía en sus extremos antes que 
permitir que los cables alcancen una deformación plástica cercana a la 
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rotura o hacer una estructura tan grande y robusta que solo se deforme 
elásticamente. 


5.1.2 . Modelamiento del experimento 2 de C. Stolle. 

El valor de las variables de entrada y los parámetros para el modelo 
matemático del experimento 2 de C. Stolle se dan en la tablas 11, 12 y 13. 
Los resultados de los cálculos se presentan en las figuras 32, 33, 34 y 35y 
las tablas 14, 15, 16 y 17. 


Velocidad de 

Masa del carro. 

. . Tensión inicial del 

Distancia entre . . , , . 

. . cable (cable nuevo sin 

anclajes. . v . .. . 

J ninguna pre-tension) 

impacto. 

6.75000E+00 m/s 

7.72000E+02 Kg 8.61100E+00m 

O.OOOOOE+OO N 

Tabla 11: Variables de entrada para el modelamiento matemático del 


experimento 2 de C. Stolle [31], 

/m mi m Mk m 


Masa del cuerpo 

Velocidad de 
impacto del cuerpo 

V oo 

Longitud del cable 
idealizado L. 

Módulo elástico E 

7.72000E+02 Kg 

6.75000E+00 m/s 

8.61100E+00 m 

1.04000E+11 Pa 

Área del cable 

Densidad del 

Velocidad del 

Masa del cable 

material del cable. 

sonido C M 

M hl 

cable 

1.54508E-04 m 2 

7.80000E+03 Kg/m 3 

3.65148E+03 m/s 

1.05746E+01 Kg 

M=M + M .. 

c cable. 

u 1 

u 2 

^C.MÁX 

7.82575E+02 Kg 

2.35529E-02 m 

2.86589E+02 m° 5 /s 

1.83612E-01 Seg 


Tabla 12: Parámetros del modelo matemático del experimento 2 de C. Stolle 

durante el proceso de carga. 


Masa del cuerpo 

M c 

Velocidad de 
impacto del cuerpo 

Voo. 

Longitud del cable 
idealizado L. 

Módulo elástico E 

7.72000E+02 Kg 

6.75000E+00 m/s 

8.61100E+00 m 

1.68410E+11 Pa 

Área del cable 

1.54508E-04 m 2 

Densidad del 
material del cable. 

7.80000E+03 Kg/m 3 

Velocidad del 
sonido C M 

4.64661 E+03 m/s 

Masa del cable 

^cable 

1.05746E+01 Kg 

M= M + M hl 

c cable. 

u, 

u 2 

^C.MÁX 

7.82575E+02 Kg 

1.85088E-02 m 

3.64692E+02 m° 5 /s 

1.83612E-01 Seg 


^D.MÁX 

t +t 

c.MAX l D.MAX 

1.43980E-01 Seg 

3.27592E-01 Seg 


Tabla 13: Parámetros del modelo matemático del experimento 2 de C. Stolle 

durante el proceso de descarga. 


76 


Esta obra ha sido publicada bajo la licencia Creative Commons Reconocimiento-No Comercial-Compartir bajo la misma licencia 2.5 Perú. 
Para ver una copia de dicha licencia, visite http://creativecommons.Org/licences/by-nc-sa/2.5/pe/ 























0.05 0.1 0.15 0.2 0.25 0.3 0.35 


BIBLIOTECA DIGITAL - DIRECCIÓN DE SISTEMAS DE INFORMÁTICA Y COMUNICACIÓN 


o o o o o 

OOVD^rCMOOOOO 
T— T— T— T— T— OOVO^fMO 



LO 

m 

o 


m 

o 


m 

Psl 

d 


cnj 

d 


m 

T— 

d 


m 

o 

d 


o 


(N>0 1 safeioue so| ua uoisuaj. 


0 ) 

c 

03 

Si 

c 

03 


o 

V) 

s- 

03 

> 

lo 

C 

O 


o 

T3 

O 

-{^ 

O 

Q_ 


JU 

-Q 

s 


03 

"O 

lo 

, 0 ) 

O vi 

o o 

c g 

° -i 

co 

o E 


c 

_ ^ 

« ___ 

o 2 


o 

u 

‘43 

'O 


LO 

03 

c 

o 

vi 

c 

cu 

V) 

J3 

c 

o 

V) 


43 vi 
O 

■s -5 

03 £ 
C 
8 


</) 

E ■§ 
Si .§ 

o "O 

2 -S 

V) 

<U 

c 
.o 

"vi 
c 
o> 


o 

<u 

T3 

O 


03 


£3 03 

i2 "o 
u ■+=. 
o LO 

o <J 

I -o 

O « 

.. -t-> 

CN C 

co <u 

¡I 

fio & 

u. oí 


c 

c /3 


'O 

o 

£ 

E 


£ -2 

b <u ■ 

E 

03 

03 

o -O O 

h- o u 

E 

c/d 

03 r- 

■o £ vas 

CN 

CN 

O 

cp 

« o E 

LÚ 

LLJ 

T3 _ 03 

h- 

LO 

O v§ OS 

LO 

LO 

CO 

O) £ 

% a> 

CN 

O) 


co 

CC (/> 

>< 

00 

co 

'S 



1 - _ . 


03 

03 

C 03 O 

2 

W 

:2 = c 

LO 


<n a> a> 

O 

0 

= «a £ 

£¡ 2 “ 
-53 03 

S c X 

| E S 

+ 

LU 

O 

O 

LO 

LO 

cq 

LLJ 

0 

co 

LO 

03 




c 



'§><-> 


03 

03 

w 

0 ) m 

</> -S 

2 

l- C ■ 

LO 

O 

cp 

Tensión 

imulació 

Stolle 

+ 

LU 

O 

O 

00 

0 

cq 

LLJ 

0 

0 

03 

CN 

03 

o w 


T— 

ved « 






2 o 

'3 O 


03 

03 

w 

33 

<13 v<C 

« E 

2 

1 - O 

LO 

O 

cp 

_ 

c <c 

+ 

LLJ 

CN 

CO 

:2 £ 

C O 

LU 

LO 

CO 

LO 

CN 

cq 

03 03 

co 

1 — "O 
■ O 

00 

£ £ 

T— 


^03 




h- 

h- 


Esta obra ha sido publicada bajo la licencia Creative Commons Reconocimiento-No Comercial-Compartir bajo la misma licencia 2.5 Perú. 
Para ver una copia de dicha licencia, visite http://creativecommons.Org/licences/by-nc-sa/2.5/pe/ 


Tabla 14: Comparación de los cálculos de la tensión en los anclajes del cable impactado del 

Exp.2 de C. Stolle. 
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1.99984E-01 Seq. 1.99600E-01 Seq. 2.08100E-01 Seq. 

Tabla 15: Comparación de los cálculos del desplazamiento del 
cuerpo que impacta del Exp.2 de C. Stolle. 
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3.27590E-01 Seg. 3.00000E-01 Seg. 2.81300E-01 Seg. 

Tabla 16: Comparación de los cálculos de la velocidad del cuerpo 
que impacta del Exp.2 de C. Stolle. 
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1.94500E-01 Seg. 1.88190E-01 Seg. 1.84094E-01 Seg. 

Tabla 17: Comparación de los cálculos de la aceleración del 
cuerpo que impacta del Exp.2 de C. Stolle. 
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5.1.2.1 .Discusión del modeiamiento del experimento 
2 de C. Stolle. 

Aquí se usan los mismos módulos elásticos que en el caso anterior al 
tratarse del mismo cable. Aún así, como se ve en la figura 32, la diferencia 
de los niveles de tensión en los anclajes, durante la carga, son mucho 
mayores a los medidos y la diferencia va mucho más halla del valor del límite 
elástico de referencia, 110 KN. Pero también se observa que el valor máximo 
de la tensión es mucho más próximo que en el modeiamiento del 
experimento 1: ahí el valor máximo calculado es 156.62 KN. frente a 165.92 
KN. medidos, aquí el valor calculado es 162.56 KN. frente a 165.5 KN. 
medidos. A pesar que la veloqidad inicial es solo un poco mayor: en 1 es 
6.625 m/s y aquí es 6.75 m/s. Sin embargo, este incremento es el que 
permite calcular un valor 


tener el mismo comporte 
calculados se parecer 


cercano al real. Pero la curva debería 
(figura 28), en cambio, los valores 
k La explicación para todas estas 


diferencias se encuentra en el mayor desplazamiento permanente del cable. 
C. Stolle, notó que apart^ de su longitud,^ólgada mayor a la distancia entre 
los anclajes, se suma la deformación plástica y el desplazamiento adicional 
debido al mayor deslizamiento de los extremos del cable en la resina de los 
sockets usados en los anclajes: 2.4 mm frente a 1.5mm del experimento 1. 
Esto absorbió un mayor porcentaje de energía durante la carga, aumentando 
la diferencia con el modelo de anclajes perfectamente rígidos y amenoró el 
pico de la tensión [31], En la figura 36, se aprecia el deslizamiento del cable 
en los anclajes 8 . 


8 Se debe notar que, a diferencia del módulo elástico variable, el fenómeno de absorción de la onda 
elástica que provoca que la amplitud reflejada no llegue a duplicar la magnitud que llega y el 
desplazamiento adicional en los anclajes por la deformación de algún absorbente de energía si se 
puede incluir en desarrollos posteriores de las ecuaciones aquí planteadas, siguiendo el mismo 
método presentado en 4.1. 
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Figura 36: Deslizamiento final de 2.4mm del cable 3x7 en la resina de los sockets 
usados como anclajes en el experimento 2 de C. Stolle [31]. 

tensión debido a un exceso de 


Además, la mayor^SÍfe|^^pí^ 


desplazamiento permanente del 



también se evidencia en como la 


aceleración calculada rfy^xima, de la figura 35, es incluso mayor que la 
medida y en toda la etapa de'$e&fcáirga. Igual pasa en la mayor parte del 


proceso de carga. La velocidad presentada en la figura 34, es menor en todo 
instante debido a la mayor aceleración. Lo mismo se da con el 
desplazamiento de la figura 33. Y al igual que en el experimento 1 el tiempo 
en que se da la velocidad cero si se aproxima a lo real pero tampoco 
coincide con el tiempo calculado para la tensión máxima. 


La simulación numérica de C. Stolle también presenta una mayor 
diferencia comparado con lo que resulta para su experimento 1. Su 
explicación es que esto sucede debido al mayor estiramiento permanente del 
cable por la deformación plástica y el deslizamiento en los sockets de los 
anclajes [31], 
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5.1.3 . Modelamiento del experimento 4 de C. Stolle. 

El valor de las variables de entrada y los parámetros para el modelo 
matemático del experimento 4 de C. Stolle se dan en la tablas 18 y 19. Los 
resultados de los cálculos se presentan en las figuras 37, 38, 39 y 40 y las 
tablas 20, 21,22 y 23. 


Velocidad de 
impacto. 

Masa del carro. 

Distancia entre 
anclajes. 

Tensión inicial del 
cable (cable nuevo sin 
ninguna pre-tensión) 

1.06300E+01 m/s 

7.69300E+02 Kq 

8.61100E+00 m 

0.00000E+00 N 


Tabla 18: Variables de entrada para el modelamiento matemático del 
experimento 4 de C. Stolle [31]. 




Masa del 
cuerpo M c 

Velocidad de 
impacto del cuerpo 

Voo 

Longitud del cable 
idealizado L. 

Módulo elástico 

E 

7.6930E+02 Kg 

1.05700E+01 m/s 

8.61100E+00 m 

1.04000E+11 Pa 

Área del cable 

Densidad del 
material del cable. 

Velocidad del 
sonido C M 

Masa del cable 

^cable 

1.54508E-04 m 2 

7.80000E+03 Kg/m 3 

3.65148E+03 m/s 

1.05746E+01 Kg 

M= M + M hl 

c cable. 

u, 

u 2 

^C.MÁX 

7.79875E+02 Kg 

3.68184E-02m 

2.87084E+02 m° 5 /s 

7.93582E-02Seg 


Tabla 19: Parámetros del modelo matemático del experimento 4 de C. 
Stolle durante el proceso de carga. 
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7.33976E-02 Seg. 1.04000E-01 Seg. 8.62000E-02Seg. 6.62788E-02 Seg. 

Tabla 20: Comparación de los cálculos de la tensión en los anclajes del cable impactado del 

Exp.4 de C. Stolle. 
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Tabla 21: Comparación de los cálculos del desplazamiento del cuerpo que Impacta del Exp.4 

de C. Stolle. 
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7.04283E+00 g 1.10490E+01g 1.05960E+01g 1.04106E+01g 

7.33976E-02 Seg. 8.79130E-02 Seg. 8.66130E-02Seg. 8.62000E-02 Seg. 

Tabla 23: Comparación de los cálculos de la aceleración del cuerpo que impacta del Exp.4 de C. Stolle. 
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5.1.3.1 .Discusión del modelamiento dei experimento 
4 de C. Stolle. 

Aquí solo hay un proceso de carga hasta la rotura, por eso se ha 
escogido el mismo módulo elástico para esta etapa que en los dos casos 
anteriores, es decir, E= 104 Gpa. La predicción matemática de: la tensión, 
aceleración, velocidad y desplazamiento tiene las mismas características 
que antes, es decir, la tensión calculada es mayor, lo mismo pasa con la 
aceleración lo que lleva a que la velocidad sea menor aunque el 
desplazamiento sea casi igual (hasta el tiempo de rotura predicho). La 
tensión calculada es mayor debido a la deformación permanente y al 
deslizamiento del cable en la resina de los sockets. Por esto el tiempo 
predicho para la rotura, es menor al medido o simulado (que incluye la 
deformación plástica)^ Esta razón llevg á^gue no se pueda predecir los 


valores de la aceleració 


El hecho que resalta acá 



d final o desplazamiento máximo. 


e da la rotura del cable 1x7, esta no 
mbres si no que primero se rompen 


se da al mismo tiempo en todos 
dos cuerdas y con la que queda se ,$éfiene el carro, que acá tiene una 
velocidad inicial de 10.57 m/s. Esta forma de rotura no se puede predecir 
con el modelo matemático de una cuerda unidimensional. En la figura 57 se 
presenta como quedó el cable 3x7 en la zona de contacto luego del impacto. 


Una medida a tener en cuenta es la tensión final con la que sucede la 
rotura de las dos primeras cuerdas, 171.26 KN., mientras que la resistencia 
máxima o de rotura que señala C. Stolle es 181 KN., esta diferencia de solo 
unos ~10 KN. (5.5% del límite rotura) no es algo significativo cuando se trata 
de aplicaciones de ingeniería porque el factor de seguridad usual es siempre 
mayor a 2.5 [76], Además, indica que la suposición de tomar la tensión en 
los anclajes como igual a la que se da en el punto de impacto si se aproxima 
a la realidad y coincide también con el hecho de que la inercia de la masa 
del cable no es una variable determinante durante los procesos de carga (no 
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aporta significativamente a la tensión) 9 . Sin embargo, esta diferencia de 
tensiones si indica una disminución de la capacidad del cable y es la razón 
del planteamiento del modelo de esfuerzos internos debido al contacto. 


5.2. Validación de ios modelos de los esfuerzos 
axiales y movimiento del cuerpo que impacta 
con las simulaciones de Tytko et al y P. 
Gospodarczyk. 

En estos trabajos se usó el cable 1x7, sus características se dan en la 


tabla 24 y en la figura 41: 


Tipo. 

Número de 
cordones. 

Número total de 
alambres. 

Diámetro nominal 

mm 

1x7 

1 

7 

15.7 mm 

Área metálica o 
real del cable. 

Peso unitario. 

Densidad del Acero. 

Límite elástico. 

1.50000E-04 m 2 

1.17200E+00 Kg/m 

7.80000E+03 Kg/m 3 

2.46000E+02 N 


Resistencia a la 
rotura. 

2.79000E+05 N 


Módulo elástico 
estático sin 
pretensión. 

1.95000E+11 Pa 


Tabla 24: Características físicas del cable de acero tipo 1x7 de las 
simulaciones de Tytko et al y de P. Gospodarczyk [33]. 



Figura 41: Sección transversal del cable 1x7. 


9 A excepción de su influencia junto a la resistencia a la flexión de la estructura de alambres que 
provoca que el cable tome una forma curva en vez de un triángulo durante su movimiento. 
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Los autores de ambos trabajos solo indican que usaron el programa 
LS-DYNA (Al igual que C. Stolle) y no dan detalles del tipo de código 
numérico que usaron. Aunque en este caso las simulaciones se hicieron con 
cables rectos y una longitud igual a la distancia entre sus anclajes definidos 
como perfectamente rígidos. 


Primero se modelará la simulación de Tytko et al, los datos y 
parámetros a tener en cuenta son las propiedades del cable dadas en la 
tabla 3, la velocidad del bloque en el instante del impacto, V 00 , su masa, 


M c , y la longitud, L , entre anclajes del cable como se presentan en figura 


42 y las tablas 8, 9 y 7: 


M =4000 Kg 


V oo =3.0m/s 


Cable 1x7. 


◄ ^ 

l // 0 

// / 

P * 


' 

L = 

◄ 

// 

1 m 

► 

X 


Figura 42: Idealización de la simulación presentada por Tytko et al.L, es la 
longitud entre anclajes del cable, M c , es la masa del cuerpo que impacta, V 00 , 
es la velocidad del cuerpo en el impacto, x, es el eje coordenado y O, es el 
punto medio del cable y origen de coordenadas. 

Tytko et al, muestra el instante de rotura del cable en la figura 43: 



Figura 43: Resultado de la simulación de Tytko et al mostrando la 
ruptura del cable 1x7 de lm. 
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A continuación se aplica el modelo desarrollado para las fuerzas 
axiales en cuerdas y se compara con los resultados de la simulación de 
Tytko et al de la tensión en los extremos del cable. El caso de Gospodarczyk 
requiere un análisis diferente es por eso que se da en otra sección. 

5.2.1 . Modelamiento de la simulación de Tytko et al. 

El valor de las variables de entrada y los parámetros del modelo 
matemático de la simulación de Tytko et al, se dan en la tablas 25 y 26, los 
resultados de los cálculos se presentan en la figura 44 y la tabla 27 [29], 


Velocidad de 
impacto. 

Masa del Distancia entre Tensión inicial del 

proyectil. anclajes. cable. 

3.00000E+00 m/s 

4.00000E+03 Kg I.OOOOOE+OOm 

0.00000E+00 N 

Tabla 25: Variables de entrada para el modelamiento matemático de la 

simulación de Tytko et al [29]. 


Jfftm 

mjm 


Masa del 
cuerpo M c 

Velocidad de 

impacto del cuerpo Longitud del cable 
v idealizado L. 

v oo. 

Módulo elástico 

E 

4.00000E+03 Kg 

3.00000E+00 m/s I.OOOOOE+OOm 

1.95000E+11 Pa 

Área del cable 

Densidad del Velocidad del 

material del cable. sonido C M 

Masa del cable 

^cable 

1.50000E-04 m 2 

7.80000E+03 Kg/m 3 5.00000E+03 m/s 

1.17000E+00 Kg 

M= M + M hl 

c cable. 

u, 

U 2 


4.00117E+03Kg 

1.75437E-02 m 1.71001E+02 m 0 5 /s 



Tabla 26: Parámetros del modelo matemático de la simulación de Tytko et 
al durante el proceso de carga hasta su rotura. 
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matemático. simulación de Tytko et al. Matemático. Tytko et al. modelo matemático. 

1.99049E-02 Seg. _ 1.98825E-02Seg. _ 2.79000E+02 N _ 2.89630E+05 N 8.74344E-01 mm/mm/s 

2.46000E+05 N 2.46000E+05 N 2.12365E-02 Seg. 3.80749E-02 Seg. 2.12153E-02 Seg. 

Tabla 27: Comparación de los cálculos de la tensión en el cable de la simulación de Tytko et al. 
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5.2.1.1 .Discusión del modelamiento de ia simulación 
de Tytko et al. 

Aquí solo fue necesario hacer los cálculos para la tensión y se 
comparan con el resultado de la simulación de Tytko et al, en la figura 44. La 
diferencia de las condiciones con los experimentos de los casos anteriores 
se da por que la predicción numérica incluye un cable recto sin holgura y con 
anclajes perfectamente rígidos al igual que lo considerado en el modelo 
matemático. Además, el módulo para la zona de deformación elástica del 
cable es 195 Gpa, que además, es muy cercano al de una barra de acero, 
207 Gpa, lo que quiere decir que no vahará demasiado durante el 


movimiento, a diferencia de .„ 
experimentos de C. Stolle [c 
elástico es muy aproxima 


de referencia, 246 KN 
menor módulo elástio 
predecir una tensión de rotura 



ue sucede con el cable 3x7 de los 

esto lo calculado hasta el límite 

A 

tiempos. Después del límite elástico 
empieza a prolongarse, por el 
en el tiempo hasta que llega a 
KN. Esta tensión de rotura mayor a 


la medida estáticamente de 279 KÑ, s^y-debe a que en la simulación se 
incluyó el efecto de endureóiftteiJao del material. Por otro lado si se incluyese 
en el modelo el comportamiento del cable en la parte plástica con un módulo 
menor al elástico haría que la curva entre el límite elástico, 246 KN y el de 
rotura (medido estáticamente), 279 KN, se allane y se prolongue 
aproximando mejor la simulación. También se debe observar como es que 
los valores si coinciden a pesar que el cuerpo que impacta tiene un frente 
plano y no puntual como se considera en el modelo matemático aunque el 
hecho predominante es la gran diferencia entre las masas del cuerpo y del 
cable impactado. 
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5.2.2. Modelamiento de la simulación de 
Gospodarczyk. 

Los datos de una simulación que realizó Gospodarczyk, y que se 
deben tener en cuenta para el modelo son los siguientes: 


Velocidad de 

Masa del 

Distancia entre 

Tensión inicial del 

impacto. 

proyectil. 

anclajes. 

cable. 

7.50000E+01 m/s 

1.80000E+04 Kg 

7.50000E+00 m 

0.00000E+00 N 

Longitud adicional 
del cable, m. 

Contacto del cable 
con la masa. 

Forma inicial del 
cable. 

Número total de 
cables de la 
plataforma. 

O.OOOOOE+OOm 

Frente plano de un 
paralelepípedo (Ver 
gráfico). 

Recta 

420 


Tabla 28: Variables de entrada para el modelamiento matemático de la 

simulación de P. Gospodarczyk [33]. 

/«§ 

En la figura 45 se 
bloque de 18000 Kg 


placa metálica (de espesor p 


y 420 cables 1x7 entrecruzados a 



n Inicial de la simulación donde un 
una plataforma que incluye una 
utor no indica cuanto exactamente) 


un patrón de 2(13*6) + 4 (11*6). El 




bloque tiene una base plar%/£|g Z5rj£.£ 1.5m y golpea en el centro. 



Figura 45: Posición inicial de la simulación hecha por P. Gospodarczyk 
de una base artificial de 07.5m para minas de pozo profundo. 1 es el 
bloque que cae, 2 son los 420 cables y 3 una placa de acero [33]. 


La primera parte de la idealización de esta simulación consiste en 


suponer que solo los cables directamente bajo el paralelepípedo soportan el 
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impacto y que la resistencia de la placa es despreciable. Así se muestra en 
la figura 46: 



Figura 46: Primera idealización de la simulación de una base de seguridad 
de minas de pozo profundo p resen tada por Gospodarczyk. 

Para poder aplicar el modelo de los esfuerzos longitudinales se tiene 

que hacer una segunda^ealización dónete el peso total del paralelepípedo 

_ 

se divide entre los (13 + 8)=126 cables 1x7, es decir, 18000 /126 = 142.86 
Kg, e impacta sobre un solo cable 1x7 a una velocidad de 75 m/s, como se 
ilustra en la figura 47. 


V 00 = 75.0 m/s 


Anclajes rígidos 


n M c =18000 Kg /126 cables 
M c =142.86 Kg 

Cable 1x7. 


j/ o 

/ 

£ ^ 

L = 

◄— 

y 

7.5 m 

—► 

" X 


Figura 47: Segunda idealización de la simulación presentada por Gospodarczyk 
de una base de seguridad de minas de pozo profundo.L, es la longitud entre 
anclajes del cable, M c , es la masa del cuerpo que impacta, V 00 , es la velocidad del 
cuerpo en el impacto, x, es el eje coordenado y O, es el punto medio del cable y 
origen de coordenadas. 
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Para calcular la tensión se usa la ecuación (39), como en los demás 
casos, pero para la velocidad se usa la derivada de la ecuación (29) que 
modela los desplazamientos pequeños. Las razones para no usar el modelo 
de desplazamientos grandes se da en la discusión de este caso en 5.2.2.1. 
El valor de las variables de entrada y los parámetros del modelo matemático 
de la simulación de Gospodarczyk se dan en las tablas 28 y 29, los 
resultados de los cálculos se presentan en las figuras 48, 49 y la tabla 30. 


Masa del 
cuerpo M c 

Velocidad de 
impacto del cuerpo 

Voo. 

Longitud del cable 
idealizado L. 

Módulo elástico 

E 

1.42860E+02 Kg 

7.50000E+01 m/s 

7.50000E+00 m 

1.95000E+11 Pa 

Área del cable 

Densidad del 
material del cable. 

Velocidad del 
sonido C M 

Masa del cable 

^cable 

1.50000E-04 m 2 

7.80000E+03 Kg/m 3 

5.00000E+03 m/s 

8.77500E+00 Kg 

Masa en 
movimiento M 

+ M 

cable. 

u, 

u 2 

Uáx@278.6 KN 

1.51635E+02 Kg 

8.53823E-02 m 

8.78402E+02 m° 5 /s 

7.33144E-03 Seg 


tc.MA x@278.75 KN 

7.33462E-03 Seg 

Tabla 29: Parámetros del modelo matemático de la simulación de 
Gospodarczyk durante el proceso de carga de los cables de la base de 

seguridad hasta su rotura. 
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_ 51.29 m/s _ 48 m/s _ 

0.0073388 Seg. 

Tabla 30: Comparación del cálculo de la velocidad final del bloque de 18 TN 
después de pasar por la base de seguridad simulada por Gospodarczyk. 
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5.2.2.1 .Discusión del modelamiento de ia simulación 
de Gospodarczyk. 

Como acá se usó el mismo cable, tipo 1x7, que en el caso de Tytko et 
al, el módulo elástico escogido también es 195 Gpa, la resistencia última o 
de rotura también se asume como la misma, 279 KN. El único objetivo es 
hallar la velocidad final del bloque luego de traspasar una plataforma 
compuesta por un entramado de 420 cables. Para esto se reduce e idealiza 
la simulación para calcular primero el tiempo en que se detiene o se llega a 
romper, como predice la simulación y como también predice el modelo 
matemático. Pero para calcular la velocidad final no solo es cuestión de usar 
los modelos desarrollados y aplicar las fórmulas, se debe tener en cuenta 
que la velocidad del blc^Ére es 7|j m/s"y^omparar sus dimensiones con la 
longitud de los cables* 


El caso es que 
tercera escala de la caracteriza 



idad que se corresponde con la 
eirauch, es decir, es seguro que se 
o se ve en el caso de Tytko et al, la 


van a dar deformaciones plástica 
deformación plástica hasta la rotura, d^la figura 44, se da sin un aumento 
considerable de la tensión. Adémás, se debe considerar que el bloque 
abarca buena parte de la longitud de 7.5m de los cables de la plataforma: el 
lado de 1.5m abarca un 20% y el lado de 2.5m un 33.3%. Esto implica que 
durante el impacto los puntos de los cables directamente bajo el bloque solo 
tendrán deformaciones en dirección horizontal y la pendiente de su tangente 
siempre será cero o paralela a la superficie plana del paralelepípedo. Por 
tanto, la respuesta a la cuestión básica de porqué los cálculos con el modelo 
de desplazamientos pequeños resultan ser tan próximos a los resultados de 
la simulación, mientras que si se usa la aproximación de desplazamientos 
grandes con las ecuaciones (62) a (65), resulta que la velocidad final es 
apenas menor a la inicial. La diferencia está en como se deforma una cuerda 
elástica cuando es impactada por un punto material. En esta situación, la 
inercia local, aunque pequeña, y por que la deformación se difunde con la 


99 


Esta obra ha sido publicada bajo la licencia Creative Commons Reconocimiento-No Comercial-Compartir bajo la misma licencia 2.5 Perú. 
Para ver una copia de dicha licencia, visite http://creativecommons.Org/licences/by-nc-sa/2.5/pe/ 







BIBLIOTECA DIGITAL - DIRECCIÓN DE SISTEMAS DE INFORMÁTICA Y COMUNICACIÓN 


onda de esfuerzos; provocan que las pendientes de los puntos de la cuerda 
sean mayores alrededor del punto de impacto y menores cerca de los 
anclajes, esta forma se puede observar en las figuras 11 y 14. También lo 
muestra N. Cristescu cuando trata el impacto en una cuerda con un extremo 
libre y C. Stolle, en las fotografías de sus experimentos con cables 
impactados transversalmente [31], [37], [53], Mientras que cuando impacta 
un paralelepípedo, una buena parte de los cables tiene una pendiente cero 
y el mismo desplazamiento vertical que el punto medio. De esta manera, los 
estiramientos de los puntos alrededor del punto de impacto no contribuyen a 
su desplazamiento vertical. Por esta razón el modelo desplazamientos 
pequeños llega a predecir una velocidad final calculada de 51.29 m/s, muy 
próxima a la velocidad final simulada de 48 m/s. Y la diferencia se produce 


porque en el modelo d 
endurecimiento del ma 
considera la deforma^ 
pequeño espesor, en 



superficie 


ent'<$s pequeños no se incluye el 
de la deformación y tampoco se 
cto de la plancha de metal, de 
lataforma la cual quedó destruida al 


final como se ve en la figura 50. 

;rve iWJ 



Figura 50: Situación final de la plataforma simulada por Cospodarczyk, 
con el bloque desacelerando desde 75m/s hasta 48 m/s [33]. 
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5.3. Comparación cualitativa de la aplicación del 
modelo de los esfuerzos internos debido al 
contacto. 

Los trabajos a tener en cuena aquí son los experimentos de C. Stolle y 
la simulación de P. Gospodarczyk [31], [33] 10 . 

En esta parte se aplica el modelo de los esfuerzos internos debido al 
contacto según los datos y condiciones de los cables usados en los 
experimentos y la simulación de los trabajos citados. La tensión en el punto 
o zona de impacto o de contacto se tomará como igual a la tensión en los 
extremos del cable, en este caso solo se usará el valor máximo, T Máx , junto 



predicción de las ecuaciones de la aceleración en la carga o descarga (69) y 
(70), el diámetro equivalente del cable se halla con la ecuación (72), los 
esfuerzos internos a lo largo del eje transversal entre dos cilindros con sus 
ejes perpendiculares se halla con las ecuaciones (76), (77) y (78). Los 
esfuerzos internos en el caso del impacto de una masa rígida con superficie 
plana y un cable se halla con las ecuaciones (80), (81), (82) y (83). Por 
último el esfuerzo equivalente que combina los esfuerzos internos por el 
contacto y el esfuerzo longitudinal se calcula con el criterio de Von Mises de 
la ecuación (84). 


10 


No se incluye la simulación de Tytko et al debido a que no indica las medidas del bloque que 
impacta. 
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5.3.1 . Aplicación del modelo de esfuerzos internos 
por contacto al experimento 1 de C. Stolle. 

Los parámetros y condiciones de los experimentos de C. Stolle se dan 
en 5.1. El frente de contacto del carro de sus experimentos fue un perfil 
cilindrico de o76mm de acero A36 dispuesto como en la figura 51 [31]. 



Figura 51: Frente de impacto del carro usado por C. Stolle en sus 
experimentos [31]. 


Entonces para la aplicación del modelo solo es necesario el valor de la 
fuerza de contacto y tensión máximas de cada caso. Los datos del cable se 
dan en la tabla 3 y la figura 25. Se usa el mismo módulo de elasticidad que 
el considerado para la etapa de carga en la validación de los modelos 
anteriores. El esfuerzo de rotura, o m MAX , se halla dividiendo la tensión de 
rotura, 181 KN, con el área metálica del cable. 

El valor de las variables de entrada para el modelo de esfuerzos 
internos debido al contacto del experimento 1 de C. Stolle se dan en la tabla 
31, el resultado de los cálculos se presentan en la figura 52 y la tabla 32. 
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Área metálica del cable. 

1.5451E-04 m 2 

Diámetro equivalente del cable de acero. 

1.4026E-02 m 2 

Diámetro del cilindro que impacta. 

7.6000E-02 m 2 

Módulo de Poisson del cable equivalente. 

3.0000E-01 

Módulo de poisson del cilindro que impacta. 

3.0000E-01 

Módulo elástico del cable. 

1.0400E+11 Pa 

Módulo elástico del acero A36. 

2.0700E+11 Pa 

Tensión máxima medida: T c . 

Expl. Max. 

6.7395E+09 N 

Tensión máxima calculada: T . . ... 

mod. Mat.Expl. Max 

6.6538E+09 N 


Tabla 31: Variables de entrada para el modelo de esfuerzos 
internos debido al contacto del Exp. 1 de C. Stolle. 


Tensión en el cable: Medida (Med.) y calculada 
(Cal.) 

T . =165.2 

Med. 

T a ,c. =156.62 KN 

Rádio equivalente de los cilindros en contacto. 

1.6325E-02 

1.6325E-02 m 

Rádio del área de contacto: a. 

2.2719E-03 

2.2430E-03 m 

Presión de contacto P 0 . 

6.7395E+09 

6.6538E+09 

Fuerza de contacto máxima medida en el Exp.2 

F =M*a 

max.Exp2 Máx.Exp2 

7.2855E+04 

7.0113E+04 N 

Esfuerzo longitudinal máximo o L Máx 

1.0692E+09 

1.0137E+09 Pa 

Esfuerzo de Von- Mises superficial a SMÍn 

2.0979E+09 

2.0365E+09 Pa 

Esfuerzo de Von- Mises máximo a,,.,. 

4.7992E+09 

4.7082E+09 Pa 

Esfuerzo de Von- Mises mínimo a.,... 

V.Min 

1.0869E+09 

1.0307E+09 Pa 

Resistencia máxima del cable ct„ ... 

R.Max 

1.1715E+09 

1.1715E+09 Pa 

% del material del cable con esfuerzos mayores a 
la resistencia máxima. 

43.00 % 

35.00 % 


Tabla 32: Resultado de los cálculos con los datos del Exp. 1 de C. Stolle 
usando el modelo del contacto entre dos cilindricos perpendiculares.. 
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Figura 52: Modelo matemático de los esfuerzos internos en el cable del Exp. 1 de C. Stolle. Las tensiones señaladas son las máximas 
calculada y medida en el experimento, Or. M áx es el esfuerzo de rotura del cable calculado con la tensión de rotura nominal y su área metálica. 
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5.3.2 . Aplicación del modelo de esfuerzos internos 
por contacto al experimento 2 de C. Stolle. 

El valor de las variables de entrada para el modelo de esfuerzos 
internos debido al contacto del experimento 2 de C. Stolle se dan en la tabla 
33 y el resultado de los cálculos se presentan en la figura 53 y la tabla 34. 


Área metálica del cable. 

1.5451E-04 m 2 

Diámetro equivalente del cable de acero. 

1.4026E-02 m 2 

Diámetro del cilindro que impacta. 

7.6000E-02 m 2 

Módulo de Poisson del cable equivalente. 

3.0000E-01 

Módulo de poisson del cilindro que impacta. 3.0000E-01 

Módulo elástico del cable. 

1.0400E+11 Pa 

Módulo elástico del acero A36. 

2.0700E+11 Pa 

Tensión máxima medida: T c , 

Exp2. Max. 

1.6550E+05 N 

Tensión máxima calculada: T . M , ... 

mod. Mat.Exp2. Max 

1.6256E+05 N 

Tabla 33: Variables de entrada para el modelo de esfuerzos 
internos debido al contacto del Exp.2 de C. Stolle. 

Tensión en el cable: Medida (Med.) y calculada 
(Cal.) 

T Me , =165.5 T ca|c =162.56 KN 

Rádio equivalente de los cilindros en contacto. 

1.6325E-02 1.6325E-02 m 

Rádio del área de contacto: a. 

2.2885E-03 2.2849E-03m 

Presión de contacto P 0 . 

6.7888E+09 6.7780E+09Pa 

Fuerza de contacto máxima medida en el Exp.2 

F =M*a 

máx.Exp2 Máx.Exp2 

7.4467E+04 7.4113E+04N 

Esfuerzo longitudinal máximo o L Máx 

1.0711E+09 1.0521E+09 Pa 

Esfuerzo de Von- Mises superficial cr s Mfn 

2.1084E+09 2.0907E+09 Pa 

Esfuerzo de Von- Mises máximo <x, ... 

4.8347E+09 4.8157E+09Pa 

Esfuerzo de Von- Mises mínimo ct„ ,, 

V.Min 

1.0893E+09 1.0702E+09 Pa 

Resistencia máxima del cable <r„... 

R.Max 

1.1715E+09 1.1715E+09 Pa 

% del material del cable con esfuerzos mayores a 
la resistencia máxima. 

44.00 % 40.60 % 

Tabla 34: Resultado de los cálculos con los datos del Exp. 2 de C. Stolle 


usando el modelo del contacto entre dos cilindricos perpendiculares. 
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Figura 53: Modelo matemático de los esfuerzos internos en el cable del Exp. 2 de C. Stolle.Las tensiones señaladas son las máximas 
calculada y medida en el experimento, íTr.móx es el esfuerzo de rotura del cable calculado con la tensión de rotura nominal y su área metálica. 
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5.3.2.1 .Discusión de la aplicación del modelo de 
esfuerzos internos debido al contacto a los 
experimentos 1 y 2 de C. Stolle. 

Las predicciones del modelo son similares para ambos experimentos . 
Cuando se usan los datos medidos resulta que un 43%, en el Exp. 1, y un 
44%, en el Exp. 2, de la sección transversal del cable, en la zona del 


impacto, debería estar rota. Mientras que al usar los datos calculados con 
los otros modelos resulta que el porcentaje disminuye a un 35% y 40% 
respectivamente. Sin embargo, nada de esto sucedió en ambos 
experimentos, como se observa en las figuras 54 y 55. Los cables 



terminaron con una parte encorvada alrededor de la zona de contacto. 
Además, ninguna de su^Oerdas^cacraí jjna de las cuales es el 33.33% del 
cable, se rompió. Es^fecir, el i^rojaelo no^redice lo que sucedió en estos 
experimentos. if\8llk 


Figura 55: Deformación del 
cable en la zona del 
impacto al final del Exp. 1 
de C. Stolle [31]. 


Figura 54: Deformación del 
cable en la zona del 
impacto al final del Exp. 2 
de C. Stolle [31]. 


Esto pasa porque no se tiene en cuenta el deslizamiento relativo entre 


los alambres y las mismas cuerdas que no se comportan como una barra 
sólida al ser aplastadas por una fuerza transversal. Los esfuerzos no se 
distribuyen como en una sección transversal continua porque la sección 
transversal del cable cambia durante el impacto. Si se observa el arreglo de 
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los alambres del cable 3x7 en la figura 25 los tres cordones, cada uno, solo 
se apoya en la superficie de dos de sus alambres dejando sin restricción su 
movimiento lateral, al menos en que la espiral de la cuerda completa un 
ciclo 11 . Esto también significa que no hay un camino estructural continuo por 
donde los esfuerzos se transmitan a toda la sección transversal al mismo 
tiempo, como pasa en mayor medida en, por ejemplo, un cable 18x7 o en 
cualquier otro con más cuerdas y sin alma sintética o natural. 


5.3.3 . Aplicación del modelo de esfuerzos internos 
por contacto al experimento 4 de C. Stoile. 

El valor de las variables de entrada para el modelamiento de los 


esfuerzos internos debida 
tabla 35, el resultado ep lo 

/«SST 

Como en el experiment 
máxima usada coinci 
autor, así se usó en la aplicá 5 ©* 

Aquí, también, se ha escogido 

% 



exp. 4 de C. Stoile se dan en la 
presentan en la figura 56 y tabla 36. 
el cable se rompe la resistencia 
de rotura medida por el mismo 
modelos anteriores y se usa acá. 

smo módulo elástico que para los 

á#" 

Vi* , 


experimentos 1 y 2 y tambí%,gl mismo esfuerzo último de resistencia. 


Diámetro equivalente del cable de acero. 

1.4026E-02 

m 2 

Diámetro del cilindro que impacta. 

7.6000E-02 

m 2 

Módulo de Poisson del cable equivalente. 

3.0000E-01 


Módulo de poisson del cilindro que impacta. 

3.0000E-01 


Módulo elástico del cable. 

1.0400E+11 

Pa 

Módulo elástico del acero A36. 

2.0700E+11 

Pa 

Tensión máxima medida en el Exp4: T c . ... 

” Exp4. Max. 

1.7126E+05 

N 

Tensión de resistencia máxima del cable medida 



estáticamente según C. Stoile: T . . _ . ... 

1.8100E+05 

N 

Tensión máxima alternativa: T . .... 

max. Alt-1 

1.7490E+05 

N 

Tensión máxima calculada: T . .... 

max. Alt-2 

1.7126E+05 

N 


Tabla 35: Variables de entrada para el modelo de esfuerzos 
internos debido al contacto del Exp.4 de C. Stoile. 


11 El movimiento relativo de las cuerdas dentro del cable por una fuerza transversal también se da en 
otros casos que resultan en el aplastamiento de la sección transversal como cuando un cable pasa 
por los labios de una polea o. en general, debido al contacto con cualquier objeto de frente angosto 
[77]. 
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Tensiones máximas en el cable: Medida en el 

Exp.4 (Med.) y según la prueba estática (Cal.) 

T„ ri =171.26 

Med. 

T calc. = 181 KN 

Rádio equivalente de los cilindros en contacto. 

1.6325E-02 

1.6325E-02 m 

Rádio del área de contacto: a. 

2.3542E-03 

2.0545E-03 m 

Presión de contacto P 0 . 

6.9837E+09 

6.0947E+09 Pa 

Fuerza de contacto máxima medida en el Exp.2 

F =M*a 

máx.Exp2 Máx.Exp2 

1.7126E+05 

1.8100E+05 N 

Esfuerzo longitudinal máximo o L Máx 

1.1084E+09 

1.1715E+09 Pa 

Esf. de Von- Mises superficial a s Mín 

2.1743E+09 

2.0703E+09 Pa 

Esf. de Von- Mises máximo 

4.9478E+09 

4.4811E+09 Pa 

Esf. de Von- Mises mínimo 

V.Min 

1.1282E+09 

1.1845E+09 Pa 

Resistencia máx. del cable ... 

R.Max 

1.1715E+09 

1.1715E+09 Pa 

% del material del cable con esfuerzos mayores a 
la resistencia máxima. 

58.00 % 

100.00 % 


Tabla 36: Resultado de los cálculos con los datos del Exp. 4 de C. Stolle 
usando el modelo del contacto entre dos cilindricos perpendiculares. 

Iculos usando los datos del modelo 

^alternativas que vendrían a ser 

las cuales el cable alcanza un 

ligeramente mayor, 


En la tabla 37, se 

j 

matemático y dos ten 
tensiones alternativas 
nivel de esfuerzos igual 



medido 


Tensiones máximas alternativas en el cable: 

T =174 9 T 

náx. Ait-2. 171.26 K N 

Rádio equivalente de los cilindros en contacto. 

1.6325E-02 

1.6325E-02 m 

Rádio del área de contacto: a. 

2.0199E-03 

1.9988E-03 m 

Presión de contacto P 0 . 

5.9921E+09 

5.9293E+09 Pa 

Fuerza de contacto máxima medida en el Exp.2 
F _ M*a 

máx.Exp2 iVláx.Exp2 

1.7490E+05 

1.7126E+05N 

Esfuerzo longitudinal máximo ct l Máx 

1.1320E+09 

1.1084E+09 Pa 

Esf. de Von- Mises superficial cr SMÍn 

2.0185E+09 

1.9873E+09 Pa 

Esf. de Von- Mises máximo or . 

V.Max 

4.3916E+09 

4.3369E+09 Pa 

Esf. de Von- Mises mínimo 

V. Min 

1.1444E+09 

1.1204E+09 Pa 

Resistencia máx. del cable cr ... 

P^.Max 

1.1715E+09 

1.1715E+09 Pa 

% del material del cable con esfuerzos mayores 
a la resistencia máxima. 

58.00 % 

45.00 % 


Tabla 37: Resultado de los cálculos con los datos del Exp. 4 de C. Stolle 
usando el modelo del contacto entre dos cilindricos perpendiculares con dos 

tensiones máximas alternativas. 
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Profundidad desde el punto de impacto (% del diámetro del cable) 

Figura 56: Modelo matemático de los esfuerzos Internos en el cable del Exp. 4 de C. Stolle. Las tensiones señaladas son las máximas 
calculada, con dos alternativas, y la tensión máxima medida en el experimento, Or. m¿x es el esfuerzo de rotura del cable calculado con la 
tensión de rotura nominal y su área metálica. 
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5.3.3.1 .Discusión de la aplicación del modelo de 
esfuerzos internos debido al contacto al 
experimento 4 de C. Stolle. 


Los resultados, de usar en el modelo los datos experimentales para la 
rotura de dos cuerdas coincide con lo observado en el experimento. En este 
caso la figura 56, muestra que el 58% de la sección del cable está sometido 
a esfuerzos mas allá del límite de rotura, en la experiencia resultó un 
porcentaje similar, dos de las cuerdas se rompieron, 66.66% cuando la 
tensión en los anclajes fue de 171.26 KN. Por otro lado, si se usa el nivel de 
tensión de rotura medido estáticamente, 181 KN., la figura indica que todo el 
cable debería estar roto y así si coincide con lo que se supone en el modelo 

.i) 

ideal. Las tensiones máximas 


considerado para ser JTsa 
que si se usa la tensió 


entonces la predicció 
rota. Pero si se usa una tensió 



Iteffef^Jtivas, de la tabla 37, se han 
atos del modelo idealizado y resulta 
experimentalmente, 171.26 KN., 
í> de la sección del cable estaría 
mayor, 174.9 KN., entonces un 58% 
ra^a y así coincide con la predicción 


de la sección del cabl^-estaría 
hecha al usar los datos Óx|i^^%iplÉ1es. Y si se quiere hacer coincidir la 
predicción del modelo con la rotura de dos trenzas o una rotura del 66.66% 
del cable se debe reemplazar en las ecuaciones una tensión de, 176.5 KN, y 
este sería el nuevo valor de resistencia del cable impactado 
transversalmente predicho con los modelos matemáticos. 


Sin embargo, si se tiene en cuenta que el modelo no predice lo que 
sucede en los experimentos 1 y 2 entonces las predicciones no son 
consistentes. Una explicación para la coincidencia de los resultados es que 
el desplazamiento de las cuerdas y alambres dentro del cable se da en los 
primeros milisegundos. Luego la nueva sección transversal soporta la fuerza 
de impacto sin mayores cambios de su forma. Pero como la velocidad de 
impacto, 10.57 m/s, es hasta 56.6% mayor que las velocidades de impacto 
de los experimentos 1 y 2 en los cuales la tensión máxima fue de 165.5 KN.; 


111 


Esta obra ha sido publicada bajo la licencia Creative Commons Reconocimiento-No Comercial-Compartir bajo la misma licencia 2.5 Perú. 
Para ver una copia de dicha licencia, visite http://creativecommons.Org/licences/by-nc-sa/2.5/pe/ 






BIBLIOTECA DIGITAL - DIRECCIÓN DE SISTEMAS DE INFORMÁTICA Y COMUNICACIÓN 


entonces alcanzar la tensión de rotura en el experimento 4 era inevitable. Lo 
particular de este experimento es que la tensión de rotura fue de 171.26 KN. 
Un valor apenas mayor a las tensiones máximas de los otros 2 
experimentos; pero lo suficiente como para que la predicción del modelo de 
esfuerzos internos coincida con lo observado en el experimento 4. 

En la figura 57 se muestra como quedó el cable al final del experimento 
4 de C. Stolle, solo con una de sus cuerdas sin romper y flexada localmente. 



Figura 57: Situación final del cable 3x7 en el frente del carro que lo 
impactó transversa ¡mente en el experimento 4 de C. Stolle [31]. 
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5.3.4 . Aplicación del modelo de esfuerzos internos 
debido al contacto a la simulación de 
Gospodarczyk. 

Las características de los cables y condiciones de la simulación se dan 
en 5.2.2. El sólido que impacta tiene un frente plano de 1,5m por 2.5m y cae 
en cables de 7.5m de largo, aquí, al igual que en la simulación se 


considerará como rígido, su masa se dividirá entre 126 cables y la 
aceleración se calculará con el modelo de desplazamientos pequeños (Ver 
5.2.2.1) [33], Para la aplicación del modelo solo es necesario el valor final la 
fuerza de contacto y tensión alcanzada. El valor de las variables de entrada 


se 

58 


dan en la tabla 38 y el resultado de los cálculos se presentan en la figura 
y la tabla 39. 


Área metálica del cable. 

1.5000E-04 

m 2 

Diámetro equivalente del cable de acero. 

1.3820E-02 

m 

Longitud del cuerpo rígido que impacta ñ-1. 

2.5000E+00 

m 

Longitud del cuerpo rígido que impacta ñ-2. 

1.5000E+00 

m 

Módulo de Poisson del cable equivalente. 

3.0000E-01 


Módulo elástico del cable. 

1.9500E+11 

Pa 

Tensión máxima calculada: T . ut _ -1 

mod. Mat.Gosp. 

2.7875E+05 

N 

Tensión máxima calculada: T . -2 

mod. Mat.Gosp. 

2.7856E+05 

N 


Tabla 38: Datos de la simulación de Gospodarczyk para el 
modelo de esfuerzos internos debido al contacto. 


Tensión en el cable: Calculada, lado de 1.5m y 2.5m. 

T . =278.56 

calc.l.5m 

T . =278.75 KN 

cale. 2.5m 1X1 

Semiancho del área de contacto b. 

2.1989E-04 

1.7033E-04m 

Presión de contacto P 0 

5.4554E+08 

1.1739E+09 Pa 

Fuerza de contacto máxima según el modelo 
matemático F . _ =M*a„. 1 _ 

max.Gosp. Max.Gosp. 

4.7108E+05 

4.7112E+05 N 

Esfuerzo longitudinal máximo <r L Máx 

1.8583E+09 

1.8571E+09 Pa 

Esfuerzo de Von- Mises superficial ct s M(n 

2.0765E+09 

2.3266E+09 Pa 

Esfuerzo de Von- Mises máximo cr, 

2.0765E+09 

2.3266E+09 Pa 

Esfuerzo de Von- Mises mínimo ct„ 

V.Min 

1.8600E+09 

1.8600E+09 Pa 

Resistencia máxima del cable (T„ 

t^.Max 

1.8600E+09 

1.8600E+09 Pa 

% del material del cable con esfuerzos mayores a 
la resistencia máxima. 

100.00 % 

100.00 % 


Tabla 39: Cálculos con los datos de la simulación de Gospodarczyk usando el 
modelo de esfuerzos internos para los cables bajo el cuerpo que impacta. 
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Figura 58: Modelo matemático de los esfuerzos internos en los cables de la base de segundad para minas simulada por Gospodarczyk. Las 
tensiones señaladas son las máximas calculadas para las cuales el cable se rompe totalmente según el modelo, <r R . Ma -x es el esfuerzo de rotura 

del cable calculado con la tensión de rotura nominal y su área metálica. 
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5.3.4.1 .Discusión de la aplicación del modelo de 
esfuerzos internos debido al contacto a la 
simulación de Gospodarczyk. 

La predicción del modelo indica que la disminución de la capacidad del 
cable por los esfuerzos de contacto no es significativa. El resultado dice que 
hay una disminución pequeña de la capacidad del cable con una resistencia 
en la rotura, medida estáticamente, de 279 KN. a una de 278.75 KN. para los 
cables bajo el lado de 2.5 m y 278.56 KN. para los cables bajo el lado de 1.5 
m del paralelepípedo. Esto coincide con los resultados de Tytko et al, donde 
la resistencia última en la simulación fue mayor, 289 KN. La explicación para 
este resultado están en la construcción del cable y de sus propiedades. El 



viACI WU\i 


de carga o descarga. Es decir, el mode^) da una predicción que coincide 
cualitativamente con lo sinY^^pp^^ el características del cable son muy 
similares a la de una barra sólida. 
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6 . Conclusiones y recomendaciones. 

El impacto transversal en cables de acero ha sido investigado usando 
las leyes generales de la mecánica y construyendo nuevos modelos 
matemáticos. 


El objetivo general de determinar los esfuerzos a lo largo de un cable 
impactado transversalmente en su punto medio se ha logrado parcialmente: 

✓ Se ha logrado construir dos modelos matemáticos: uno de la fuerza 
reacción en los extremos y otro de la cinemática y dinámica del cuerpo 
que impacta. Ambos para ser aplicados cuando el cuerpo que impacta 
lleva a que el cable tenaa desplazamientos grandes o pequeños. 
También se planteópm mod e)p parados esfuerzos internos del cable en 
la zona de impacto. 


Además se ha 
simulados y experiment 



dad de los modelos con datos 
abajos anteriores, que no habían 
lítico de sus problemas. Aunque las 


considerado un planteamie 
predicciones del modelo de los esftórzos internos son no concluyentes 
y no debe usarse tal y como está planteado. 


✓ Como parte de la construcción de los modelos matemáticos se tuvo 
que idear un método, llamado aquí de “movimientos delimitantes”, que 
permite establecer y determinar las condiciones de frontera de una 
cuerda impactada transversalmente: Tensión en sus extremos y 
movimiento de su punto impactado. 


✓ También se encontró que la relación entre un proceso de carga y uno 
de descarga en una cuerda ideal elástica es simétrica respecto al 
tiempo. 
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6.1 . Conclusiones sobre la validación de los 
modelos de los esfuerzos axiales y movimiento 
del cuerpo que impacta. 


Tomando en cuenta la validación con los experimentos de C. Stolle y la 
simulación de Tytko et al, se puede concluir que el modelo matemático de 
las tensiones en los anclajes predice con exactitud, siempre y cuando, el 
módulo elástico del cable se mantenga constante para la etapa de carga y la 
de descarga (si es que llega a darse). También, debido a que el cable usado 
por C. Stolle, y el simulado por Tytko et al, tienen diseños diferentes; se 
puede decir que, con los modelos desarrollados, si es posible predecir el 
movimiento de un cuerpo que impacta en el punto medio a un cable de acero 
y la tensión de este sin ijgjfJórtar $u trp&'ppnstructivo. Siempre y cuando sea 


estabilizado para que 

/ 

principalmente su módu 
como sea posible. 

Y con respecto al modelam 



, como su límite de rotura y, 
lo mas cercanos a una constante 


e la simulación de la plataforma de 


Gospodarczyk. no se puede afirmar qu^fos modelos desarrollados puedan 
modelar el impacto sobre una lí niflfe hecha de cables de acero. Segundo, de 
lo expuesto, queda claro que la predicción de la velocidad se aproxima al 
resultado de la simulación debido a las condiciones particulares de la misma 
y las propiedades del tipo de cable simulado. Aún así, los modelos si se 
pueden ampliar para modelar estructuras más complejas, aunque es 
necesario contar con más datos experimentales. 


6.2 . Conclusiones sobre la aplicación del modelo de 
los esfuerzos internos debido al contacto. 

Ya que las predicciones no coinciden con los resultados de los 
experimentos 1 y 2 de C. Stolle, entonces el modelo desarrollado no es 
adecuado para estos casos. Aunque tampoco se puede decir que no sirva 
para otros casos y cables por que falta analizar la influencia de los 
parámetros seleccionados para representar las propiedades del cable. Aquí 
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se han escogido como se indica en el trabajo de la referencia [72], Pero está 
claro que en un cable de acero, en general, el módulo de Poisson debería 
ser mayor a 0.3 y el módulo elástico debería ser medido aplicando una 
fuerza transversal ya que, respecto a esta propiedad, el comportamiento de 
los cables es anisotrópico. Una mejor aproximación sería considerarlo 
ortotrópico como el material de sus alambres. 


Pero de acuerdo a la predicción con los datos del experimento 4 de C. 
Stolle y con los de Gospodarczyk, si se puede decir que el modelo de los 
esfuerzos internos por contacto durante el impacto transversal en cables 
puede arrojar predicciones aproximadas siempre y cuando la sección 
transversal del cable no sufr§ .c6ftíbtes.significativos durante el impacto. 


Aún así lo anterior 
esfuerzos internos funci 
sin alma sintética, 
desarrollar o mejorar 



tivo de que el modelo de los 
una sola cuerda o de más de 3 
ás pruebas experimentales y 
pueda predecir el comportamiento 


de una cable sometido^ aplastamiento por una fuerza transversal a su eje 


longitudinal por que con 'eí«fc|S^ri[plj^ao aquí no se puede llegar ha hacer 
alguna predicción concluyente sobre el estado de esfuerzos en la zona de 
impacto de un cable sometido a impacto transversal. 
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6.3 . Resumen de resultados 

Se ha hecho lo siguiente: 


✓ 


✓ 


Con los modelos se ha podido calcular la fuerza de reacción en los 
extremos o anclajes del cable y deducir que se aproxima a la tensión 
en el punto medio que es donde se produce el impacto de un sólido 
modelado como punto material, además también se ha modelado, con 
buena aproximación, y analizado la dinámica del cuerpo que impacta. 
Y también se encontró que mientras un cable tenga anclajes rígido y 
que su módulo elástico sea más cercano a una constante (cable 
estabilizado) entonces los modelos matemáticos se aproximan mejor. 


TjV-H-- 

También se ha determinado que modelar los esfuerzos internos que 

un cable debido al contacto con un 


surgen en la sección 

/ 

sólido que lo impa 
cierta geometría 
usando la teoría de esfuerz 



concluyente debido a que no tiepe 

. '(/Ir.. -sááÉF* 


el contacto durante el impacto de 
a que idealiza al cable de acero, 
contacto de Hertz no llega a ser 
en cuenta los deslizamientos 


relativos de los alam 


que todo tipo de cable presenta al 


ser sometido a una fuerza transversal. 


✓ La importancia del método presentado también radica en que puede 
ser usado junto al “método de características” para hallar la geometría 
y esfuerzos en todos los puntos de un cable impactado luego de 
solucionar el sistema de ecuaciones diferenciales parciales parabólicas 
de segundo grado que representan el movimiento general de una 
cuerda en un plano. 

✓ Además, se ha determinado que la ecuación de la tensión durante la 
carga (aumento de la tensión) de una cuerda impactada 
longitudinalmente por una masa puntual, es una función simétrica de la 
ecuación que modela a una cuerda idéntica o la misma cuerda con el 
mismo módulo elástico durante su proceso de descarga. 
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En suma, a parte de dar un método y modelos matemáticos se ha 
logrado una herramienta ingenieril para poder diseñar barreras de protección 
y así dar más y mejores soluciones a los accidentes de tránsito, a los 
accidentes por caída de objetos en minas y las caídas de rocas. 


6.4 . Sugerencias para trabajos futuros. 

✓ En cualquier otro desarrollo o mejora que se haga de los modelos aquí 
planteados debe de tener como base la filosofía del cambio, de otra 
manera lo desarrollado será solo un algoritmo más. 


✓ 


✓ 


De la misma forma que en este trabajo se ha desarrollado y usado el 


método de “movimie^ÍPSIbí^limitantes” para una cuerda ideal 

- su punto medio por un punto material; 

•£ 

método para cuando el impacto es 


impactada transversal me 
también es posible ,- 


oblicuo o se da e 


estudiar los trabajos 


manteniendo velocidad del 



s el medio. Aunque antes se debe 
ara cuando el impacto se da 
que impacta constante, estos se 


pueden encontrar solamente en la referencia [361. 


Otro punto importante es tener en cuenta que el factor de seguridad a 
la hora de diseñar un estructura con cables sometidos a impacto 
transversal u oblicuo, variará de acuerdo a la capacidad de los 
absorbentes de energía que se use en los extremos. 


✓ Como uno de los supuestos bases de los modelos desarrollados fue 
que los anclajes sean perfectamente rígidos, esto supone que el 
modelo no se puede aplicar en los casos donde hay absorbentes de 
energía; por lo que si se quiere incluir estos elementos es necesario 
modificar la ecuación (22) haciendo un nuevo análisis geométrico 
como el de la figura 17 e incluir una reflexión que no duplica el nivel de 
tensión que llega a los anclajes. 
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✓ Otra sugerencia es tratar de mantener la respuesta de los cables en la 
zona de esfuerzos elásticos de manera que los modelos matemáticos 
desarrollados aquí se aproximen lo más posible a la realidad. Es cierto 
que se puede usar el software LS-DYNA u otro parecido pero los 
trabajos hechos con el software no revelan los códigos ni la totalidad 
de los parámetros que han usado. Y si se quiere hacer experimentos 
se tiene que tener en cuenta que el estándar ASTM A741 -11 (2016) 
sugiere otros 3 tipos de cables pero estos tipos son denominaciones 
genéricas de varias construcciones que varían según el fabricante lo 
que hace difícil usar solo experimentos para desarrollar diseños de 


barreras de pro tíir ' r ' ifSn 
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8 . Anexos. 

8.1 . Anexo 1. 

Se trata de modelar el impacto entre un punto material y una cuerda 
desde el Instante en que se contactan. Para esto se presenta la figura 59. 


M i voo 



Modelo aproximado 
con perfil triangular 


Idealización del cable 


▼y 


L 


Figura 59: Modelo del impacto transversal de un cable aproximado con un 
perfil triangular y desde el Instante de contacto.M, es la masa del cuerpo, V 00 , 
es su velocidad inicial, V, es su velocidad actual, O, es el origen de 
coordenadas de ejes x, y, Y 0 , desplazamiento del punto medio de las dos 
cuerdas. L, es la longitud del cable, F, es la fuerza de reacción. Cu, es la 
velocidad del frente de la onda de esfuerzos longtudinal. X* es la distancia 
recorrida por el frente de onda transversal. 

Aquí se representa ^'Gl3í^§y$£8éTideal en los primeros instantes del 

movimiento que adquiere por el impacto del punto material de masa M y 
velocidad inicial, V 00 . La onda transversal aún no ha alcanzado los 
extremos y la onda longitudinal tampoco. La línea curva que es la 
idealización del cable mientras que el perfil triangular es su aproximación. 
Entonces se puede formular las relaciones matemáticas son las siguientes: 


F—2 Tsin a 


(85) 



( 86 ) 



(87) 
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F=Usina (88) 

El extremo x* del triángulo se halla usando la velocidad de la onda 

transversal: 


x ‘ = i* tóüW 


T ,» 


A cable P 


(89) 


En la ecuación (89), B , es una constante y T es la tensión de la 
cuerda en el punto de impacto y viene a ser: 


E — A cable E 


Luego reemplazan^ 
distancia x* : 


Y reemplazando 
del punto material vendría 
Fowler: 



(90) 


obtiene una aproximación para la 


y Y (91) 

(87), la ecuación del movimiento 
, una ecuación del tipo Endem- 





y 0 _ 2 5/4 Á 7l4 pE 314 


d t t 


M 




(92) 


Ahora, si se considera la suposición de Yu. A. Dem’yanov, el 
desplazamiento del punto material tiene que ser (como la ecuación (11)) [67]: 


y(So=0,t) = y 0 (t)=B 1 t n ;B 1 >0 (93) 

Si esta relación se reemplaza en la ecuación (92), se puede dar una 

solución cuando n=2 y B 1 =(b 13 /2) 2 . Hasta acá todo parece lógico pero si 
se examina el coeficiente B 13 , la masa está dividiendo a las otras cifras 
indicando que mientras más masa tenga el punto material menor aceleración 
o menor será la fuerza F sobre el punto material durante los primeros 
instantes del impacto. Esto no tiene sentido físico por que mientras más 
masa tenga el objeto que impacta mayor será su cantidad de movimiento lo 
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que debe resultar en que la cuerda, con una masa mucho menor, ideal y 
hecha de una continuidad de puntos materiales, acelere mucho más o con 
mayor fuerza mientras la masa que impacta es mayor. Coloquialmente 
hablando: mientras el objeto que choca tenga más masa, más fácil será que 
lo lleve de un plumazo. 

Entonces, la otra alternativa es aplicar directamente una ley general del 
movimiento de Euler, la ley general del momentum lineal para una cuerpo 
unidimensional en un espacio bi-dimensional, como la cuerda. Luego la 
ecuación del movimiento es: 


s~r^-=MV 


Pero acá de nuévo 
cuerda que contacta 
primer término a la izquierd 



(94) 


ntét Si la velocidad del punto de la 
s inicialmente cero entonces el 
igual de la ecuación (94) es cero, 
a onda aún no se ha propagado. Sin 


Pero también lo será la integral p 
embargo, el lado derecho no puedef^ser cero porque la cantidad de 


' r Ao ifis 1 

movimiento inicial en el sistema no es nula. 


Entonces, ¿Qué pasa?. Y la respuesta se encuentra en los siguientes 
hechos: 


1. El suponer que en el impacto con una cuerda ideal unidimensional 
formada por una continuidad de puntos materiales (de masa 
infinitesimal) que no tiene resistencia a la flexión provoca que se 
forme una discontinuidad en la tangente del punto donde el punto 
material de masa, M , impacta. 

2. La ley general del movimiento de Euler, para la conservación de la 
cantidad de movimiento, se puede aplicar a cualquier cuerpo 
unidimensional en un espacio bidimensional (o tridimensional) 
siempre y cuando su geometría sea una curva continua de clase C3. 
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Los dos punto anteriores explican el porque no funciona la ley de Euler 
pero no explican el primer desarrollo, de las ecuaciones (85) a (93), y la 
respuesta está en que ahí se usa la tercera ley de Newton, que solo es 
aplicable a puntos, literalmente solo a puntos materiales con masa discreta y 
no a la continuidad de estos como la cuerda. Y es que además, Newton, 
nunca formuló esta ley. En realidad enunció tres leyes que luego se 
trasladaron a fórmulas matemáticas por sus seguidores y solo lo aplicaron a 
masas puntuales. Fue en realidad, L. Euler, quien dio las leyes de 
movimiento generales y que se pueden aplicar a cualquier cuerpo ya sea 
sólido rígido o deformable [78], [79], En la cita que sigue se explica esto: 


"Every body perseverp§\iW ií ít§*%ta}e of rest or of uniform motion straight 

A) 

ahead, unless it be compellejy^CafajTge that State by forces impressed upon 


it." In the generalityíína 
always valid, for a bod 
expressed in terms of a 



erTP mechanics, this axiom is not 
internal or external constraints not 
. For example, a rigid body subject 
xis through its center of mass; its 

mass 


to no applied forcé spins about 

parts, which also are bodies. move in such a way that their centers of 
describe circles about that axis. Newton himself did not specify any 
mathematical properties of bodies or forces, and so his intentions must be 
inferred by the reader, and in the course of time different readers have read 
different meanings into his words.” [79] 


El desarrollo anterior también quiere decir que un modelo de cuerpos 
unidimensionales (puntos y cuerdas) no sirve para modelar el impacto entre 
cuerpos ya que no se puede tener en cuenta la velocidad relativa que estos 
presentan durante el intervalo de tiempo entre el primer contacto y entre el 
instante en que adquieren una velocidad común. Además, este último valor 
se da entre los puntos de la parte de sus superficies en contacto y no 
necesariamente entre los centros de gravedad de los cuerpos. Por esto, el 
modelo usual de dos esferas rígidas que impactan, que al final se simplifica 
al de dos puntos que chocan, no es adecuado para captar todos los 
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fenómenos tribológicos que suceden en la superficie de contacto de 
procesos reales como en un molino de bolas. También, la ley de la 
conservación de la cantidad de movimiento aplicada a cuerpos puntuales 
puede aproximar lo que sucede antes y después del intervalo de tiempo 
donde se da el impacto o contacto entre dos esferas pero no da detalles 
sobre lo que pasó durante ese intervalo y así no se sabe cuales fueron los 
verdaderos movimientos de los puntos de los cuerpos. Por lo que tampoco 
se puede determinar los esfuerzos dentro del cuerpo durante la colisión. 


8.2 . Anexo 2. 



En la figura 60 se presentan los resultados de otros 3 modelos 
numéricos desarrollado&ppor varios años para cables impactados y que 
fueron usados por C. Jftolle para contrastadla efectividad se su simulación. 
Aquí se presenta el modero/matemático de la tensión en los anclajes 
aplicado al Exp. 1 de C. Stolle. También, para compararlo con estos otros 3 
modelos numéricos. 


135 


Esta obra ha sido publicada bajo la licencia Creative Commons Reconocimiento-No Comercial-Compartir bajo la misma licencia 2.5 Perú. 
Para ver una copia de dicha licencia, visite http://creativecommons.Org/licences/by-nc-sa/2.5/pe/ 







0.05 0.1 0.15 0.2 0.25 0.3 0.35 


BIBLIOTECA DIGITAL - DIRECCIÓN DE SISTEMAS DE INFORMÁTICA Y COMUNICACIÓN 


O O O O O O O 

olooloolooo 

^■rororsirsiT-T-uno 



ooooooooo 

omomoinoLn 

'í m ro c\j rvj r- t- 


(N>l) ± safepue so| ua upisuaj. 


Esta obra ha sido publicada bajo la licencia Creative Commons Reconocimiento-No Comercial-Compartir bajo la misma licencia 2.5 Perú. 
Para ver una copia de dicha licencia, visite http://creativecommons.Org/licences/by-nc-sa/2.5/pe/ 


Tiempo (S) 

Figura 60: Comparación del Modelo Matemático de la tensión T(KN) en los anclajes del cable ¡mpactado transversal mente en el experimento 
1 de C. Stolle con la simulación de C. Stolle y otros modelos numéricos que usó para comparar su simulación [31]. 




